Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,875,521	0	T	C	22.2%	73.1 / 11.8	47	A28A (GCT→GCC)	yidX	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (20/15);  new base C (10/0);  total (30/17)
Fisher's exact test for biased strand distribution p-value = 1.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

AGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGAT‑TTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACG  >  NC_000913/3875236‑3875810                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                   aGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATAcc                                                                                                                                                                                                                                                                                     >  2:244621/1‑301 (MQ=255)           cGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTACAGGATTAATCATGACGTTGAATTTTAAGGGATTTTTTAAGGCTTCCGGTTTATTCCCACTGGCGTTGATGCTTTCAGGCGGTATCTCTGAGGCTCTGGTTTCCATTACCGCaacaggt                                                                                                                                                                                                                                                                            >  1:347304/1‑295 (MQ=255)            gTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGATGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTGCTTTATCGCAAAGAGtt                                                                                                                                                                                                                                                                           >  2:214166/1‑299 (MQ=255)              tCGCTTCTGCCGTGCGTGTAATTCAGGCAAAACTTATATTGATGTATTAAAAGTTAGACAAAAAAAAAAAAAATGAATAGAAAAATGAAATGCGACAGGCAGAAATCAGAAAAGGGAAATCAGGAACTAAAAAGGAGGGAAAGAAAAGAAAAATATAAAAAGATTAGAAAATCGTTTAAAAGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCAAAATGTCGCTGATGCTTTAAGGATGTATATAGTATGATATGGTTTCCCATACCGCAAAGGGTAGt                                                                                                                                                                                                                                                                         <  2:90713/300‑1 (MQ=255)                              tgtAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGTTTTCAACGGCTTAGCTACTCGTTTAAAGGATTTATCATGTAGTTGAATTTTAAGGAATTTTGTACGGCTGCCGCTTTATTCCCATTGGCGCGGATGCTTTCAGGTTGTATCTCGTATGCTCGGGTCGCCCTTTCCGCAACGGGTAGTTCAGGCAAGTAGCAAtt                                                                                                                                                                                                                                                        >  1:94964/1‑300 (MQ=255)                              ggtAATGAAGTACTAACGTTTATTGTCGTACTACAAATTATATCACAGTAAGATCAATTCAAAAAAAAAGGACATGCGTCGTGCAGATAAATGAAAAGTGAAATCAAGAGCTAAACAGGAGATAAAGCAGAGTATGCGTTCAACTGCGTAGTTCCTCGTTTAAATGATTAATCATGAAGTTGAAATTTAAGGCATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCATGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAAGAAAGTATCAATc                                                                                                                                                                                                                                                        <  2:418929/300‑1 (MQ=255)                                            aacagatatTGTTGTACTACAAGTTAGAGCACAAAAATAACAATGCATAAAAAAATGACATGGGACGGGCAGAAATATGAAAAGGGATATAAGTCGCTAAACAGAAGGAAAAAAAGAGAAAGATTTGAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCa                                                                                                                                                                                                                                          <  2:86971/296‑1 (MQ=255)                                                    ttgttgTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTc                                                                                                                                                                                                                                                                                            >  1:276926/1‑243 (MQ=255)                                                    ttgttgTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAAGGACATGAGTCGTGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCGTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTc                                                                                                                                                                                                                                                                                            <  2:277008/243‑1 (MQ=255)                                                         gtACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCCCTGATGCTTTCAGGCTGTACCTCGTATGCTCTGGTTTCCCATACCGTAAAGGGTAGTTCAGGAAAGTACCAATCGCAGCCATACACCATCACTGGGCTATc                                                                                                                                                                                                                             >  1:50511/1‑301 (MQ=255)                                                         aaaacaaaaattataacacaAAAAAAAAAAAGCAAAAAAAAATGACATACGTCCGGCAGAAAGCTGAAAAGGGAAATCAGGCGTGAAACAGGAGGGAAAGAAAAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATc                                                                                                                                                                                                                             <  1:244548/285‑1 (MQ=255)                                                                          acaAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTcc                                                                                                                                                                                                                                                                                           >  1:45442/1‑222 (MQ=255)                                                                          acaAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAAGCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTcc                                                                                                                                                                                                                                                                                           <  2:45454/222‑1 (MQ=255)                                                                            aaaaaGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACAACATCACTGGGCTATc                                                                                                                                                                                                                             >  1:113246/1‑282 (MQ=255)                                                                            aaaaaGAACAATGCATAAAAAAATGACATGCGTCAGGCAGAAAACTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAATATGCTTTCAACGGATTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATc                                                                                                                                                                                                                             <  2:113274/282‑1 (MQ=255)                                                                              aaaGAACAATGCATAAAAAAATGACAGGATTCGGACAGAAATCTGTAAAGGGATATCAGGCGCTAAACAGAATGGAATAAGTAATATGCATTAAACGGCTAAGCTAATAGTATAAAAGATAAATCATACAGTTGAAGTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATgg                                                                                                                                                                                                        <  2:94988/301‑1 (MQ=255)                                                                                   aCAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGATAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCCGGCAAAAGATAGTAATGGaatc                                                                                                                                                                                                    >  2:82161/1‑298 (MQ=255)                                                                                    cAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGTTTTGCTATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACTCCTTTACTGGGCTATCGCAGGCCAAAGATAGTAATGGAACaa                                                                                                                                                                                                   >  2:172991/1‑300 (MQ=255)                                                                                        gCAGAAAAAAAGGACAGGGGTAGGGGAGAAAGCGGAAAAGGGATAACAGGGGATAAACGGGAGGGAAAGAGGACTATCCTTCCAACGGATTAGCTACACGTTAAAAGGATTAATCATTAAGTTGAATTTTAAGGGATTTTTGAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGc                                                                                                                                                                                              <  2:50527/301‑1 (MQ=255)                                                                                                                          gAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCttt                                                                                                                                                            >  2:255514/1‑300 (MQ=255)                                                                                                                          gAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTATAAGGGATTTTTTAAGGCTGCCGGTTTATTCACACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTATGGTTTCCCTTACCGCAAAGGGTAGTTCAGGAAAGTATCCATCGTAGACAGATACCCTCACTGGGCTATCGCTGGTTAATTATAGTAATGGCTAAATAGGCTATGtctgttatggtgcattgacggattactttt                                                                                                                                                            >  1:158164/1‑272 (MQ=255)                                                                                                                                       aGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTTGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTTACCGAAAAAAGTAGTTTAGGAAAGTATCCATCGCAATCAGACACCAACACTGCGCTATCGCAGGCCAGAGTTAAAAATGGCACAAAAGGGTATATTTTAGTAGGGGTAAAGCTGAGGTACTTACTCACTGATgggac                                                                                                                                               >  2:447009/1‑298 (MQ=255)                                                                                                                                        ggCGCGAAACAGGAGGGAAAGGAGAGAGTGCTTTCAACGGTGTATCTACTCGTTGAAAGAGTTTATCATGAAGTTTAATTTTAATGGCTTTGTTAAGGCGGCCGGTGTATTCCCACTGGCGCTAATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGGAAAGGGTAGTTCAGGAAAGTATCAAGCGCAGTCAGACACCATCACTCTGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGcc                                                                                                                                              <  1:172941/301‑1 (MQ=255)                                                                                                                                          cgcTAAAGCGGAGTGAAGGAAGAGTTAGCTTTAAACGGTCTAACTTCTCGTTTAAAGGATTAATCATGAAGTTGAATGTTAAGGGATTTTTTACGGCTTCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGa                                                                                                                                            <  1:255438/301‑1 (MQ=255)                                                                                                                                               atcaggactgaccgaagaGTGAGATTTTACCGGATTAACTGATCGTGTAATTGGTTAATCATGTAGTGAAATTTTAAGGCGTTGTAGTAGGCTTCAGGATTATCCCCCCTGACGCCATTTCTTTCCGACTTTATCTAGTAAGCTTGGGTTTCCCATACCGCAAAGGGTAGGTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGTGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGTATTACCTTATCACTGATGGTGCCGATGACa                                                                                                                                       <  2:122465/288‑1 (MQ=255)                                                                                                                                                aCAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGCATCGCTGGATTAACTTATCAATGAGGGTGCCGATGACAt                                                                                                                                      >  2:283727/1‑301 (MQ=255)                                                                                                                                                   ggaggGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAGCAAAGGGTAGTTCAGGAAAGTTTCAATCGCAGTCAGACACCAGCAATGGGCTATCGCAGGCAAAAGCTAATAATGGACCAAAAGGCTATGTTTTTGTATGTGAATTGCTGGATTACATTATCACTGATGGTGACGATGACATCGt                                                                                                                                   >  2:231515/1‑301 (MQ=255)                                                                                                                                                      gggAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGAAGGAAAAAGATAGTAATGGAACAAAATGATATGTTATTGTAAGGGAAGCACTGGGTTACCTTATATATGATGGAGccagtgaagtcgtca                                                                                                                                 >  2:284264/1‑287 (MQ=255)                                                                                                                                                        gAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCAGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCTTACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCATTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGAAATCGCTGGATTAACTTATCAATGATGGTGCCGATGACATCGGTacg                                                                                                                               >  2:316471/1‑298 (MQ=255)                                                                                                                                                             aagatgtagCTTGCAAGGGGTGAGCTACTCGTTTAAATGATTTATTATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTAATCCCACTGGAGATGATGCTTACAGGCTGTATCTAGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAAACGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTc                                                                                                                         <  2:318978/293‑1 (MQ=255)                                                                                                                                                                     gCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTTTCCTATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACATTATCATTGATGGTGCTGAGGCCATCGTTAACATGCTCCATAATcc                                                                                                                 >  2:278743/1‑301 (MQ=255)                                                                                                                                                                          acacttgtgatgttctagTTTCAAAGATTAATCATGGATTTGTATTTTAATGGATTTTTTAGGTCTGCCGGTTTATTCCCAATGGCGCTGATTCTTTCAGTCTGTATTTCTTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGTGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCAc                                                                                                            <  1:283644/284‑1 (MQ=255)                                                                                                                                                                                        cTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCCTGACCAAAAAGGGTAATTCAGGAAAGTAACAATAGAAGTCAGACACCATCACTGGGCTATCGCCGGCAAAAGATATTAATGGAACGAAAGGCAATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGACACaa                                                                                               >  2:167865/1‑300 (MQ=255)                                                                                                                                                                                          cGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTCTTCTGCCAGACCACAAAGCGTAGTTCAGGGAAGTATCAATCTCAGTAAGTCACCATTAATGGGCTATCGCTGGTAACAGATAGTAATGGAAGAAAAGGATATGTTTCTGTAGGGGAATGGCTGTATTACCTTAGCAGTAGTGGTGCCGTTGACATCGATACGATGCTCACAGACCCAGCAccaacgcggacaagcat                                                                                            >  2:323692/1‑285 (MQ=255)                                                                                                                                                                                           gTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGACATGTTTTTGTAGGGGAATTGCTGGATTAACTTAGCCATGGTTGTGCCGATGAACTCGTTAAGTCGGTCAAAGATCTAGCACTTAACCGGTCCAattat                                                                                           >  2:73823/1‑298 (MQ=255)                                                                                                                                                                                                            tGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTTTCCCATACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCCGCACTTAACCGGCACAATATTCAGGTTGCCTATGACg                                                                           >  2:256594/1‑300 (MQ=255)                                                                                                                                                                                                                          taGGGATTTTTGTTGGCTCCCTGTTTATTCCCACTTGGTCTTTTGCTTTCCTCCTTTATCTCGTATGTGTTGTTTTCCCATACCGCAAATGATAGTTCAGAAAAGAATCAATTGCATTCAGACACCATCACTGGGCTATCGCAGGCCAAAGATAGTAATGGAACAAAATGCTATGTTTTTGTAGGGAAATCGCTGGATTACCTTATCACTTATGGTGCAGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGAT‑TTGTTTTa                                                            <  1:316382/299‑1 (MQ=255)                                                                                                                                                                                                                            gggATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGCTGGTGCCGATGCCATAGTTAAGATTCTCAATGAACACGCAGTTAACCGGCCAAATATTCGTGTTGCCTACGACGCAAAGT‑TTGTTTTAAAt                                                         >  1:230094/1‑301 (MQ=255)                                                                                                                                                                                                                            gggATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTTTCCCTGACAGAAAAAGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATAACTGGGCTATCGCAGGCAAAAGATAGAAATGGAACAAAAGGCCATGATTTTGTCGGGGAAACGCTGGATTACCTTCACACTGATGGTGCCGATGACATCGTTCAGATGCTCAATGATCAAGGACTTAACGCGCCCAATCCTATGGTTGACGATGACGCAACATGttgttgaaa                                                           >  2:386526/1‑296 (MQ=255)                                                                                                                                                                                                                             ggatatgttactgtttccGATTTATTCACACTGGCACTGATGCTGTCATGCGTTAGCTCGTCTGTGCTTGTTTCCCCTACCGTAAAGGGTAGTTCAGTAAAGAATCATTCGCAGTCAAACGACCTCCCTGGGCTATCCCAGGAACAACATATCAATGTAACAACATGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCAAGCACTTAACCGGCACAATATTCAGCTTGCCGATGACGCAAGAT‑TTGTTTTAAATg                                                        <  1:231445/285‑1 (MQ=255)                                                                                                                                                                                                                               aTTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTCTGCCATACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCCCTTAACCGGCACAATATTCAGGTTGACGATGACGCAAAAT‑TTTGTTTAAATgtg                                                      >  2:267706/1‑299 (MQ=255)                                                                                                                                                                                                                                        ggcgggcGGGTTATTCCAACTGGCGCTGATGCTTTCAGGGTGGATTTTGTATGCTCTGGATTCCCATACAGCAAAGGGTAGTTCAAGAAAGTATCAAACGCAGTCAGACACCATCACTGGGATATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGAT‑TTGTTTTAAATGCGGGGAAAAAg                                             <  2:347399/295‑1 (MQ=255)                                                                                                                                                                                                                                                    aTTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTGAAGATGCTCAATGATCCAGCACTTAACCGGCACAGTAATGTGTTTGCCGATGCCGCAAGAT‑TTGTTTTAAATGCGGGGAAAAAGAAATTTACAGGc                                 >  1:270659/1‑301 (MQ=255)                                                                                                                                                                                                                                                       cccACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGACGATGACATCGTTAAGATGATCAATTATCCAGCACTTAACCGGCAAAATATTCAGGTTGCCGATGACGCAAGAT‑TTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGcac                               >  1:247141/1‑300 (MQ=255)                                                                                                                                                                                                                                                                      gCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCTCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAATGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGTTGCCGATTACATCGTTAAGATGCTTAATGTTCCAGCACTTAAACGGGACAATACTAAGGTAGACGATTACGAAAGAT‑TTGTTTGTATTGCAGGTAAAAAGACATGTAcagtcaaccattcgttttccttt               >  2:369429/1‑279 (MQ=255)                                                                                                                                                                                                                                                                         ttCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGAAACCATCACTGGGCTATCGCAGGCAAAAGATAGGAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTGCCTTATCACCGCTGGTGGCGATGACCTCGTTAAGATGCTCAATGATACAGCACTTAACCGGCACAATCATAAGGTTGCCGAATAGCCAAGAT‑ATGTTTTAAATGCGGGGAAAAAGAAATTTAAAGGCTCAAGATCGCTTtattattac            >  1:350245/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                    tctcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGAT‑TTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACg  >  2:340445/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                   AGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGAT‑TTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACG  >  NC_000913/3875236‑3875810

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