Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,256,061	0	C	T	30.8%	53.7 / 10.0	39	E136K (GAA→AAA)	plsB	glycerol‑3‑phosphate 1‑O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (7/20);  new base T (12/0);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 1.29e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TATTAAACAGATCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGCCCCGGCGCGCGACCAAACA‑TCACCGACA‑C‑TGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTG‑TAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTC  >  NC_000913/4255812‑4256305                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         tattaAACAGATCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGACAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTTCCCTCCCCC‑TTTTT‑c                                                                                                                                                                                                                                                          >  2:429383/1‑251 (MQ=255)     aaaCAGATCCTGACGCGCAGGCAAGCGGGGGCCTACGGCTGCCTGACGTGGACGGTCAAAGTGCAAACGCGCCACGCGCGCCAGTATCGGAGCTATAGTTTTATCCGTGCCTTGTTCATCCGCCATACGGCGCAGCTAAACTGGCGGCGAGAAACGCACAAAACTTTCTCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGGTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGcc                                                                                                                                                                                                                                                     <  1:601472/251‑1 (MQ=255)     aaaCAGATCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGcc                                                                                                                                                                                                                                                     <  1:399841/251‑1 (MQ=255)             ccTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑GTCCCTTCGCCCTTTTT‑C‑CGCCCCCgcgcc                                                                                                                                                                                                                                             >  1:376888/1‑250 (MQ=255)                          cAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTTACCTTGCC‑TTTTT‑TTCGCCCCcgcgcggggcccaaaatc                                                                                                                                                                                                                                 >  1:284699/1‑238 (MQ=255)                              cGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑TCCGCCCCCGCGCGCGCCCCAACA‑TCAcccc                                                                                                                                                                                                                           >  2:178571/1‑249 (MQ=255)                              cGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑GTCAATTCTCC‑TTTTTTCA‑GCCCCCGCGCCCGACCCAACA‑TCAccc                                                                                                                                                                                                                            >  1:115828/1‑249 (MQ=255)                                 gggCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑GTTACTTCGCC‑CTTTT‑TTCGCCCCCGCCCGCCGCCAAACA‑ACAcccccac                                                                                                                                                                                                                        >  2:153467/1‑246 (MQ=255)                                            ggccgaaCGTTGAGGGGAAAAGTGCATACGCGCCACGCGCGCCATTTTCTGAGCGATAGTTTTACCCGTGCCGTGTTCATCCGCCATACGGCGCAGCTAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGTCCCGGCGCGCGACCAAACA‑TCACCGACA‑C‑TGGCACCATCt                                                                                                                                                                                                           <  1:41553/245‑1 (MQ=255)                                             gCCAGACGTTGACGGGCAAAGTGCATCCGCGCCACGCGCGCCATTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTTTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGCCCCGGCGCGCGACCAAACA‑TCACCGACA‑C‑TGGCACCATCTg                                                                                                                                                                                                          <  1:124953/251‑1 (MQ=255)                                             gCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑CGGACTTCGCC‑CTTTT‑CTTACCCCGGCGCGCGGCCACACA‑TCaacccccctcggtcccctct                                                                                                                                                                                                            >  2:498557/1‑232 (MQ=255)                                                        aCGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTCCGGAGCGATAGTTTTATCCGTGCCGGGTTCATCCGCCATCCGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGCCCCGGCGCGCGACCAAACA‑TCACCGACA‑C‑TGGCACCATCTGCACATCCAGAt                                                                                                                                                                                               <  1:284802/251‑1 (MQ=255)                                                                      tACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑GTCACTTCGCC‑TTTTT‑TACCCCCCCGCGCGCGGCCCAACA‑TCACCCACA‑C‑TCGCACCCTCTGCCCATCCCGATTTTGGGTGCTAccg                                                                                                                                                                                 >  1:18689/1‑249 (MQ=255)                                                                                      cAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACGGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGCCCCGGCGCGCGACCAAACA‑TCACCGACA‑C‑TGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGt                                                                                                                                                                 <  1:521744/251‑1 (MQ=255)                                                                                         tttCTGAGCGATAGTTTTAGCCGGGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGCCCCGGCGCGCGACCAAACA‑TCACCGACA‑C‑TGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCg                                                                                                                                                               <  1:125914/250‑1 (MQ=255)                                                                                             ggAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGCCCCGGCGCGCGACCAAACA‑TCACCGACA‑C‑TGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGaa                                                                                                                                                          <  1:429383/250‑1 (MQ=255)                                                                                               aGCGATGGTTTTATCCGTGCCGTGTTCATCCGCCATGCGGCGCAGCGAAACTGACGGCGAGAACCGCACAAAACTTTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGCCCCGGCGCGCGACCAAACA‑TCACCGACA‑C‑TGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAGTCGAGATAGTCGTGGAACa                                                                                                                                                        <  1:498557/251‑1 (MQ=255)                                                                                                 atagagttttacccgggcctgtttaaTCCGCCATACGGCGCAGCGAAACTGACGGAGAGAAACGAACAAAACTGTGGCGCCCGGGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGCCCCGGCGCGCGACCAAACA‑TCACCGACA‑C‑TGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGc                                                                                                                                                      <  2:149150/226‑1 (MQ=255)                                                                                                  gATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGCAACTGACGGCGAGAACCGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑GTCACTTCGCC‑CTTTT‑TACGCCCCCCCGCGCGGCCCAACA‑TCACCCACA‑C‑TGGCCCCCTCTGCACCTCCCGATTTTGGTTGCCACCGTGCAAAACGCGGTAGTCGTGGAACAGCt                                                                                                                                                     >  1:625621/1‑251 (MQ=255)                                                                                                                 gccgtgTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCTAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑TACCCCCCCGCCCGCGACCCAACA‑CCCACGACAAC‑TGCCCCCA‑CTTCACACCAACAATTAGGTTTCTACCGTGCCAATCGAGAAAGTCCGGGAACAACTTAATAGACTcttctt                                                                                                                                      >  1:480604/1‑251 (MQ=255)                                                                                                                         caaccGCCATCCGGCGCAGCGAAACTGCCGGCGAGAAACCCACAAACCGGTCGCGCCCGAGCCACAGTACAGCGAAAAATTTCTGTCCGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGCCCCGGCGCGCGACCAAACA‑TCACCGACA‑C‑TGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCgtg                                                                                                                              <  2:492849/248‑1 (MQ=255)                                                                                                                                       cgcAGCGAAACGGACGGCGAGAAACGCACAAACCTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGCCCCGGCGCGCGACCAAACA‑TCACCGACA‑C‑TGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTG‑TAATAGGTGAacac                                                                                                               <  2:18689/251‑1 (MQ=255)                                                                                                                                         cAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGCCCCGGCGCGCGACCAAACA‑TCACCGACA‑C‑TGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTG‑TAATAGGTGAacacac                                                                                                             <  1:298046/251‑1 (MQ=255)                                                                                                                                            cGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGCCCCGGCGCGCGACCAAACA‑TCACCGaca‑c                                                                                                                                                                                                                       <  1:96290/144‑1 (MQ=255)                                                                                                                                              aaaCTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑ATCCCTTCGCC‑TTTTT‑TACGCCCCCGCGCGCGACCAAACA‑TAACCGACA‑C‑TTGCACCCTCTGCACATCCCGATTTTGGGTGCTACGGTGCCAAAAGCGATAGTCGTGGGACAGCTTTATATACTCTTCTTTTCGCGTG‑TAAAAGGTGAAAACACGCGGc                                                                                                        >  1:122735/1‑251 (MQ=255)                                                                                                                                                 cTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTTGCC‑TTTTT‑CACGCCCCCGCGCGCGACCAAACA‑TCACCGACA‑C‑TGGCACCCTCTGCACATCCAGATTTGGGTTGCTACGGGGCCAATCGAGATAGTCGTGGGACAGCTTAATAGACTCTTCTTTTGGCGTG‑TAATAGGTGAAAAAAAGCGGcccc                                                                                                     >  1:215505/1‑250 (MQ=255)                                                                                                                                                   gcgggcgagaaccccaaaaacgtGTCGCGCCCGAGCCGCAGTCCAGCGAAAAATTTCTGTACGCCGGTAAGCATACGCAGCGGCGGG‑TTCACTGCGCC‑TTTTT‑CACGCCCGGGCGCGCGACCAAACA‑TCCCCGACA‑C‑TGGCACCATCTGCACATCCAGAGTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACATCGGAAGAGACTCTGCGTTGGGCGTG‑TAATAGGTGCACACACTCGGCccg                                                                                                     <  2:286532/227‑1 (MQ=255)                                                                                                                                                          agaAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑TACGCCCCCGCGCGCGACCCAACA‑ACACCCACA‑C‑TGGCACCATCTGCACATCCAGATTTTGGGTGCTACCGGGCAAATCGAGAGAGTAGTGGAACAGCTTAATAGACTCTTCTTTTTGCGTG‑TAATAAGTGAAACCACACGGCCCCCCCtgaatg                                                                                            >  1:159079/1‑251 (MQ=255)                                                                                                                                                                gCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTTACTTTGCC‑TTTTT‑TACGCCCCGGCGCGCGACCCAACA‑ACACCGCCA‑C‑TCGCCCCCTCCTCCCATCCCGATTTTGGGTGCTCCGGGGCCAAACGAGATAGTCGGGGGACAACTTTATAAAACCTTTTTTCGGCGTG‑TAATAGGGGAAAACACCCGGCCCCCCGgggatgaaaac                                                                                       >  1:474343/1‑247 (MQ=255)                                                                                                                                                                   aaaaacggccgcgACCGAGCCACAGTACAGCGAAAAATTTCTGTCCGCCGTTAAGCATCCGCAGGGGCGGG‑TTCACTTCGCC‑TTTTT‑CCCCCCCCGGCGCGCGACCAAACA‑TCCCCGACA‑C‑TGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGCCTCTTCTTTCGGCGTG‑TAATAGGTGAACACCCGCGGCCCGCCGTGAATGAACACATAg                                                                                   <  1:178571/242‑1 (MQ=255)                                                                                                                                                                    aaaaCTGTCGCGACCGAGCCACAGTACAGCAAAAATTTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGCCCCGGCGCGCGACCAAAAC‑TCACCGACA‑C‑TGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGAAACAGCTTAATAGACTCTTCTTTCGGCGTG‑TAATAGGGGAACACACGCGGCCCGCCGTGAATGAACACATAgc                                                                                  <  1:357581/251‑1 (MQ=255)                                                                                                                                                                                     cccACAGTGCAGCGAAAAATTTCTGGCCCCCGTAAAGCACTCGCGGCGGCGGG‑TTCCTTTCCCC‑TTTTT‑CCCCCCCCGCCCCGCAACAAACCC‑ACACCGACA‑C‑TGGCACCATCTGCACACGCATTTTTGGTTTGGTACGGGGTAAATCAAGATATTAGTGTAACATCTTAATATACTCTTCTTTCGGCGTG‑TAATAGGTGATCCCACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCg                                                                 <  1:616771/250‑1 (MQ=255)                                                                                                                                                                                         caGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑GTTACTTTGCC‑TTTTT‑TACGCCCCCGGGCACGGCCAAACAACCCCCGCCA‑‑‑AAGCACCCACTCCCCATCCAGATTTTGGGTGGTACGGTGCAAATCGCGATAGTAGTGGGACAGCTTTATATACCCCTCTTTCGTCGGGGTAAAAGTGGAAAACAC‑CGGCCCCCCGCGAATGGACACATAACGCGGCAGCAGCGTGGCGtggg                                                             >  2:254559/1‑248 (MQ=255)                                                                                                                                                                                          aGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCCCCTCGCC‑TTTTT‑TACGCCCCCGCGCGCGGCCCAACA‑TCACCCACA‑CCTGGC‑CCATCTGCACATCCAGATTTTGGGTGCTACCGTGCAAAACGAGATAGTCGGGGGAAAGCTTAATAAACTCTTTTTTCCGCGGG‑TAATAGGGGAAAACACCCCGCCCCCCCTGGATGAACACAAAACGCGGCAGAAGCGGGcccgcggt                                                            >  1:176865/1‑245 (MQ=255)                                                                                                                                                                                            aaCAGCAAAAAATTTCTGTACGCCGTTAAGCAAGCGCAGCGGGGGT‑TTCCCTTCCCC‑TTTTC‑CCCCCCCCGGCGCGCGCCAAAACA‑CCCCCGACA‑C‑TGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGCCTCTTCTTTCGGCGTG‑TAAAAGGTGAACCCCCGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGAttt                                                          <  1:473177/250‑1 (MQ=255)                                                                                                                                                                                                              tACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑TACGCCCCCGCGCGCGACCAAACA‑ACCCCGACA‑C‑TCGCACCCCCTTCCCATCCCGATTTGGGGTTCTACGGGGCCAAACGAGATAGTCGGGGAACAACTTAATAAACTCTTCTTTCCGCGGG‑TTATAGGGGAACACACGCGGCCCCCCCGGTATGGACACAAAAAACGGCCGTAGCGTGCCGCCGATTTCCCGCgggtccaacaggg                                        >  2:384252/1‑240 (MQ=255)                                                                                                                                                                                                                      aaggCATCCGCAGCGGCGGT‑TTCCCTCCCCT‑TTTTT‑CCCCCCCCGGCGCGCGACAAAACA‑TCCCCGACA‑C‑GGCCACCATCGGCCCATCAAGTTTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTG‑TAATAGGTAACCCCACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAg                                <  1:500118/248‑1 (MQ=255)                                                                                                                                                                                                                                       ggg‑TTCACTTCGCC‑TTTTT‑CACGCCCCGGCGCGCGACCACACA‑TCACCGCCA‑C‑TGGCACCATCTGCACATCCAGCTTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTG‑TAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTg               >  1:278177/1‑251 (MQ=255)                                                                                                                                                                                                                                       ggg‑TTCACTTCGCC‑TTTTT‑CACGCCCCGGCGCGCGACCAAACA‑TCACCGACA‑C‑TGGCACCATCTGCACATCCAGATTTGGGCTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTG‑TAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCCCTAACGGGTCAGGCAAGTCATGTGCCAGACACTg               >  1:350255/1‑251 (MQ=255)                                                                                                                                                                                                                                                        tttt‑CACGCCCCGGCCCGCGACAAAACA‑TCACCGACA‑C‑TGGCCCCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTG‑TAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTc  >  2:95392/1‑249 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         TATTAAACAGATCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCC‑TTTTT‑CACGCCCCGGCGCGCGACCAAACA‑TCACCGACA‑C‑TGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTG‑TAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTC  >  NC_000913/4255812‑4256305

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