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.4%	78.9 / 12.1	51	A28A (GCT→GCC)	yidX	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (16/23);  new base C (10/0);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 7.86e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

CTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAG  >  NC_000913/3875224‑3875816                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                        cTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCtgtg                                                                                                                                                                                                                                                                                                      >  2:516779/1‑301 (MQ=255)               cGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATAcc                                                                                                                                                                                                                                                                                          <  1:77096/299‑1 (MQ=255)                gTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAGAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCa                                                                                                                                                                                                                                                                                       >  2:379858/1‑301 (MQ=255)                    tcacgtgtgagTTTGTCGTGCGTGTAATGTAGTATAAACTTATATTGGTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGAAGAAATCTGAAAAGGGAGATCAGGCGCTAAACCGGAGGGAAAGAAGAGTATGCTTTCAACTGATTAGCTACTCGTTTAAAGGATTAATCATGAAGTGTAATTGTAAGGGATTTTTTAAGGCGGCCGGTTTATACCCACTGGCGCGGATGCTTTCAGGCTTTATCTCGTATGCTCTGGTTTCCCATACCGCAAAgg                                                                                                                                                                                                                                                                                   <  1:450733/291‑1 (MQ=255)                        gTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGTTTTCCCTTACAGCAAAGGTTAg                                                                                                                                                                                                                                                                               >  2:182808/1‑301 (MQ=255)                              agttGCGGTGCGGTGAAAGTAGTACAAAGTTATATTGTTGTAATACAATTTAGATCACAAAAAGAACAATGCATAAAAGAAAGGCATACGACAGAGAGAAATCTGGAAAAGGATATAAGGAGATAAAGAAGAGGAAAAGAAGAGTATGATGTCAAGGACTTAGCTACACGTTTAAAGGATAAATCATGAAGTGGAATTTTAAGGGATGTATTAAGGCGGCCGGTTTATTCCCACTGGCGCTGATGCATTCAGGCTGTAGCTCGTATGCTCTGGTATCACATGCCGCAAAGGGTAGTTCAg                                                                                                                                                                                                                                                                          <  1:625435/298‑1 (MQ=255)                                ttGCCGAGCGTGTAATTTAGTACAAACTTATATTGTTGTACTAAAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAAAAATCTGAAAAGGGATATCAGGCTCGAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAAACATGAAGTTGAATTTTAAGGTATTTTTTAAGGCTGCCGGTGTATTCCCACTGGCGCTGATTCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCGAAGGGTAGTTCAGGaa                                                                                                                                                                                                                                                                       <  1:485075/301‑1 (MQ=255)                                   ccGGGCGTGTAATTTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAg                                                                                                                                                                                                                                                                     <  1:485531/300‑1 (MQ=255)                                       gcgtgTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCTTACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCGTTTGCTAGACAGAAAACGTTAGTTCAGTCAAGTATc                                                                                                                                                                                                                                                                 >  2:395956/1‑300 (MQ=255)                                         gtgtAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTTTCCCTTACCGAAAAGGGTAGTTCAGGAAAGTATCAAt                                                                                                                                                                                                                                                              >  2:419092/1‑301 (MQ=255)                                          tgtAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCGCTGGTTTCCCATACCGCAATAGGTAGTTCAGGTAAGTATCAAt                                                                                                                                                                                                                                                              >  2:445895/1‑300 (MQ=255)                                              aTGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAg                                                                                                                                                                                                                                                         >  2:133297/1‑301 (MQ=255)                                                   gtaCAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGACGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTGGTTTGCCTTAACGCAAAGAGTAGTTCAGGACAGTACCAATCGCATTCaaa                                                                                                                                                                                                                                                    >  2:629730/1‑299 (MQ=255)                                                   gtaCAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCCTGAAAAAAAAGAGTTGTTTAGGAAAGTAGTAATCGCAGTCCGa                                                                                                                                                                                                                                                    >  2:633227/1‑301 (MQ=255)                                                         aataacatTGTAGTACTATAATTTAAATCACAAAAAGAACAAGTCATAAAAGAATGTAAGGCGGTACGCAGAAATATTGTAAGAGATAACAGTCGATAAACAGGAGAGAAAGAAAAGTATGAAGTCAACGGATGAGCTACACGTTTAAAGGAGTAATCATGAAGTTAAATTTTAAGGGATTTTTGAAGGCTGCCGGTTTATTCCCACTGGCGCTTATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACAcc                                                                                                                                                                                                                                                <  2:521475/293‑1 (MQ=255)                                                              tatagaaagactacAATTTAGATCACAAGAAGGACAATGCATAAGAAAATGATATGTGTCGGGCAGAAAACTGTTAAAGTATTGCAGGAGATAAACAAGAGGGAAAGAAGAATAGGCTATCAACGGCTTAGCTACTCGTTTAAAGGATTAATAATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTg                                                                                                                                                                                                                                         <  2:280506/292‑1 (MQ=255)                                                                          acAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGTCTATCGCAgg                                                                                                                                                                                                                             >  2:26541/1‑301 (MQ=255)                                                                          acAATTTAGATCACAAAAAGAACAATGAATAAAAAAATGACATGCGTCGGGGAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAgg                                                                                                                                                                                                                             <  2:368315/301‑1 (MQ=255)                                                                                tAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCCTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAg                                                                                                                                                                                                                       <  2:108890/301‑1 (MQ=255)                                                                                  gATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTTTGCTTTAACACAAAGAGTAATTAAGGAAAGTATCCATCGCAGTCAGACACCATCACTGGGCTAACGCAGGCAAAAGAt                                                                                                                                                                                                                     >  2:572470/1‑301 (MQ=255)                                                                                     cacaAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTTAAAGCAAAAGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGTCTATCGCAGGCATAAGATAga                                                                                                                                                                                                                  >  2:222828/1‑300 (MQ=255)                                                                                         aaaaGAACAATGCATAAAAAAATGACATGCGTCGGGCAGATATCTGGAAAGGGATATCAGGCGCTAAACAGGAGGAAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGAGTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATg                                                                                                                                                                                                              <  2:316576/301‑1 (MQ=255)                                                                                               aCAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACaa                                                                                                                                                                                                        >  2:236759/1‑301 (MQ=255)                                                                                                       tAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTAt                                                                                                                                                                                                >  2:373631/1‑301 (MQ=255)                                                                                                         aaaaaaTGACATGCGTCGGGCAGAAATCTGTAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGA‑TTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTCTCCCATACTTCAAAGCGTAGTTAAGGATAGTATCAATCGCAGACCGACACCTTCACTGGGCTATCACAGTCAATAGCTATGAATAGTACAAAACGCTATGt                                                                                                                                                                                              >  2:293421/1‑300 (MQ=255)                                                                                                         aaaaaaTGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGt                                                                                                                                                                                              >  2:233780/1‑301 (MQ=255)                                                                                                              aTGACATTCGGCGGGCAGAAATCTGAAAAGGGATATCAGGAGCTAAACAGGAGAGAGAGAAGAGTATGCGTTCAACGGCTTAACGACTTGTTTAAAGTATTAATCATGAAGTTGAATTATAAGGGTTTTTTTAAGGCTACTGGTGTACTCCCACTGGCGCTGATGCTTTCAGGATGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAGAAGATAGTAATGGAACAAAAGGCTATGTTTTTg                                                                                                                                                                                         <  1:379853/301‑1 (MQ=255)                                                                                                                       gTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTCCTTGAAAAAAAAAAGTAGTTCAGCAAAGTATCAATAGCAGTCAGACACCATCACTGTGCTCTCGCAGGCAACAGATAGTAATTGAACAAAAGGCTATGTTTTTGTAGGGGAAt                                                                                                                                                                                >  2:668227/1‑301 (MQ=255)                                                                                                                                         aaGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGTTGGATTACCTTATc                                                                                                                                                               >  1:355347/1‑300 (MQ=255)                                                                                                                                            ggagagaaGGCGCGAGAAAGGAGGGAAAGAAGAGTATGCTTTTAAGGGAGTAGATACTCGTGTAAAGGATTAATCATGAAGTTGAATATTAAGGGATTTTTGAAGGATGCCGGTTTATTCACACTGGCACTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTGTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCAc                                                                                                                                                             <  1:588830/292‑1 (MQ=255)                                                                                                                                             tatatCACGCGATTATCAGGAGGGTAAGAAGAGTGTGCTTTCTGCGGCTTAGCTATTAGTTTAAAGTATTAATCATGAAGTTGAATTATAAGTGATTTTTTAAGGCTGCCGATTTATTCCCAATGGAGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATAACGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTTAGACACCATCACTGTGATATCGCAGGCAAAAGTTAGTAATGTAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTg                                                                                                                                                           <  1:163691/299‑1 (MQ=255)                                                                                                                                                 tCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTTCCATACCGCAAAGGCTAGATCAGGAAAGTCTCAATTGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAATGATAGTACTGGAACAAAAGGCTATGTTTTTGTAGGGGAATCTCTGGATCACCTTATCACGGATGGt                                                                                                                                                      >  2:115742/1‑301 (MQ=255)                                                                                                                                                    ggCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAATGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGAAACCACCACTGGGCTATCGCAGGCAAAAGATAGTGAGGGAACAAACGGCTATGATTTTGTAGGGGAATCGCTGGATTACCTTATCAATGATGGTGc                                                                                                                                                    >  1:283828/1‑300 (MQ=255)                                                                                                                                                                   ggAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAATGCTGCCGGTTTATTCCCACTGGGGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTaa                                                                                                                                     <  1:236755/300‑1 (MQ=255)                                                                                                                                                                            agTATTCTTTCAACGGCTTAGCTACTCGTTTAATGGATTAATCATGAAGTTGAATTTTAATGGATTTTTGTATGCTACCGGTTTCTTCCCACTGTCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCa                                                                                                                             <  2:283832/299‑1 (MQ=255)                                                                                                                                                                                  ccttCGACGGCTTAGCTAATCGTTTAAAGGAGTAATCATGAAGTTGAATTTTGGGGGATTTTTTAAGGCTAACAGGTTATGCCCACTGGCGCTGATGCTATCAGGCAAGAACTCGAATAAACTAATTGCCCATAACGAGAAGGGTAATTCAGAAAAGTATGAATCACAGTCGGACACGAACAATGAGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTGATCACTGATGGTGCCGATGGCAGCGTTAAGATGCTCAATGATcc                                                                                                                      <  1:115741/298‑1 (MQ=255)                                                                                                                                                                                            tttAGCTACTCGTTTAAACGATTACTCATGAAGTTGAATTTTAATGGAATTTTTAAGGCTGCCAGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAAcc                                                                                                           <  1:233776/300‑1 (MQ=255)                                                                                                                                                                                                   aCTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACaa                                                                                                    >  1:1239/1‑301 (MQ=255)                                                                                                                                                                                                    cTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGCTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTCGTATGCCTTACCGAAAAGAGCAGTTCAAGACATGATCAATCGCAGTTAGTCTCCATCTCTTGGCTATCGCCGGCAAAATAAAGTATTTTAATAAAAGTCTATGTTTTCGTAGTGGAATGGTCGGCTTTCCCTATTACTACTGTTTCCGCTGACATCGTCACTACGCCCAATCATCTCGCACTTAACGGTcacac                                                                                                    >  2:94508/1‑299 (MQ=255)                                                                                                                                                                                                      cGTTTAAAGGATTAATCATGAAGTTGAATCTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGTTAGTTCAGGAACGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGCACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGTTCAATGATCCAGCACTTAACCGGCACAatat                                                                                                 >  1:417649/1‑301 (MQ=255)                                                                                                                                                                                                       gTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTCGTTTGCCTTACCAAAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCCCAATAtt                                                                                                >  2:556564/1‑301 (MQ=255)                                                                                                                                                                                                         ttAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCa                                                                                              >  2:153473/1‑301 (MQ=255)                                                                                                                                                                                                         ttAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCa                                                                                              <  2:1239/301‑1 (MQ=255)                                                                                                                                                                                                         ttAAAGGATTAAGCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATAACGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCa                                                                                              <  1:133296/301‑1 (MQ=255)                                                                                                                                                                                                                               gAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCACGTTGCCGATGACGCAAGCttt                                                                        >  2:121507/1‑301 (MQ=255)                                                                                                                                                                                                                                 aTTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCCTGAAAAAAAAAGGTAGTTCAGGAAAGTACCAATCGCAGTCAGACACCATCACTGGGCTATCGCCGGCACAATATAGTAATGGAATAACAGGCTCTGTTTCTTTAGGGGAATCCTTGGCTTACTTTATCACCGATGGTTCCGTTCACTTCGTTAATTTGTTCAATGTTCTTGCACTTACCCGGTACATTATTCAGGTTCCCgtatcctcaactgttgt                                                                      >  2:306786/1‑285 (MQ=255)                                                                                                                                                                                                                                      aaGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCTATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCCTGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACAGGCACAATATTCAGGTTGCAGATGACGCTAGATTTGTTTTaa                                                                 >  1:630922/1‑301 (MQ=255)                                                                                                                                                                                                                                                 tAAGGTTGCTGGTTTATTCCCACTGGCTCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTACCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGAATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGaaa                                                      <  1:373626/301‑1 (MQ=255)                                                                                                                                                                                                                                                      atgccagatgaTTCCCACGGGCGATTGTGCTCTCAGGCTGTCGATCGTATGATCTGGTTGTCAATAACGATAAGGATAGTCAAGAAAAGTCACAAACGCAATCAGAAAGCAACAGGGGGATATCGCAGGCAAGAGATAATAATGAAACAAAAGGCTATGTATGTGTAGGGGAAGCGCTGGATTACCTTAACACTGTTGTTGCCGTTGCCATCGTTAAGATACTTAATGATCAAGGACGTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGTATTGTGTTTAATGGGGGGAAAAAGaa                                                 <  1:222824/291‑1 (MQ=255)                                                                                                                                                                                                                                                                             cTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCttt                          <  1:153472/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                     ttCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTtactactac                  <  2:355352/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                     tATGCTCTTGGTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAAGCGCAGTCAGAAACCATCAATGTGCTATCGCAGGCAAAAGATAGTAATGAAACAAAAGGCGATGTTTTTGAAGGGGAATCGCTGGATTACCTTATCACTGATGGTGACGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCTGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCTGCACAATATCGCTTTACTACTACGGGAATAACgaagaag  <  1:26541/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                        CTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAG  >  NC_000913/3875224‑3875816

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