Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,475,333	0	C	A	47.1%	5.8 / 74.9	51	*129L (TGA→TTA)	tusD	sulfurtransferase for 2‑thiolation step of mnm(5)‑s(2)U34‑tRNA synthesis
Reads supporting (aligned to +/- strand):  ref base C (14/13);  new base A (14/10);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01

CGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGAAGCCGCACAAACCCAGCACTGTTCAATGTCGTACAGACCCAACAATTTAAAAGTGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCTGCTTAACTC  >  NC_000913/3475044‑3475622                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                  cGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGAAGCCGCACAAACCCAGCACTGTTCAATGTCGTACAGACCCAACAATTTAAAAGTGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAACCGCACTTCGTTTCATCAGAACTGTac                                                                                                                                                                                                                                                                                         >  1:175498/1‑300 (MQ=255)    caacaaAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGAAGCCGCACAAACCCAGCACTGTTCAATGTCGTACAGACCCAACAATTTAAAAGTGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTaccac                                                                                                                                                                                                                                                                                      >  2:400582/1‑300 (MQ=255)         aaGGGTATCTGCGGATCTAACCCGCGTTCGCGCAGTGAAGCCGCACAAACCCAGCACTGTTCAATGTCGTACAGACCCAACAATTTAAAAGTGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTc                                                                                                                                                                                                                                                                                <  1:367077/301‑1 (MQ=255)               gtCTGCGGATCTAACCCGCGTTCGCGCAGTGAAGCCGCACAAACCCAGCACTGTTCAATGTCGTACAGACCCAACAATTTAAAAGTGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGt                                                                                                                                                                                                                                                                          >  1:338772/1‑301 (MQ=255)                       aTCTAACCCGCGTTCGCGCAGTGAAGCCGCACAAACCCAGCACTGTTCAATGTCGTACAGACCCAACAATTTAAAAGTGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAgg                                                                                                                                                                                                                                                                  >  2:399746/1‑301 (MQ=255)                                  tgtgctcttccgatctcGCACAAACCCAGCACTGTTCAATGTCGTACAGACCCAACAATTTAAAAGTGGCAATGTAATCACGCGCCAGCGCTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAg                                                                                                                                                                                                                                                        <  2:338989/284‑1 (MQ=255)                                                  cGCACAAACCCAGCACTGTTCAATGTCGTACAGACCCAACAATTTAAAAGTGGCAATGTAATCACGCGCCAGCGCTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAgagatcggaagcgcgtcg                                                                                                                                                                                                                                       >  1:338614/1‑284 (MQ=255)                                                  cGCACAAACCCAGCACTGTTCAATGTCGTACAGACCCAACAATTTAAAAGTGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCtt                                                                                                                                                                                                                                        >  1:124882/1‑300 (MQ=255)                                                          cccAGCACTGTTCAATGTCGTACAGACCCAACAATTTAAAAGTGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCCATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCCAGTCCGCTTCAGGTaaaa                                                                                                                                                                                                                               >  1:356877/1‑300 (MQ=255)                                                           ccAGCACTGTTCAATGTCGTACAGACCCAACAATTTAAAAGTGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAAt                                                                                                                                                                                                                               >  2:372061/1‑300 (MQ=255)                                                           ccAGCACTGTTCAATGTCGTACAGACCCAACAATTTAAAAGTGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATc                                                                                                                                                                                                                              <  2:367311/301‑1 (MQ=255)                                                           acgctcttccgatctggcctACAGACCCAACAATTAAAAAGTGGCAATGTAAGCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAAt                                                                                                                                                                                                                               <  1:44967/281‑1 (MQ=255)                                                                          gTCGTACAGACCCAACAATTTAAAAGTGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAAtagatcggaagagcac                                                                                                                                                                                                                >  2:45019/1‑285 (MQ=255)                                                                              tACAGACCCAACAATTTAAAAGTGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACg                                                                                                                                                                                                            <  2:125003/300‑1 (MQ=255)                                                                                    cccAACAATTTAAAAGTGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAagccag                                                                                                                                                                                                     >  2:263976/1‑301 (MQ=255)                                                                                           aTTTAAAAGTGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGACGCCAGCCCCAg                                                                                                                                                                                               >  1:125564/1‑300 (MQ=255)                                                                                              tAAAAGTGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCtt                                                                                                                                                                                           >  1:94226/1‑301 (MQ=255)                                                                                                   gtgtCAATGTAAGCACGCGCCAGCACGGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAAGGCGGAAGTTGCCAGTAAAGCATCTAAACCTGCCCGGCCTTCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCgg                                                                                                                                                                                       <  1:59844/300‑1 (MQ=255)                                                                                                    tGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGc                                                                                                                                                                                      <  2:125685/300‑1 (MQ=255)                                                                                                          tGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGtt                                                                                                                                                                               <  2:235411/301‑1 (MQ=255)                                                                                                            tAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGTCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTtc                                                                                                                                                                             <  2:358893/301‑1 (MQ=255)                                                                                                              aTCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTtcat                                                                                                                                                                           <  2:339147/301‑1 (MQ=255)                                                                                                                              gCATCGGGCTTTTGTCCTTGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGGCAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCCAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCCGCCCCAGTCTTCCGGCCTCCGTTTCATCCACCACGCCCCgcc                                                                                                                                                            >  1:336267/1‑298 (MQ=255)                                                                                                                              gCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCg                                                                                                                                                           >  2:283387/1‑301 (MQ=255)                                                                                                                               acgctcgtcagatctcTGGCAGCAGCAGAAAAAAGCCATCAGCTATAAAGAAGACCGCCAGATCGTCAGGTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGACCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCg                                                                                                                                                           <  1:103637/285‑1 (MQ=255)                                                                                                                                             ccTGGCAGCAGCTGAAAAACGCCATAAGCTATAAAGAAGACAGCCAGATCGTCAGTTAAGGCGGAAGTTGCCAGGAAAGCAACGAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTAc                                                                                                                                            <  1:399288/301‑1 (MQ=255)                                                                                                                                              cTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCgagatcggaagagcaca                                                                                                                                           >  2:103744/1‑285 (MQ=255)                                                                                                                                                 gcagcaGCTGAAAAACGCCATCAGTTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCa                                                                                                                                         <  1:283087/300‑1 (MQ=255)                                                                                                                                                       gcTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATAtt                                                                                                                                   >  2:298047/1‑300 (MQ=255)                                                                                                                                                        cTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTCAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCCGATATTc                                                                                                                                  >  2:254058/1‑300 (MQ=255)                                                                                                                                                        cTGAAAAACGCCATCAGCTATAAAGAAGAAAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACAGTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGatat                                                                                                                                    <  1:400124/298‑1 (MQ=255)                                                                                                                                                           aaaaaCGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAgcg                                                                                                                              >  2:300769/1‑301 (MQ=255)                                                                                                                                                             ccccctactccccgctcttccgatctgCCAGCTCTTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGCACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCCGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGcc                                                                                                                            <  1:99388/275‑1 (MQ=255)                                                                                                                                                             aaaaGCCACCAGCTGTAAAGAAGACAGCCAGCTCGTCAGTTAATGCGGAAGTTTCACCTAAAGCACCTAAACCTTCCCGGCCTGCTGCTGTACCATGAGGTGCAGTAGACAAAACAAACGCAATTCGTTTCATCATAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGcc                                                                                                                            <  1:131408/301‑1 (MQ=255)                                                                                                                                                               aCGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCcac                                                                                                                          <  1:371648/301‑1 (MQ=255)                                                                                                                                                                                       gCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCcagatcggaagagcacacgtctgaa                                                                                                    >  2:99487/1‑276 (MQ=255)                                                                                                                                                                                              cgcaaGTGAATGCGGAAGTAGCCAGAAAAGCATCTAAACCCTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGcc                                                                                           <  1:297727/297‑1 (MQ=255)                                                                                                                                                                                                                 tGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCAGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACCAAGATATTCAGCGCCACACCATGTTGCGCATTCCGGTGTTGCCAGGCCCGTACGAGGTCAAATTCGt                                                                        >  1:242725/1‑301 (MQ=255)                                                                                                                                                                                                                                tAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAgg                                                         >  2:352649/1‑301 (MQ=255)                                                                                                                                                                                                                                tAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAAGTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCCCTTGCCGGAGAgg                                                         >  1:20771/1‑301 (MQ=255)                                                                                                                                                                                                                                      ttCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCCTTCAGTTGTTGCCAGGCCCGATCGAGGTCAAATTCGTCACTTGCCGGAGAGGTCaagt                                                   >  2:309086/1‑299 (MQ=255)                                                                                                                                                                                                                                          cGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGtt                                               >  2:322780/1‑301 (MQ=255)                                                                                                                                                                                                                                                         ccATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATCGAccc                                 >  1:260520/1‑300 (MQ=255)                                                                                                                                                                                                                                                              aGGTGCAGTAGAAAAAACAACCGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTcc                            <  2:357277/300‑1 (MQ=255)                                                                                                                                                                                                                                                               ggTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTccc                           <  2:94325/300‑1 (MQ=255)                                                                                                                                                                                                                                                               ggTGCAGTAGAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGCTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCCATTGGTTAGCGTTATCGACCCCTTccc                           >  1:165214/1‑300 (MQ=255)                                                                                                                                                                                                                                                                        gAAAAAACAAACGCAATTCGTTTCATAAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGaaaaa                 >  2:398736/1‑301 (MQ=255)                                                                                                                                                                                                                                                                           aaaaCAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTTACTTTCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAATGa               >  2:336/1‑300 (MQ=255)                                                                                                                                                                                                                                                                            aaacacacGCAATTCGTTTCACAAGACCTGTACCACCCTGTCACAGGTCAGCGAAGCTTCCGCCCACGCCCCAAGTCCGCTTAATGTAAATCCCTGCGGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCAAGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTAGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACg             <  1:253797/301‑1 (MQ=255)                                                                                                                                                                                                                                                                            aaacaaagCCACGTCGTTTCAGAAGACCAGCGCCACCCTGTCACAGGGCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTACGGTACATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTACCCGCCTCCGTTTCGTCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACg             <  1:308745/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                       aTTCGTTTCATCAAAACTGACCCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTGGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCTGCTTAACTc  <  1:263701/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                  CGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGAAGCCGCACAAACCCAGCACTGTTCAATGTCGTACAGACCCAACAATTTAAAAGTGGCAATGTAATCACGCGCCAGCACTGCATCGGGCTTTTGTCCTGGCAGCAGCTGAAAAACGCCATCAGCTATAAAGAAGACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCGTAATGCTGCCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCTGCTTAACTC  >  NC_000913/3475044‑3475622

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