Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	55,591	0	T	A	43.6%	24.6 / 11.7	39	N36K (AAT→AAA)	BCHEFKBB_00050	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/3);  new base A (17/0);  total (36/3)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCGCAAAGGTAACGCGCCCCTCTGACGGGGTGCAATGTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑TTGCCGCCTGCGCGGCAATGGCGGG‑CCGGGGAGGCGGCATTCTCCAAGTTGCGGCAACCTTTCCCAACAACT‑TAACCGCCTAAATACAAGCCTTGCGGCTTGTTGCCCTCTCTCCAGCTCTCTCCCACAGGGAGAGAGGACGGGTCGGCTGTTGGGGTTGAGGTT  >  Exported/55450‑55737                                                                                                                                                 |                                                                                                                                                     tCGCAAAGGTAACGCGCCCCTCTGACGGGGTGCAATGTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑AATAGCC‑CCAGCGCGGCCATGGCggggcg                                                                                                                                             >  1:453764/1‑148 (MQ=255)        ggTAACGCGCCCCTCTGACGGGGTGCAATGTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCATCCGGTTGCCCGGCCCG‑CC‑A‑ATGCCGCCTGCGCGGCAATGGCGGG‑CCggggggg                                                                                                                                      >  2:1038540/1‑147 (MQ=255)        ggTAACGCGCCCCTCTGACGGGGTGCAATGTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑TTGCCGCCAGCGCGGCCATGGCGGG‑CCggggggg                                                                                                                                      >  2:299349/1‑147 (MQ=255)              gcgcCCCTCTGACGGGGTGCAATGTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑ATGCCGCCGGCGCGGCCAAGGCGGG‑CCGGGGAggcggcctt                                                                                                                               >  2:1000975/1‑148 (MQ=255)                   ccTCTGACGGGGTGCAATGTATAGGGTTGCAGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑ATGCCGCCAGCGCGGCAATTGCGGG‑CCGGGGAgggcgcaatagcca                                                                                                                          >  2:839362/1‑139 (MQ=255)                   catctGACGGGGAGCAATGAATAGGGTTGACGCATCATTCTCCTTGCGCGCTTATCAATGGTACAGCGGTTTATGGCAATCCGCAACCGATTGCCCGGCCCG‑CC‑A‑ATGACTCCTTCGCGGCAATGGCGGG‑CCGGGGGGGCGCCAGTCTCCa                                                                                                                          >  2:31504/3‑151 (MQ=255)                    ctctGACGGGGTGCAATGTATAGGGTTGCGGGGTCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑C‑CTGGCGCCGGCGCGGCAAAGGGGGG‑GCGGGGGggcggccgacaacca                                                                                                                         >  1:205687/1‑142 (MQ=255)                              gTGCAATGTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CCCA‑TT‑CCGCCGGCGCGGCAAAGGCGGG‑GCCGGGAGGCGGCAATCTCCaacgttcggcca                                                                                                               >  2:1289721/1‑141 (MQ=255)                              gTGCAATGTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑TTTCCCCCAGCGCGGCAAAGGCGGG‑CCCGGGAGGGGGGATTCTACAAGTTTCGGCaa                                                                                                               >  2:43250/1‑151 (MQ=255)                                     gTATAGGGTTGCGGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CCCA‑TT‑CCGCCGACGCCGCAATGGCGGG‑CCGGGGATGCGGCCATCAACAAATTGCGGCAACCATTc                                                                                                         >  2:884762/1‑150 (MQ=255)                                     gTATAGGGTTGCCGCGTCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑AATTGCC‑CCAGGGCCGCAATTGCGGG‑GCGGGGGGGCGGCAATATACAAAATGCGGCAAccaata                                                                                                         >  2:884756/1‑146 (MQ=255)                                     gTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCGGCA‑A‑TTGCC‑CCCGCGCGGCCAAGGCGGGGCCGGGG‑GGCGGCCGTGTCCAAGTTTAGGCAAccaata                                                                                                         >  1:1193438/1‑146 (MQ=255)                                     gTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑GC‑A‑ATTCCGCGAGCGCGGCAATGGCGGG‑CCCGGGAGGCGGGATTAAACAAGGTGCGGCAACCATTc                                                                                                         >  2:953054/1‑150 (MQ=255)                                     gTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑TTTCCGCCTGCGCGGCAATGGGGGG‑CCGGGGAGGCGGCAATCAGCAATTTGCGGCAACCATTc                                                                                                         >  2:159259/1‑150 (MQ=255)                                     gTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑TTGCCGCCTGCGCGGCAAAGGCGGG‑CCGGGGAGGCGGCATTATCCAAGTTGCGGCCACCGTTc                                                                                                         >  2:643428/1‑150 (MQ=255)                                     gTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑TTGCCGCCAGCGCGGCAAATGCGGG‑CCCGGGCGGCGGCCATATCCAACTTGAGGCAAccgtta                                                                                                         >  2:373375/1‑146 (MQ=255)                                     gTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑ATTCCGCCAGCGCGGCCCAGGGGGG‑CCGGGGAGGCGGCAATCTCCCACTTGCGGCAAccatta                                                                                                         >  2:800139/1‑146 (MQ=255)                                     gTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑ATGGCGCCAGGGCGGCAAAGGCGGG‑CCGGGGGGGCGGCAATAAACAAGTTTCGGCAACCGTTc                                                                                                         >  2:1031376/1‑150 (MQ=255)                                     gTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑ATGCCGCCAGCGCCGCCAAGGCGGG‑CCCGGGGGGGGGCAGTCTCCAAATTGCGGCAAccgtta                                                                                                         >  2:768435/1‑146 (MQ=255)                                     gTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACAGTTTGCCCGGCCCG‑GC‑C‑TTGGCGCCTGCGCGGCAAATGCGGG‑CCGGGGAGGCGGCATTATACAAGTTGCGGCAAccgtta                                                                                                         >  2:358155/1‑146 (MQ=255)                                     gTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGGTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑ATGCCCCCTGCGCCGCAAAGGCGGG‑CCCGGGAGGCGGCCAACACCAAGTTGCGGCAACCATTc                                                                                                         >  2:1273226/1‑150 (MQ=255)                                     gTATAGGGTTGCCGCATCATTCGCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑C‑TTGCCCCCTGCGCGGCAAAGGCGGG‑CCCGGGCGGCGGCATTCTCCAAGTTGCGGCAACCttta                                                                                                         >  2:1147017/1‑149 (MQ=255)                                          gggTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑GC‑A‑TTGCCGCGGGCGCGGCAATGGCGGG‑GCGGGGAGGCGGCATTCTCCAAGGTGCGGCAAACTTTGCcaaca                                                                                                   >  2:35271/1‑151 (MQ=255)                                                  gCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑ATGCCCCCAGCGCCGCCAAGGCGGG‑CCCGGGAGGCGGCCTTATCCAAGTTGCGGCAACCTTTCCCAACAACA‑Taaa                                                                                           >  2:1350481/1‑149 (MQ=37)                                                     tcatTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑ATGACGCCAGCGCGGCCAAGGCGGG‑CCCGGGAGGCCGCAATATACAAGATTCGGCAACCGTTACcaaaaaaataaaacc                                                                                         >  2:90031/1‑138 (MQ=255)                                                            ccTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑ATTCCGCCAGCGCGGCAAATGCGGG‑CCGGGGAGGCGGCATTAACCAATTTTCGGCAACCTATACCAAAAACT‑AAAACGccccaataa                                                                                >  2:1283569/1‑144 (MQ=255)                                                                   cgcTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑TTGCCCACTGAGCGGCGATTGCCGG‑CCGGGGAGGCGGCCTTCTCCCAGGTGCGGCAACCGTTGCCAACCACT‑TAAACGccgtaaaaaaacccgt                                                                         >  2:1100428/1‑137 (MQ=255)                                                                          ccATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑TTTCCGCCTGCGCGGCAAAGGCGGG‑CCGGGGAGGCGGCATTCTCCAAGTTGCGGCAACCTTTCCCAACAACA‑TAACCGCCTAAATTCaaaccatgaggggg                                                                  >  2:848158/1‑139 (MQ=25)                                                                                                                       cg‑CC‑A‑TTGCCGCCTGCGCGGCAATGGCGGG‑CCGGGGAGGCGGCATTCTCCAAGTTGCGGCAACCTTTCCCAACACCT‑TAACCGCCTAAATACAAGCCGTGCGGCTTGTTGCCCTCTCTCCAGCTCTCTCCCACAGGGAGAGAGGACGGGt                      >  1:733593/1‑150 (MQ=255)                                                                                                                          cc‑A‑TTGCCGCCTGCGCGGCAATGGCGGG‑CCGGGGAGGCGGCATTCTCCAAGTTGCGGCAACCTTTCCCAACAACT‑TAACCGCCTAAATACAAGCCTTGCGGCTTGTTGCCCTCTCTCCAGCTCTCTCCCACAGGGAGAGAGGACGGGTCgg                   >  1:766355/1‑151 (MQ=255)                                                                                                                          cc‑A‑TTGCCGCCTGCGCGGCAATGGCGGG‑CCGGGGAGGCGGCATTCTCCAAGTTGCGGCAACCTTTCCCAACAACT‑TAACCGCCTAAATACAAGCCTTGCGGCTTGTTGCCCTCTCTCCAGCTCTCTCCCACAGGGAGAGAGGACGGGTCgg                   >  1:464109/1‑151 (MQ=255)                                                                                                                             a‑TTGCCGCCTGCGCGGCAATGGCGGG‑CCGGGGAGGCGGCATTCTCCAAGtt                                                                                                                      >  1:174732/1‑51 (MQ=255)                                                                                                                             a‑TTGCCGCCTGCGCGGCAATGGCGGG‑CCGGGGAGGCGGCATTCTCCAAGtt                                                                                                                      <  2:174682/51‑1 (MQ=255)                                                                                                                             a‑TTGCCGCCTGCGCGGCAATGGCGGG‑CCGGGGAGGCGGCATTCTCCAAGTTGCGGCAACCTTTCCCAACAACT‑TAACCGCCTAAATACAAGCCTTGCGGCTTGTTGCCCTCTCTCCAGCTCTCTCCCACAGGGAGAGAGGACGGGTCgg                   >  1:474267/1‑149 (MQ=255)                                                                                                                             a‑TTGCCGCCTGCGCGGCAATGGCGGG‑CCGGGGAGGCGGCATTCTCCAAGTTGCGGCAACCTTTCCCAACA‑CTCTAACCGCCCCAATACAAGCCTTGCGGCTTGTTGCCCTCTCTCCAGctc                                               <  2:874187/121‑1 (MQ=11)                                                                                                                             a‑TTGCCGCCTGCGCGGCAATGGCGGG‑CCGGGGAGGCGGCATTCTCCAAGTTGCGGCAACCTTTCCCAACA‑CTCTAACCGCCCCAATACAAGCCTTGCGGCTTGTTGCCCTCTCTCCAGctc                                               >  1:874423/1‑121 (MQ=11)                                                                                                                                           cgGCACTGGCGGG‑CCGGGGAGGCGGCATTCTCCAAGTTGCGGCAACCTTTCCCAACAACG‑TAACCGCCCCAATACAAGCCTTGCGGCTTGTTGCCCTCTCTCCAGCTCTCTCCCACAGGGAGAGAGGACGGGTCGGCTGTTGGGGTTGAgg    <  2:742635/151‑1 (MQ=11)                                                                                                                                           cgGCAATGGCGGG‑CCGGGGAGGCGGCATTCTCCAAGTTGCGGCAACCTTTCCCAACACCG‑TAACCGCCTAAATACAAGCCTTGCGGCTTGTTGCCCTCTCTCCAGCTCTCTCCCACTGGGAGAGAGGACGGGTCGGCTGTTGGGGTTGAgg    >  1:742840/1‑151 (MQ=255)                                                                                                                                             gCAATGGCGGG‑CCGGGGAGGCGGCATTCTCCAAGTTGCGGCAACCTTTCCCAACAACT‑TAACCGCCTAAATACAAGCCTTGCGGCTTGTTGCCCTCTCTCCAGCTCTCTCCCACAGGGAGAGAGGACGGGTCGGCTGTTGGGGTTGAGGtt  >  2:610682/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                     TCGCAAAGGTAACGCGCCCCTCTGACGGGGTGCAATGTATAGGGTTGCCGCATCATTCTCCTTGCGCGCTTATCCATGGTACAGCGGTTTTTGGCAATCCGCAACCGTTTGCCCGGCCCG‑CC‑A‑TTGCCGCCTGCGCGGCAATGGCGGG‑CCGGGGAGGCGGCATTCTCCAAGTTGCGGCAACCTTTCCCAACAACT‑TAACCGCCTAAATACAAGCCTTGCGGCTTGTTGCCCTCTCTCCAGCTCTCTCCCACAGGGAGAGAGGACGGGTCGGCTGTTGGGGTTGAGGTT  >  Exported/55450‑55737

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑