Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,512,423	2 bp→TT	coding (283‑284/1239 nt)	mntH ←	Mn(2(+))/Fe(2(+)): H(+) symporter MntH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,512,423	0	C	T	95.6%	130.3 / ‑1.7	46	W95K (TGG→AAG)	mntH	Mn(2(+))/Fe(2(+)): H(+) symporter MntH
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (24/19);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-01
*	NC_000913	2,512,424	0	A	T	93.5%	130.5 / ‑0.3	46	W95K (TGG→AAG)	mntH	Mn(2(+))/Fe(2(+)): H(+) symporter MntH
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base T (23/20);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.88e-01

TGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACCACACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCT  >  NC_000913/2512180‑2512655                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                         tGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGACCGCGCCGTTGCAGCATTAAAATCAGGCAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCAGCACCAATAAATTCCGCCAGGTAGGCTGCCATCGCAATAATTTCTGCCTGAAACCAATAGAACTTCACTAc                                                                                                                                                                                                                                    >  1:295393/1‑251 (MQ=255)  gggAGAAAATCAACCCGACAATGCAAGCCGCGGCAACAAACAACATTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTAc                                                                                                                                                                                                                                    <  1:308702/250‑1 (MQ=255)   ggAGAAAATCAACTCGACAATGTAAGCCGCGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACgg                                                                                                                                                                                                                                  >  1:357222/1‑251 (MQ=255)                    cAATGTAAGCCGCGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATcg                                                                                                                                                                                                                 <  1:12555/251‑1 (MQ=255)                     aaTGTAAGCCGCGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATcgc                                                                                                                                                                                                                >  2:274643/1‑251 (MQ=255)                      aTGTAAGCCGCGGCAACAAACAACAGTAACCCGCCAATCTCTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATcgcg                                                                                                                                                                                                               <  1:506108/251‑1 (MQ=255)                              cccgccaaaaaacaaCAGAAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAAgctgct                                                                                                                                                                                                       <  2:30696/242‑1 (MQ=255)                                gggCAACAAACAACGGCAACCCGCCATGCACTTTCTCCAGCGGTTTTTGCCCGCGACGGTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCTCCCTGCAACAACGAAACACCAAGAATGAGTTTAAATCCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAAGAGAACTTCACTACGGGACGCGGATAGTGATCGCGAAGCTGCTcc                                                                                                                                                                                                     <  2:498374/250‑1 (MQ=255)                                 ggCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGACCCCAATATAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCg                                                                                                                                                                                                    >  2:362378/1‑251 (MQ=255)                                                            cACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCt                                                                                                                                                                          <  1:232788/250‑1 (MQ=255)                                                                  ctcCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTgg                                                                                                                                                                   >  1:157946/1‑251 (MQ=255)                                                                           tttttGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAgg                                                                                                                                                           >  1:231326/1‑250 (MQ=255)                                                                                 cccGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGa                                                                                                                                                     >  1:63618/1‑250 (MQ=255)                                                                                    gcgactcttcagcaTCAAAAACAGGAAAGTCGCGATCCCAGTCAGCGTCGCGCCCTGCAACAACGAAACACTAAGAATGAGTTTAAAGCCGAACGCCGCACCAACAAATTCCGCCAGGGCGGTTGCCATCGCAATAATTGCGGCCTGAACCCAATAGAACTTCACTACGTGACGCGGATAGTGATCGCGAATCTGCTCCGGCAGATTCTTACCGGTGGCAATCCCTAGTTTGGCAGAGATGATCTGACTCa                                                                                                                                                 <  2:126701/242‑1 (MQ=255)                                                                                        cGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCAt                                                                                                                                             >  2:195144/1‑251 (MQ=255)                                                                                                    aaaaTCAGGAAAGGCGCGAGCCCCGTCAGCACCGCGCCCTGCAACAACGAAACTCCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTTCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCAGCAGGt                                                                                                                                  <  1:107572/250‑1 (MQ=255)                                                                                                    aaaaTCAGGAAAGACGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGACACACCAAGAATGAGTTTACAACCGATCGCCGCACCAATAAATTCCGCCCGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGTACTTCACTAAGGGCCGCGGATAGTGATAGAGAATCTGTTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAAtaggcatagcactcagttt                                                                                                                                 >  2:433205/1‑233 (MQ=255)                                                                                                              aaGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGATTCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACg                                                                                                                       >  2:376991/1‑251 (MQ=255)                                                                                                              aaGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCAATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATAGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACg                                                                                                                       >  1:233943/1‑251 (MQ=255)                                                                                                                      aTCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAAccc                                                                                                                >  1:213567/1‑250 (MQ=255)                                                                                                                           cGTCAGCACCGCGCCCTGCAACAACGCAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGAGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCCGTCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCCAACGACAACCCA‑CAGt                                                                                                          >  1:438377/1‑250 (MQ=255)                                                                                                                                  aCCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTTGCCCAAACGACAACCCA‑CAGTAGCTGATa                                                                                                  <  1:158124/251‑1 (MQ=255)                                                                                                                                    cgcgCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGc                                                                                                >  2:361867/1‑251 (MQ=255)                                                                                                                                         ccTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGa                                                                                             >  1:6197/1‑249 (MQ=255)                                                                                                                                             caacaaCGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCt                                                                                         >  1:269532/1‑249 (MQ=255)                                                                                                                                                caacGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCTAGCAc                                                                                    <  1:53407/251‑1 (MQ=255)                                                                                                                                                    gAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CATTAGCTGATAGCCGAAGCTAGCACCCGc                                                                                >  2:216522/1‑251 (MQ=255)                                                                                                                                                                 tGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTc                                                                   >  2:1401/1‑251 (MQ=255)                                                                                                                                                                 tGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCAAGCACCAGCCTGAATATTGGTc                                                                   >  2:293608/1‑251 (MQ=255)                                                                                                                                                                  gAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCg                                                                  <  1:361027/251‑1 (MQ=255)                                                                                                                                                                       tAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCaaa                                                              <  1:195144/250‑1 (MQ=255)                                                                                                                                                                          gacctaTCGCCGCACAAATAAATTAAGCCAGGCAGGATGCCATCACAATAATTTATGCCTGACCCCAATAGAACTACACTACGGGACGCGGATAGTGATCGCGAATCTGATCCGCCAGCTTTTTGCCGGTGGCCATACCTCGTTTGACAGAGAGGATCTGGATCAGCATCGCCATCAGGTTGGCCCCAACGA‑AACCAAACAGTAGCTGCTAGCCGAAGCTAGCTCCCGCCTggcaattgatcgcaaactt                                                           >  2:352678/6‑232 (MQ=255)                                                                                                                                                                             cGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACCACACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCg                                                       >  2:321131/1‑251 (MQ=255)                                                                                                                                                                                 cgccgcACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACCACACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGCATCTGAATCAGCATCGCCATCAGATTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGAt                                                   >  1:27410/1‑251 (MQ=255)                                                                                                                                                                                 cgccgcACCAATAAATTCCGCCAGGTCGGTTGCCAGCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGAt                                                   <  2:189334/251‑1 (MQ=255)                                                                                                                                                                                            aaaaTTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCGATATAACCa                                        <  1:1401/250‑1 (MQ=255)                                                                                                                                                                                                      cAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCa                               >  2:14107/1‑250 (MQ=255)                                                                                                                                                                                                        ggTCGGTTGCAATCGCAATAATGGCCACCTGAACCCAAAAGAACTTCACTCCGGTACGCAGATAGGGAGCGCGAAACTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCACCGCCATCAGGTGGGCCCAAACGACAACCCA‑CAGTAGCTGAGAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCTATATAACCAATCGCCGCAAtt                            <  2:283120/251‑2 (MQ=255)                                                                                                                                                                                                           cGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAAc                         <  1:491948/251‑1 (MQ=255)                                                                                                                                                                                                                gCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAgg                    <  1:189797/251‑1 (MQ=255)                                                                                                                                                                                                                       cAATAATTTCTTCCTGAACCCAATAGAACTTCACTACGGGACGCGAATAGTGACAGCGAATCTGCCCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCAGCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCTAGCCCCCGCCTGAATAGTGGTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGGCCCAtt             <  2:219084/251‑1 (MQ=255)                                                                                                                                                                                                                            atTTCTGCCTGAACCCACTAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGAGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCACACGACAACCCA‑CAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGc        >  2:30420/1‑251 (MQ=255)                                                                                                                                                                                                                            atTTCTGCCTGAACCCAATAGAACTTCACTACGGGAAGCGGATAGTGATCGCGACTCTGCTCCGCCAGATTTTTACCGGAGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCTAGCACCCGCCTAAATATTGGTAGCAAAGTTAACGAGATCGATATAACCAATCGCCGCAATGTACGCCGGACACATTAATGc        >  1:152807/1‑251 (MQ=255)                                                                                                                                                                                                                                  gCCTGAACCCAATAGAACATCACTACGGGACGCGGATAGTGATCGCGAATCGGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAGCCCA‑CAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTATTGCGAGCCt  <  1:493708/251‑1 (MQ=255)                                                                                                                                                                                                                                     tGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCca                                                                                                               <  2:23943/140‑1 (MQ=255)                                                                                                                                                                                                                                     tGAACCCAATAGAACTTCACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCca                                                                                                               >  1:23943/1‑140 (MQ=255)                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                         TGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACCACACTACGGGACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCA‑CAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCT  >  NC_000913/2512180‑2512655

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