Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,102,348	G→A	L169L (CTG→TTG)	trmB ←	tRNA m(7)G46 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,102,348	0	G	A	97.5%	121.3 / NA	40	L169L (CTG→TTG)	trmB	tRNA m(7)G46 methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/14);  minor base C (1/0);  total (26/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

TTTACGCAGACGACGGCTACGGTTCTTTGCCATTATTTCACCCTCTCGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGA  >  NC_000913/3102101‑3102581                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                          tttACGCAGACGACGGCTACGGTTCTTTGCCATTATTTCACCCTCTCGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAAttt                                                                                                                                                                                                                                        >  1:141383/1‑251 (MQ=255)             acgGCTACGGTTCTTTGCCATTATTTCACCCTCTCGAACATTAAGTCACATACACCGTGACCAAGACGATGACCACGTTGTTCAACTTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTAcc                                                                                                                                                                                                                              >  1:23138/1‑249 (MQ=255)                     ggTTCTTTGCCATTATTTCACCCTCTCGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCCAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCCGCTGCAATTTGCTTTTTACCAGTTCGGCa                                                                                                                                                                                                                     >  1:362265/1‑250 (MQ=255)                      gTTCTTTGCCATTATTTCACCCTCTCGAACATTAAGTCCCATACTCCGTGACCAAGACGATGAACACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTCCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGcca                                                                                                                                                                                                                    >  1:218355/1‑248 (MQ=255)                           ttGCCATTATTTCACCCGCTCGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACgg                                                                                                                                                                                                                <  1:302720/249‑1 (MQ=255)                            tGCCATTATTTCACCCTCTAGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTAAAATTTCGTAACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATAACTTAAAGCATATGTTCCGAATAAGGTTCCCAGTCGGTCGCCATATGGAATACGACACCCAGATGCAATTTGATTTTTACCAGTTCGGCAAACGGCAc                                                                                                                                                                                                             >  2:386888/1‑251 (MQ=255)                             gCCATTATTTCACCCTCTCGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGTTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCAcc                                                                                                                                                                                                            >  1:273465/1‑251 (MQ=255)                                       cACCCTCTCGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATAc                                                                                                                                                                                                  <  2:302394/251‑1 (MQ=255)                                               cGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGtt                                                                                                                                                                                            <  2:140868/249‑1 (MQ=255)                                                      tAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATgg                                                                                                                                                                                    <  1:390807/250‑1 (MQ=255)                                                      tAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGgc                                                                                                                                                                                   >  1:352805/1‑251 (MQ=255)                                                       aaGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACAGCGTTTATTATGgcg                                                                                                                                                                                  >  1:202060/1‑251 (MQ=255)                                                                 cTCCGTGACCAAGAAGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGATACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGc                                                                                                                                                                        >  2:285244/1‑251 (MQ=255)                                                                         ccAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTc                                                                                                                                                                <  2:156150/251‑1 (MQ=255)                                                                          cAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCttcccgccggtgaccg                                                                                                                                                               >  2:316296/1‑237 (MQ=255)                                                                              accaTGACCACCTGGCTCAAATAACGTCACCGGACGTTATGACGGCAGCGGTCCGAAATCATTGCTCTCTGACAGGTTTTTATAACAGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGACGGTCGACATATGGCATACGCCCCCAAGCTGCTACTTGCTTTTAACCAGTTTGGCCAAGGGCAACTGACCGAAACGGCGTTTATTCTGCAGCGCTTTGTGCCACTGGTCAggg                                                                                                                                                            >  2:229644/4‑250 (MQ=255)                                                                                  tGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGaaaa                                                                                                                                                        >  1:385157/1‑250 (MQ=255)                                                                                  tGACCACGTTGTTCAAATTTCGTCACCGGACGTGAGGGCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTGGTACCAAGTGGCAGGGaaaa                                                                                                                                                        >  2:139200/1‑250 (MQ=255)                                                                                      cACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGc                                                                                                                                                   >  1:354967/1‑251 (MQ=255)                                                                                                       gTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCCATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATg                                                                                                                                   >  1:92752/1‑250 (MQ=255)                                                                                                           ccGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTg                                                                                                                              >  2:1085/1‑251 (MQ=255)                                                                                                                     tGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATc                                                                                                                     >  2:27097/1‑250 (MQ=255)                                                                                                                          gACGCGGTGCGTAATCATTGCTCGCTGACAGGTTTTTATAACCGGCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTGTTTGCCAGTTCGGCAGACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTa                                                                                                               <  1:357523/251‑1 (MQ=255)                                                                                                                             gcgGTACGTAATCATTGCTCTCTGACAGGGTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGTTAAAAAGAGCTGCACCATGCGCAATGAATTGGCAGGAAGCATTTTATgc                                                                                                            <  2:371584/251‑1 (MQ=255)                                                                                                                                         cATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAAcc                                                                                                >  2:212752/1‑251 (MQ=255)                                                                                                                                         cATTGCTCTCTGACAGGAATTTATAACCGTCAATAGAGGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCCGCTGCACTTTGCTTTTTACCAGATCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGAGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGAACCATGCGAAATGAATTGTCAGGAATCATTTTATACAGAACTTCAAcc                                                                                                >  2:387265/1‑251 (MQ=255)                                                                                                                                            tgttctctGACAGGTTTTTATAACGGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTCTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAATAGCTGCTCCATGCGCAATGAATTGTCAGTAATCATTTTATGCAGCACTTCAACCGc                                                                                              <  2:388907/247‑1 (MQ=255)                                                                                                                                             gCTCTTTGACCGGGTTTTATAACCGTAAATAGAAGACATCACTTCAAGCATAGGTTCCGCATAAGGCTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTAGGCACACGGCACCTGAACGATACGGCGTGTATTATGGCGAGCTTTGGGACACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGCATCATTTTATgctggcatccaaccgaac                                                                                            >  2:49881/1‑235 (MQ=255)                                                                                                                                                            tttATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCAcgc                                                                             >  1:336199/1‑251 (MQ=255)                                                                                                                                                                       cAATAGAAGACATCACTTCAAGCATATGTTACGCATAAGGTTCCCAGTCGGTCGCAATATGGAATACGCCCTCCAGCTGCAATTTGCTTTATACTAGTTCGGCAAACGGCACATGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCAATTTATGCAGAACTTCAACCGAATCGTGACACATCCAGTGCAGGTTGCtt                                                                  >  1:223243/1‑251 (MQ=255)                                                                                                                                                                        aaTAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTa                                                                 <  1:12334/251‑1 (MQ=255)                                                                                                                                                                                 cATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCAATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTATCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACAGCATCGTGACACATCACGCGCAg                                                                         >  2:255902/1‑234 (MQ=255)                                                                                                                                                                                 cATCACTTCAAGCATATGTTCCGCATAAGGTCCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCTCCATGCGCAATGAATTGTCAGGAATCATTTTATGCACCACTTCAACCGCATCGTGACACATCACGCGCAg                                                                         <  1:255902/234‑1 (MQ=255)                                                                                                                                                                                             cATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAg                                            <  1:269983/251‑1 (MQ=255)                                                                                                                                                                                                              gTNCCCAGTCGGTCGTCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTCCCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGCTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCCGGTTGCTTAAGCCTTCTTCACGCGCAGAAGCCAGGCACCCGCCAAc                            <  2:354967/250‑1 (MQ=255)                                                                                                                                                                                                                  ccAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCTCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCg                       <  1:220025/251‑1 (MQ=255)                                                                                                                                                                                                                    aGTCGGTCGCCAGATGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGt                     <  2:385157/251‑1 (MQ=255)                                                                                                                                                                                                                                tatGGAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGGGAATGCACTTc          >  2:299672/1‑250 (MQ=255)                                                                                                                                                                                                                                    gAATACGCCCCCCAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGGCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGc     >  1:286787/1‑251 (MQ=255)                                                                                                                                                                                                                                       tacgcacaccAGCTGCAATTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTACCACGGGTCACGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCCTGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGNCgc  >  2:356472/9‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                          TTTACGCAGACGACGGCTACGGTTCTTTGCCATTATTTCACCCTCTCGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGA  >  NC_000913/3102101‑3102581

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑