Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,286,870	0	C	A	27.6%	97.5 / 29.1	58	S54S (TCC→TCA)	ychS	putative uncharacterized protein YchS
Reads supporting (aligned to +/- strand):  ref base C (20/22);  new base A (10/6);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCT  >  NC_000913/1286580‑1287159                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                  gTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCtttt                                                                                                                                                                                                                                                                                          >  2:660327/1‑300 (MQ=255)      gggTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGACTTTTGaaa                                                                                                                                                                                                                                                                                      >  1:550097/1‑299 (MQ=255)      gggTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCAGGGACTTTTGAAAg                                                                                                                                                                                                                                                                                     <  1:336678/300‑1 (MQ=255)      gggTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGATTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGaaa                                                                                                                                                                                                                                                                                      >  1:276093/1‑299 (MQ=255)      gggTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATT‑AAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTACGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCCCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCCTCTCATCCTTCCCCGCTTGGGCCGAATGTTTAATTGCGAACTCGTTTGGGAA‑TTCAGGGGCTTTCGAAAGt                                                                                                                                                                                                                                                                                    >  2:48300/1‑300 (MQ=25)                          ttcgcttATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGTAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                                                                                                                                                                                                 <  2:684710/300‑1 (MQ=255)                                tATCCGACCTACAGTGTAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                                                                                                                                                                                                                           <  1:272881/300‑1 (MQ=255)                                tATCCGACCTACAGGGGAGGATATTGTATGCCCGGTAAACGTGATGCCATCGGGTAGAATTAGCTCGGAGTAACAGGTTTTGATTATTGAAGTAAGCGATGTTGGTGGGGGAAGGATTACTCATCGCTGCGCGCTTCGCCCTTCTGTTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAATTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                                                                                                                                                                                                                           <  2:199217/300‑1 (MQ=21)                                tATCCGACCGACAGGGGAGGATATTGTATGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGATTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                                                                                                                                                                                                                           <  2:460793/300‑1 (MQ=255)                                     gACCTACAGGGGAGGATATTGTAGGCCCGTTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                                                                                      <  1:532334/300‑1 (MQ=255)                                     gACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCTCAgcgc                                                                                                                                                                                                                                                      >  1:277870/1‑300 (MQ=255)                                     gACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                                                                                      <  1:532691/300‑1 (MQ=255)                                     gACCTACAGGGGAGGATATGGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCGCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAG‑TTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                                                                                      <  2:655170/300‑1 (MQ=17)                                      aCCTACAGGGGAGGATATTGTAGGCCCTGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTGGATTATTTAAGGAAGCGATGTTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                                                                                      <  2:22245/299‑1 (MQ=255)                                      aCCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCANAATATCTGATTGCGGATTCGTTTGAGAA‑TTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCt                                                                                                                                                                                                                                                     >  1:531604/1‑300 (MQ=255)                                          aCAGGGGAGGATATTTTAGGCCCGGTAAACGTGTTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTT‑GGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg                                                                                                                                                                                                                                                                  <  2:348443/283‑1 (MQ=14)                                          aCAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTT‑GGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg                                                                                                                                                                                                                                                                  >  1:348443/1‑283 (MQ=17)                                          aCAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgc                                                                                                                                                                                                                                                   >  1:348277/1‑298 (MQ=255)                                                  ggaTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTGATTGCGGATTCGTTTGGGAA‑TTcc                                                                                                                                                                                                                                                                                                   >  2:136567/1‑242 (MQ=255)                                                  ggaTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAG‑TTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                                                                                      <  2:425928/287‑1 (MQ=14)                                                  ggaTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAG‑TTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                                                                                      >  1:425928/1‑287 (MQ=14)                                                       ttGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAA‑TTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                                                                                                                                                                                                                     >  1:440878/1‑299 (MQ=255)                                                       ttGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTT‑GGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcaccccctctacgcccctcacccccaa                                                                                                                                                                                                                                    >  1:716701/1‑273 (MQ=14)                                                       ttGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTAATCCGTGCTGCGCGTTTCGCACTTCGGGTCACTTCCTGCGGCAACGCCCTCTCGCTGGTGCTCACGTCGAAACTTGGTCGAAGCTTACCCTCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGGGGA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGACCtt                                                                                                                                                                                                                                      >  1:167999/1‑298 (MQ=21)                                                       ttGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATAGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCAGCTTGGGCAGAATATTTAATTGCGGATTCGCTCGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                                                                                                                                                                                                                     >  1:441155/1‑299 (MQ=255)                                                            ggCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCg                                                                                                                                                                                                                               <  1:280552/300‑1 (MQ=255)                                                                cGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGc                                                                                                                                                                                                                           >  1:119261/1‑300 (MQ=255)                                                                 ggTAAACTTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAG‑TTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGACTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGcc                                                                                                                                                                                                                          <  1:239238/300‑1 (MQ=255)                                                                 ggTAAACGTGGTGCCAGCGTGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGTTGGTGTGGGAAGGACTACTCAGCCCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAG‑TTCAGGGACTTTTGAAAGTGATGGTGGTGGTGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGcc                                                                                                                                                                                                                          <  2:380497/300‑1 (MQ=11)                                                                             gCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGc                                                                                                                                                                                                               <  1:648469/299‑1 (MQ=255)                                                                                                   gAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcaccccttcga                                                                                                                                                                                                                                                    >  2:603907/1‑229 (MQ=255)                                                                                                   gAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTG‑GAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcacaccttcgacgt                                                                                                                                                                                                                                                 >  1:709962/1‑229 (MQ=14)                                                                                                                       ttAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCCTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAATCCgcgc                                                                                                                                                                                                                                                      >  2:261298/1‑218 (MQ=255)                                                                                                                       ttAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTAATTGCGGATTCGTTGGAAAG‑TTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGAttctcccattc                                                                                                                                                                                                                                                      >  2:257354/1‑209 (MQ=11)                                                                                                                                 gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTcccc                                                                                                                                                           <  1:348085/299‑1 (MQ=255)                                                                                                                                 gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCGCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTcccc                                                                                                                                                           <  1:347744/299‑1 (MQ=255)                                                                                                                                  aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCACCGCTCTCTCGCTGACGCTCGAGc                                                                                                                                                                                           >  1:130912/1‑266 (MQ=255)                                                                                                                                             ggAATGATTACTCAGCGCCGCGCACTTCGCCCTTTGTTTCGTTGCCTGCGTCAACGCTCTCTCGCTGGCGCTCAAGTCAAACCTTGGTCGAAGCTTCTCACCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGACGATGGTGGGGGAAGGATTAccc                                                                                                                                                                                                                                                           <  2:199341/191‑3 (MQ=255)                                                                                                                                                 ggATTACTCAGCGCTTCGTGCTT‑GCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCNTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTg                                                                                                                                         <  2:220686/300‑1 (MQ=255)                                                                                                                                                    ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAt                                                                                                                                       <  1:13302/300‑1 (MQ=255)                                                                                                                                                      aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGCTGCCTGCGTCAACGCTCTCTCGCTGACGCTCAAGCCGAACATTAGTCGAAGCTTCTAATCCTTc                                                                                                                                                              >  2:448036/1‑275 (MQ=255)                                                                                                                                                      aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTAATCCTTCCCCGCATGGGCAGAATATTTGAtt                                                                                                                                      >  1:564707/1‑299 (MQ=255)                                                                                                                                                      aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGTGCTGCGCGCTTAGccc                                                                                                                                                                                                                                        >  2:103820/1‑201 (MQ=255)                                                                                                                                                                   cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCTCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                                                                                      <  2:713433/174‑1 (MQ=255)                                                                                                                                                                   cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCTCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                                                                                      >  1:713433/1‑174 (MQ=255)                                                                                                                                                                                          cTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCacc                                                                                                                                                                  >  1:588985/1‑233 (MQ=14)                                                                                                                                                                                          cTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCacc                                                                                                                                                                  <  2:588985/235‑3 (MQ=14)                                                                                                                                                                                                     ctctctCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggt                                                                                         >  1:403047/1‑297 (MQ=255)                                                                                                                                                                                                             cTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAACGCTACGCGCCTAGCCCTTCAGGTCACTGCCTACGGCAACGCTCTCTCGCTGACGCTCGAGCAGGACCTTAGTCGAAGCCTCTCATCCTTCCCCGCATGGGCAGAATACTCGATCACAGATTCGTTTGAGAGTTCAGGGATTTTTGAAAGTGAtggtggt                                                                                         >  1:342195/1‑289 (MQ=17)                                                                                                                                                                                                              tGGCGCTCAAGTCAAACCTTGTGCGAAGCTCCTCATCCTTCCCCGCTTGGTCAGAATATTTGATTGCGGATTCGTGTGAGAA‑TTCCGGGGCTTGTGAAAGTGATGGTGGTGGGGTAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTTCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGTGGAAGGAt                                                                              <  2:161276/299‑1 (MQ=17)                                                                                                                                                                                                                                        aGCTTCTCATCCTTACCCGCGTGGGCAGAATATGGGATTGCGGAGTCGTTGGAGAA‑TTCCGGGGCTTTGGAAATTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAATTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGccc                                                      <  2:456586/297‑1 (MQ=18)                                                                                                                                                                                                                                           ttCTCATCCTTCCCCGCTTGGGCAGAATATTTAATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGc                                                                                                                                                                                                                       >  2:685702/1‑133 (MQ=255)                                                                                                                                                                                                                                                                gCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGTGCTTTTGAAAGTGATGGTGGTGGGGGAAGGACTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTTATGGTGGTGGGGGAAGGATTACTc                                                                         <  2:593283/254‑1 (MQ=255)                                                                                                                                                                                                                                                                     atatTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAAcc                                                                                                                                                                                     <  2:421040/141‑1 (MQ=14)                                                                                                                                                                                                                                                                     atatTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAAcc                                                                                                                                                                                     >  1:421040/1‑141 (MQ=14)                                                                                                                                                                                                                                                                     atatTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTAATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGGGAGTTCAGGGACTTTTGAGAGTGACGGTGGTGGGGGAAGGAtt                                                                             >  2:356774/1‑245 (MQ=255)                                                                                                                                                                                                                                                                     atatTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCAAAGCTTCTCATCCTTCCCCGCATAGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGTGCGCTTAGCACTTCGGGTCGTTGCCTGCGGCAACGCtctctc                       >  1:22216/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                 ggATTCGTGTTAGAA‑GTCCGTGACTTTCGAAAGGGATGTTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGACGCTCGAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgc                                                               <  2:124505/247‑1 (MQ=14)                                                                                                                                                                                                                                                                                         ttGAGAG‑TTCCGGTGCTTTTTAAATTGATGTTGGTGGGGGAAGGACTACTCATCGCTGCGCGCTTCGCCCTTCGTGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAATTGATGGTGGTGGGGTAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGTTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGATTCGAACCt  <  1:567184/300‑1 (MQ=17)                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                  GTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAA‑TTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCT  >  NC_000913/1286580‑1287159

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