Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,875,521	0	T	C	23.8%	113.3 / 15.1	65	A28A (GCT→GCC)	yidX	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (24/24);  new base C (15/0);  total (39/26)
Fisher's exact test for biased strand distribution p-value = 4.24e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAA  >  NC_000913/3875274‑3875808                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                gtaCAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGa                                                                                                                                                                                                                                            >  1:60306/1‑301 (MQ=255)        cTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGTTTTCCCATACCGCAAAGGTTAGTTCAGGACAGCATCAATCGCAGTCCGACACCAt                                                                                                                                                                                                                                      >  1:271328/1‑300 (MQ=255)        cTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTTTCCTTTACAGAAAAAGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCAt                                                                                                                                                                                                                                      >  1:976609/1‑300 (MQ=255)              ttgttgTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCGTCTGCTTGACAGAAAAAAGTAGGTCAGGAAAGTATCAATCGCAGTCAGACACGATAACgggg                                                                                                                                                                                                                               >  1:523076/1‑301 (MQ=255)                gttgtACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTGTTCTGCTATACAGCAAAGGGTAGTTCAGGAAAGTATCAATAGCAGTCAGACACATTCACTggg                                                                                                                                                                                                                               >  1:388590/1‑299 (MQ=255)                   gtACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCCCTGGGCTAt                                                                                                                                                                                                                           >  1:194128/1‑300 (MQ=255)                      ctacAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTGCCATACCGAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGGCAGACACACTCACTGGGCTATGGCa                                                                                                                                                                                                                       >  1:189752/1‑301 (MQ=255)                                     aaaaaagaaaaagGCAGAAAAAAATTACATGCGACGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAAAAGTAGGGAAAGAGGAATATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTaa                                                                                                                                                                                                        <  1:859224/300‑1 (MQ=255)                                             aCAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGTTTTGCCATACCGCAAAAGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAATATAGTACTGGAACaa                                                                                                                                                                                                >  1:971920/1‑301 (MQ=255)                                                aTGCATAAAAAAATGACATGCGTTAGGCAGAAGTCTGAAAAGGGATATGAGGCGCTAAAAATGAAGGAAAAAATTATATGCTTTTAACGGAATAGATACTCGATTAAAGAATTAATAATAAAGTTTAATTTTAAGGGATTTTTAACGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATAGAACAAAAg                                                                                                                                                                                             <  1:690854/301‑1 (MQ=255)                                                  gCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAANGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTATAACACAAAGAGTAGTTCAGGAAAGTATCAATCGCATTCAGACACCATCACTGGGCTATCGCAGTCAAAAGATGGTAATTGAACAAACgg                                                                                                                                                                                            >  1:1016735/1‑300 (MQ=255)                                                           aaTGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGtn                                                                                                                                                                                                                                                                     >  2:874473/1‑217 (MQ=255)                                                            ntGACAGGCATCGGGCAGAAATCTGAAGAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCTCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTc                                                                                                                                                                                                                                                                    <  2:689714/217‑1 (MQ=255)                                                              gACGTGCGGCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAAGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGt                                                                                                                                                                                <  1:827349/300‑1 (MQ=255)                                                                   gCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTgg                                                                                                                                                                                                                                                                                             >  2:645753/1‑186 (MQ=255)                                                                   gCGTCAGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGAAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTgg                                                                                                                                                                                                                                                                                             <  1:645776/186‑1 (MQ=255)                                                                                   tggaataggatatcactttgtaaaaccgaggGAATGAAGTTCATGTTTTCAATTGCTTAACTATACTTTTAAAGTATTATTTATGAATTTAATTTTGAAAGCATATTTTAAGGCTGCCGATTTATTCCCACTGGCGCTGCTGCTTTCAGGCTGTATCTCGTATGCTCTAGTTTCCCGTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAAACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCtt                                                                                                                                                          <  1:101566/274‑1 (MQ=255)                                                                                     aaaaGGGATATCAGGGGCGAAACAGGAGGGAAAGAAGAGTATGCGTTCAACGGCTTAGCTACTAGTTTAAAGGATTAATCATGAAGTTGAATTATAAGGGATATTTTAAGGCTGAAGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTAcc                                                                                                                                                            <  1:600774/297‑1 (MQ=255)                                                                                         gagtagacAGGCGCTGAACGGGAGGGAACGGAGAGTATGCTTTCAACGGCTTAACTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCGGCCGGTTTATTCCCCCTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCACATACCGAAAAGGGTAGTTAAGGAAAGTATCAATCGCAGTCAGACACAATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCAc                                                                                                                                                     <  1:450111/293‑1 (MQ=255)                                                                                              atCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTGCCATACCGCAAAGCGTAGTTCAGGAAAGTATCAATCGCAGTCCGACACCATCCATGGCCTATAGCAGGAAACAGATAGTGATGGAACAAAAGGCTATGTTTTTGTCGGGGAGTCGCTTGAATCCCATATCACTGATgg                                                                                                                                               >  1:879692/1‑301 (MQ=255)                                                                                                   gcgcTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTGCTATACCGCAAAGGGTAGTTCAGGAAAGTATCAATc                                                                                                                                                                                                                                                     >  2:304956/1‑194 (MQ=255)                                                                                                   gcgcTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTggg                                                                                                                                                                                                                               >  2:450095/1‑216 (MQ=255)                                                                                                   gcgcTAAACAAGAGGGAAATAAGAATATGCTGTCAACGGATTAGCTACGCGTTTAGAAGATTAATTATGGAGATGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGATTTCAGTCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATc                                                                                                                                                                                                                                                     <  1:304964/194‑1 (MQ=255)                                                                                                               aggGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTCGTCTCCTATAACGCAAAGGATAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACACACGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGttt                                                                                                                              >  1:44845/1‑300 (MQ=255)                                                                                                                 ggAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTCCCAGAAAGCAAAAGGTAGTTAAGGAAAGTATCAATCGAAGTTACACACCCTCACTGTGCTAACAGATGAAAAAGTTAGCAATGGCACAACAGGCCAAGTTTTTGCATGGGAATAACTGTCTTACCTTATTACTTATggatcagaacgcatagtttaa                                                                                                                            >  1:755078/1‑282 (MQ=255)                                                                                                                  gAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGGAGTTCAGGAAAGTATCAATCGCAGTCAGACACAATAACTGGGCTATCGCAGGCAAAAgn                                                                                                                                                                                                              >  2:767969/1‑217 (MQ=255)                                                                                                                  gAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACAGCAAAGGGTAGTTCAGGAAAGTATCAATCTCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATTGCTGGATTACCTTATCACTGATGGTGCCGATGAATTCGTTAAGa                                                                                                                           >  1:355613/1‑301 (MQ=255)                                                                                                                  gAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTCTCCCATACCACAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAACGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCAATGATGGTGCCGATGACATCGTTAAGa                                                                                                                           >  1:431065/1‑301 (MQ=255)                                                                                                                  gAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCATATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAAGCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAgn                                                                                                                                                                                                              >  2:731884/1‑217 (MQ=255)                                                                                                                   aaaGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTCTGCCATACCGCAAAGGGTAGTTAAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGAAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTCGGGGAATCGCCGGATTACCTTATCACTGATGGTGGCGATGACATCGTTAAGAt                                                                                                                          >  1:725444/1‑301 (MQ=255)                                                                                                                   aaaGAAGAGTATGCTTCCAACGGCGTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGAt                                                                                                                                                                                                             <  2:422975/218‑1 (MQ=255)                                                                                                                          agTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGTTTTCCCATACAAANAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGCACAAAAGGCTATGTTTTTTTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGAACTGGTTAAGATGCTCa                                                                                                                     >  1:1011566/1‑299 (MQ=255)                                                                                                                            tatacgttCAACTGTTGATCTACTCGAATAAAGGATTAATCAGGTAGTTGAATTTTAAGGGATTTTGTAAGGCTTCCAATATATTCCCACTGGCGCTAATACTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCAGACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGa                                                                                                                 <  1:400890/295‑1 (MQ=255)                                                                                                                               gCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGaa                                                                                                                                                                                                   >  2:794191/1‑216 (MQ=255)                                                                                                                                     naatgCTTAGATGTTCGATTAAAGGATTAATCGTGAAATTGAATTTGAAGGTATTTTTTAAGTCTGCCGGTTTATGCCCACTGGCGCTGATGCGTTCAGGCTGTAGCTCGTATGCTCTGGTTTCCCTTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAATGATAGTGATTGAACAAAAggg                                                                                                                                                                                           <  2:1030907/214‑2 (MQ=255)                                                                                                                                        ggCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCACAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTAt                                                                                                                                                                                        >  2:246974/1‑218 (MQ=255)                                                                                                                                        ggCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTCTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTaa                                                                                                     >  1:559924/1‑301 (MQ=255)                                                                                                                                          cTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTCGTCTGCTTTACCGAAAAGGGTAGTTCACGAAAGTATTAATCGTAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATCGAACAAAAGGCCATTTTTTTGTAGGGGAATCGCTGGATTACCTTATTACTGATGCTGACGATTATATCGTTAATATGCTCAATGTTCCAGctacttcc                                                                                                    >  1:64674/1‑293 (MQ=255)                                                                                                                                           ttAGCTAGTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGtt                                                                                                                                                                                     <  2:125126/218‑1 (MQ=255)                                                                                                                                               atACTCGTTTAAAGGATTAATCAAGAAGTTGAATTTGAAGGGATTTTTTAAGGCTGCCGGTTTATTACCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTAAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTg                                                                                                                                                                                 <  2:194124/217‑1 (MQ=255)                                                                                                                                                    cGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAatat                                                                                         >  1:758252/1‑301 (MQ=255)                                                                                                                                                    cGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAg                                                                                                                                                                                                           >  2:627901/1‑187 (MQ=255)                                                                                                                                                    cGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAg                                                                                                                                                                                                           <  1:627924/187‑1 (MQ=255)                                                                                                                                                    cGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTCGTCTGCCATACCACTAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATAGTTAAGATGCTCAATGAACCAGCACTTAACAGGCACAatat                                                                                         >  1:768991/1‑301 (MQ=255)                                                                                                                                                         aaaaGATGAATCATGAAGTTGAATTTTAAGGGATTGTTTAAGGCTGTCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGAAAAGGGTAGTTCAGAAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATTACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAgg                                                                                    <  1:146013/301‑1 (MQ=255)                                                                                                                                                            ggATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGttt                                                                                 >  1:191150/1‑300 (MQ=255)                                                                                                                                                             taTTAATCGAGAAATTGAAGTTTAAGGTATGTTATAAGGCGGCCGTTTTATTCCCACGAGCACTGATGCTTTCAGGCTATATCTAGTAAGCTCTGGTTTCCAATACCGCAAAGGGTAGTACAGAAAAGTATCAAACACAGACAGACACCATCAATGGGCTATCGCAGGAAAAAGATAGTAATGAAACAAAAGGCTATGTAGTTGTAGGGGAATCACTGGATTAACTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAgg                                                                                    <  1:972116/296‑1 (MQ=255)                                                                                                                                                             gATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCACAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCgatn                                                                                                                                                                   >  2:972068/1‑215 (MQ=255)                                                                                                                                                                          gTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCGTTTGCCTTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACAGGCACACTATTCAGGTTTCCGATGACGCACGa                                                                   >  1:604507/1‑301 (MQ=255)                                                                                                                                                                                            ttttAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCGTCTGCCTGACAAAAAAAGAAAGTTCAGGAAAATTACAATCGCACTCATACACCATCAACGTGCTATCGCAGGCAAAAGATAGTATTGGAAAAAAATGCTCTGTTTTTCTAGTGGAATCGCTGGATTACCTTATCACAGATGGTGCCGATGACATCGGCAAGATCCTCAATGTTCGAACACTTCACCGCTACAACCTTTAGTTTGCCAATGACACACGAtgacttttcatgaggtga                                                 >  1:729204/1‑284 (MQ=255)                                                                                                                                                                                                       ccGGTGTATTCCCACTTGCGCTGATGCTTTCAGGCTGGATCTCGTATGCGCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGAAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGa                                                                                                                           <  2:971872/216‑1 (MQ=255)                                                                                                                                                                                                         ggTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTCTTTTCCCTTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCACTGATCCACCACTTAACCGGCACAATATTAAGGTTGCAGATGACGCAAGATGTGTTTTAAATGCGGGTAAAAAGAACTTTa                                    >  1:251132/1‑301 (MQ=255)                                                                                                                                                                                                                     cTGTCGCTGATGCTTTTAGGCTGTATCTTGTATGCTCTGGTTTCACATACCGCAAAGGGTAGTTCAGGAAAGTATCCATCGCAGCCAGACACCATTGGTGGGCGATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTGAACAGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAatat                        <  1:874518/301‑1 (MQ=255)                                                                                                                                                                                                                        gcgcTGATGCTTGCAGGCTGTATCTCATATGCTCTGGTTACCCATAACGCAAAGGGTAGTTCAGGAAAGTATCAATCACAGTCAGACACCATCACGGGGCTATGGCAGGCAAAAGATAGTAAGGGAACAAAAGACTATATTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGGTCCAGCAc                                                                                                         <  2:60305/217‑1 (MQ=255)                                                                                                                                                                                                                         tggtGATGCTGTCGTGCTTTATCTAGGATGCTGTGTTTTCGCATACCGCAAAGGGTATTTCATGAAAGTATCAATAGCAGTCAGATATCATCACGGGGCTATGGCAGGCAAAAGAGAGTAATGTAACAAAAGGATATATTGTTGGAGGGGAATCGCTGGATTACCTTATCACTGATTGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCat                                                                                                         <  2:1011504/213‑2 (MQ=255)                                                                                                                                                                                                                              aTGATTTCAGGCTGGATATCGTATGCTCTGGATTACAATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGAAACCATCACTGGGCTAGCGCAGGCAAAAGATAGTAATGTAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAAc                                                                                                    <  2:1016673/216‑1 (MQ=255)                                                                                                                                                                                                                              aTGATTTCAAGCGGAATCTCTGATGCTATGGTTTCCCATACCGCAAGGGGTAGTTCGGGAAAGTATCAATCGCATTCAGAGACCATCACTGGGCTATCATATGCAAAAGTTATTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAAc                                                                                                    <  2:1003847/216‑1 (MQ=255)                                                                                                                                                                                                                                 cTTTCAGGCTGTATCTCGTATGCCGTGGTTTGCCTTAACACAAAGGGTAGTTCAGGAAAGTATCAATCCCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAACAGGCTACTCTTTTGTTGGGGAATCGCTGGATTACCTTATCACTGATGTTGCCGATGACACCGTTAAGATGCTCAATGATTCATCACTTAACCGGTACAATATTTCAGTTGCCGATGACGCACGACTTGATTTAAATGCTGGGAAAAAGAACTCTATCTGCACAACACCCCTTtacta              >  1:281522/1‑299 (MQ=255)                                                                                                                                                                                                                                 cTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAACTGCGGGGAAAAAGAAATTTACCGGCACACTATCGCTTtacta              >  1:320723/1‑299 (MQ=255)                                                                                                                                                                                                                                       ggCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACGGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAatan                                                                                         >  2:730845/1‑217 (MQ=255)                                                                                                                                                                                                                                        gCTGTAGCTCGTATGCTCTTGTTTCCCATACCGCAAAGGGTATTTCAGGAAAGTATCAATCGGAGTCAGACACCATCTCTGTGCTAGCGCAGGCAAAGGAGAGTAATGTAACAAAAGGCTATGTTTTTGTAGGGGAATCGATGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTATCCGGCACAata                                                                                          <  2:431051/216‑1 (MQ=255)                                                                                                                                                                                                                                         cTGTATCTCGTATTCTCTGGTTTCCCATACCGGAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGaa     <  1:730873/300‑1 (MQ=255)                                                                                                                                                                                                                                           gTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAg                                                                                     >  2:547348/1‑218 (MQ=255)                                                                                                                                                                                                                                            tATCTAGGATGCGCTGGGTTCCCATACCGCAAAGGGTAGTTAAGGAAAGTATGAATCGCAGTCAGACACCATCACGGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGGAGGGGAATCGCTGGATTACCTGATCACGGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATaa  <  1:820093/300‑1 (MQ=255)                                                                                                                                                                                                                                                    aTGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGAt                                                                            >  2:217164/1‑218 (MQ=255)                                                                                                                                                                                                                                                       ctctGGTTTCCCATACCGCAAAGGGTAGTTCAGGCAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAACAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGGTAAGATGCTCAATGATCCAGCACTTCACCGGCACGATATTCAGGTGGCCGATGa                                                                          >  2:84206/1‑217 (MQ=255)                                                                                                                                                                                                                                                       ctctGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGa                                                                          >  2:332736/1‑217 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                GTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAA  >  NC_000913/3875274‑3875808

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