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

TGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAG  >  NC_000913/3875227‑3875821                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                             tggagtaCCAGCGTGATGACGTTTGCGTTTGCCGTGGGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAAAAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGtt                                                                                                                                                                                                                                                                                                        <  2:427069/298‑1 (MQ=255)     gtaCCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGATGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTcc                                                                                                                                                                                                                                                                                                     >  2:70771/1‑300 (MQ=255)      taCCAACGTGATGACGTTCGCGTTTGGCGTGCGTGTAATGTAGTTTAAAGTTGTCTTGTTGTACTACAATTTAGATCACAAAAAGAACAATTCATAAAAAAATGACAAGAGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGTTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATTAAGTTGAATTTTAAGGGATTTTTTAAGGATGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTcc                                                                                                                                                                                                                                                                                                     <  1:750365/299‑1 (MQ=255)              tgatgaCGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCaa                                                                                                                                                                                                                                                                                           >  2:364279/1‑301 (MQ=255)                   aCGTTCGCGTTTGCCGTGCGTGTAATGTAGTTCAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAAAACAATTCCTAAAAAAATAATATGCGTCGGGCAGAAATCTAAAAAGAGATATCAGGCGCTAAACAGGAGGGAAAGAAAAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAAGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGt                                                                                                                                                                                                                                                                                      <  1:364281/301‑1 (MQ=255)                                   tgcgtgTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAAGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAatat                                                                                                                                                                                                                                                                       >  1:202629/1‑300 (MQ=255)                                   tgcgtgTAATGTAGTACAAACTTATATGGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATc                                                                                                                                                                                                                                                                      <  2:189268/301‑1 (MQ=255)                                                gtaCAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCCCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGa                                                                                                                                                                                                                                                         >  2:663090/1‑301 (MQ=255)                                                    aaaCTTATATTGTTTTACTACAATTTAGCTCACAAAAAGAACAATGTATCAAAAAATGACATGCATCGGGAAGAAATCTGAAAAGAGATATCTGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAAGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGGTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACAcc                                                                                                                                                                                                                                                     <  1:70772/301‑1 (MQ=255)                                                            atTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAAGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAATATCAATCGCAGTCAGACACCATCACTgg                                                                                                                                                                                                                                             <  2:202628/301‑1 (MQ=255)                                                                 tgtACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCTGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCCAGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTAt                                                                                                                                                                                                                                        >  2:584996/1‑301 (MQ=255)                                                                       acAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAgg                                                                                                                                                                                                                                  >  1:499120/1‑301 (MQ=255)                                                                                  cacaAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAACAGATAGt                                                                                                                                                                                                                       >  2:175899/1‑301 (MQ=255)                                                                                   acaAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTa                                                                                                                                                                                                                      >  1:30182/1‑301 (MQ=255)                                                                                     aaaaaGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTaaa                                                                                                                                                                                                                    >  2:468580/1‑300 (MQ=255)                                                                                      acgctcttccgatctAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAAt                                                                                                                                                                                                                    <  2:460695/286‑1 (MQ=255)                                                                                       aaaGAACATAGCATAAAAAAATGACAGGAGTCGAGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACGGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATgg                                                                                                                                                                                                                  <  2:547166/301‑1 (MQ=255)                                                                                           aaCAATGCATAAAAAAATGACATGCGTCGGGCAGAAAACTAAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAATGTTCTTTCAACGGCTTAACTACACGTTTAAAGGATTAATCATGAAGATGAATGTTAAGAGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGaa                                                                                                                                                                                                                <  1:175900/299‑1 (MQ=255)                                                                                                    tAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTAt                                                                                                                                                                                                     >  1:122241/1‑301 (MQ=255)                                                                                                     aaaaaaaTGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATAGatcggaagagcaca                                                                                                                                                                                                    >  1:460697/1‑288 (MQ=255)                                                                                                      aaaaaaTGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGt                                                                                                                                                                                                   >  1:20318/1‑301 (MQ=255)                                                                                                                       gggCAGAAATCTGAAATGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGGTTTTGTAGGGTAATCg                                                                                                                                                                                   >  1:23686/1‑300 (MQ=255)                                                                                                                             aaaTCTGAAAAGGGATATCAGGCGATAAACAAGATGGAAAGAAGAGTATGCTTTCAACGGATGAGCTACTCGTTTAAAGGATTAATCATGAAGTTAAATTTTAAGGGATTTTTTAAGGCTGCGGGTTGATTCCCACTGACGATGATGCATTCAGGCTGTATCTCGTATGCACTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGa                                                                                                                                                                              <  1:663094/299‑1 (MQ=255)                                                                                                                                 cTGAAAAGGGATATCAGTCGCTAAACAGGAGTGAAAGAAGAGTAGTCGTTCAATGGCTTAGCTACTCGTTTAAAGGATTAATCATGATGTTGAATTTTAAGGGAATTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCt                                                                                                                                                                        <  1:713246/301‑1 (MQ=255)                                                                                                                                  tGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTGCCTGAACGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCtt                                                                                                                                                                       >  1:71941/1‑301 (MQ=255)                                                                                                                                  tGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAGAAAAAGATAGTCCAGGAAAGTATCAATTGATGCCAGATACCATCAATGGGCTATCGCCGTCAAAAGATATCAATGGAACAAAAGTCTATGATTTCGTACGGAATTCGCttgctcagt                                                                                                                                                                         >  1:76630/1‑291 (MQ=255)                                                                                                                                    aaatGGGATATCATGCGCTAAAAAGGAGGGAAAGAAGATTAAGCTTACAACGGCTTAGATACGCGTTTAAAGGATTAATCATGAAAGTGAATTTTGAGGGATTTTTTGAGGCTGCCGGTTTATTCCCACGGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGTAATCGCTGGATTACCTTa                                                                                                                                                                      <  2:23685/300‑1 (MQ=255)                                                                                                                                           atatCAGGCGATAAACAAGAGGAAAAGGAGAGTATGCTTTCAACGGCTTAGCTACTAGTTTAAAGGAGTAATCATGAAGTTGAATTTTAAGGTATTTATTAAGGATGACGGTTTATTCCAACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTAGCGCAGGCAAAAGATGGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGAt                                                                                                                                                              <  1:584999/301‑1 (MQ=255)                                                                                                                                             atCAGGCGCTAAAAAGGAGGTAAAGAAGCGTATGCGTTCAAAGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTGAGGGATTTTTTAAAGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTAACGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATg                                                                                                                                                             <  2:18251/300‑1 (MQ=255)                                                                                                                                               catgcgcTAAACAGGAGGGAAAGAAGGGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTAAAGGCTGACGGTTTATGCCCAATGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTATGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATgg                                                                                                                                                            <  1:460366/296‑1 (MQ=255)                                                                                                                                                    gcTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTTCCGAt                                                                                                                                                     >  2:520000/1‑301 (MQ=255)                                                                                                                                                     cTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGAt                                                                                                                                                     <  1:468582/300‑1 (MQ=255)                                                                                                                                                      tAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCATTGGGCTATCGCCGGCAAAAGTTAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCCCTGATGGGGCCGATGa                                                                                                                                                   >  2:711960/1‑301 (MQ=255)                                                                                                                                                       aaaCAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGATACTCGTTTAAAGGATTAATCATGAAGGTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGAc                                                                                                                                                  <  1:711965/301‑1 (MQ=255)                                                                                                                                                              aggGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTCTCCCATACCGCAAAGGGTAGTTCAGTAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGTGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATTGTTa                                                                                                                                           >  1:96762/1‑301 (MQ=255)                                                                                                                                                                                                 cTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTCGTTTGCCTGAAAACAAAAGCTAGTTCTGGAAAGTATCAATCGCAGTCAGTCCCCATCACTCGGCTATCGCAGTCAAAAGCTAGTAATGCAACAAAAGGCTATGCTTTTGTAGGGGAATCGCTGTATTAACTAATCACTGACGCTGCCGTTAACATCGTTAAGATGCTCAATCATTCATCACTTCACCGGCACaaa                                                                                                        >  1:199092/1‑300 (MQ=255)                                                                                                                                                                                                     tttAAAGGAGTAAGCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTTATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTGCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTc                                                                                                    <  2:122240/301‑1 (MQ=255)                                                                                                                                                                                                      ttAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGTAGTTCAGGAAAGTATCAATCGCAGTAAGACACTATCACTGGGCTATCGAAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGTTGACATTGTTAAGATGCTCAATGATCCTGCCCTTCACCTTCACAATATTCa                                                                                                   >  1:743856/1‑301 (MQ=255)                                                                                                                                                                                                         aaGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTCGTTTCCCAGACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGt                                                                                                >  1:274738/1‑301 (MQ=255)                                                                                                                                                                                                          aGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTTAAAACAAAAGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGCCATCGTTTAGATCCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGtt                                                                                               >  1:607087/1‑301 (MQ=255)                                                                                                                                                                                                              ttAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCCTTACAGAAAAAGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAACGGCTATGTTTTTGTAGGTGAATCGCTGGATTACCTTATCACTGCTGGTGCCGATTATATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGTACAATATTCAGGTTTCCg                                                                                           >  1:585088/1‑301 (MQ=255)                                                                                                                                                                                                                aaTCATGAAGTTGAATATTAAATGATTTTTTAAGGCTACCGGATTATTCCAACTTGCGCTGATGCATTCAGGCTATATCTCGTATGCTCTGGTTTCCCATACCAAAAAGGGTAGTTCAGGGAAGTATCAATCGCAGTCAGACACCATCACTGGACTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTGATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGACACAATATTCAGGTTGCCGAt                                                                                         <  2:96761/301‑1 (MQ=255)                                                                                                                                                                                                                     tGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTCTTCTGCCTTAAAACAAAACGCAGCTCAGGAAAGTACCAATCGCAGTCAGACACCATTCCTGGGTTATCGCAGACACCAGCTACTACTGGAACAAAAGGCTATGTTTCTGTAGTGGACTCCCTGTATTACCTTATCACTGACGGTGCCGATGACATCGTTAAGATGACCCATGATCCATCACTTAACCGGCTCACTATTAACGTCGCcgacgaag                                                                                     >  1:459558/1‑298 (MQ=255)                                                                                                                                                                                                                                 ttAAGGGATTTTTTAAGGCTGCCGGTGTATTCCCACTGGCGCTGATGCGTTCAGGCTGTATCTAGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGTGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGGTGCCGATGACGCAAGATTTGtttt                                                                        <  2:274736/301‑1 (MQ=255)                                                                                                                                                                                                                                                  gCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGTTTTCCTATACAACAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTTCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGa                                                       >  1:54482/1‑301 (MQ=255)                                                                                                                                                                                                                                                                cccACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTTTGCCTTACATCAAAGGGTAATTTAGAAAAGCATCAATCTCACTCCGACACCATCACTGAGCTATCGCAGGCAAAAGATAGCACTTGAACCAAAAGCACCGTTATTGTAGGGGAATAGCTGGATTCGCTCATCACTAATCGACCCTATCACCTCGTTAACTGTCTCAATTATCCATCTAATACCCGGCACCATATTCAGGTTGCCGATGATGTAAGATTCGTTCTAATTTCGGGGGTAAAGAAATTTACCGAcaca                                         >  1:496548/1‑301 (MQ=255)                                                                                                                                                                                                                                                                 ccACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATTGCAGTCAGACACCATCACTGGGCTACCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGGAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGACGACATCGTCAAGATGCTCAATGATCGAGCACTTAACCGGCACAATATTCAGGTGGACGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACAGGCACaa                                        >  2:527510/1‑301 (MQ=255)                                                                                                                                                                                                                                                                      ggCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATc                                    >  2:236700/1‑300 (MQ=255)                                                                                                                                                                                                                                                                        agttGATGCTGTCAGACAGTAGCTCGTATGTGATGGTTTACAAAACCGCAAAGGGTAGTTCAGGAAAGTAAAAATCGCAGTCAGACAACATCACTGGGGTAACGAAGGCAAAAGATAGGAATGGAACAAAAGGCTATGTTTGTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGAAATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGc                                  <  2:499117/297‑1 (MQ=255)                                                                                                                                                                                                                                                                              tGCTTTCAGGCTGTATCTCGTATGCCGTGTTTTTCCACACCGCAAAGGGTAGGTCAGGGAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAATTAGTCCGGGAACAAAAGGCTATGTTTTTGTAGGGGACACGCTGGATTACCTTCTTACTGATGGTGCCGACGACATCGTTAATCTGCCCAATGAACCCGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGCTTTGTTTCAAATTAGGGGAAAAAGAAATTTACCGGATCAATATTGCTTtact                            >  1:6471/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                ttttCAGGCTGTATCTCGTATGCTCTGGTTTGCAATAACGGAAAGGGTAGTTCAGGAAAGTATCAATCACACTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGTTCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGAAGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATTTCGCTTtactact                         <  1:697802/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                cGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACgaaga         >  2:745323/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                       tctGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGaaaaa  >  2:547082/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                       tctGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAg  <  1:547085/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                             TGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAG  >  NC_000913/3875227‑3875821

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