Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,512,423	C→T	W95* (TGG→TAG)	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	94.9%	102.5 / ‑1.8	40	W95* (TGG→TAG)	mntH	Mn(2(+))/Fe(2(+)): H(+) symporter MntH
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (20/17);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

AGGCTGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACG‑TTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACCACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCG  >  NC_000913/2512176‑2512664                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   aGGCTGGGAGAAAATCAACCCGACAATGTAAGCCGCGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACG‑TTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACtaca                                                                                                                                                                                                                                                 <  1:131868/251‑1 (MQ=255)                 aacccgaaaaTGTAAGCCGCGGCAACAAACAACAGTACCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACG‑TTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATa                                                                                                                                                                                                                                <  1:22785/243‑1 (MQ=255)                         aaggTAAGCCGCGGCAACAAAAAACAGAAACCCGCCAATCACTTTCTCCAGCGGTGTTTGCCCGCGACG‑TTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATcgc                                                                                                                                                                                                                        <  1:133560/248‑1 (MQ=255)                              aaGCCGCGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACG‑TTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAAt                                                                                                                                                                                                                    <  1:609323/250‑1 (MQ=255)                                  cgcgGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACG‑TTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCt                                                                                                                                                                                                               >  1:453352/1‑251 (MQ=255)                                        aacaaacaACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACG‑TTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGcc                                                                                                                                                                                                          >  1:305530/1‑250 (MQ=255)                                          aaacaaacaGTAACCGGTCAATCCCTTACGCCAGCGGTTTTTTCCCGCGACT‑GTGTAGCATCAATATCAGGAAAGTAGCGATCCCCGTCAGCACCTCGTCCTGCAACAACGAAACTCCAAGAATGAGCTTAAAACAGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAACAGAACTAAACGACGG‑GACGCGTATAGTGATCTCGAATCTCCTCCGCCAGa                                                                                                                                                                                                       <  2:72222/246‑1 (MQ=255)                                                acaGTAACCCGCCAATCCCTTTCTCCAGGGGTTTTTGCCCGCGACG‑TTGCAGCATTAAACTCAGGAACGCCGCGCTCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCTGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑AACGCGGATAGTGATCGCGACTCTGCTCCGCCGGAtttttc                                                                                                                                                                                                 >  2:117065/1‑250 (MQ=255)                                                     aaCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACG‑TAGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACAACAAGAATGCGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTCAGG‑CACGCGGATCGTGATCGCGAATCTGATCCGCCAGATTTTTACCggg                                                                                                                                                                                            >  2:205706/1‑250 (MQ=255)                                                           ctaaacACTTTCTACTGCGGTCTGTTCCCGCGACGGTTGC‑GCACTAAAAACAGGAAAGTCGGGAGCCCCGTCAGCACCGCGACCTGCAACAACGAAACACCAAGAATGAGTATAAAACCAATCGCCGCACCAATAAAGTACGCCAGGTCGGTTGCCAGCGCAATAATTGCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTTCCGGTGGCa                                                                                                                                                                                        <  2:573510/244‑1 (MQ=255)                                                             aaTCAATTTCTCCAGCGGTTTTTGCCCGCGACG‑TTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GCCGCGGATAGTGATCGCGAATCTGCTCCGCCCGATTTTTACCGGTGGCAATcc                                                                                                                                                                                    >  2:630628/1‑251 (MQ=255)                                                                         cAGCGGTTTTTGCCCGCGACG‑TTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGc                                                                                                                                                                          >  2:244494/1‑249 (MQ=255)                                                                           gCGGTTTTTGCCCGCGACG‑TTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACCACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCag                                                                                                                                                                        >  1:230051/1‑249 (MQ=255)                                                                             ggTTTTTGCCCGCGACG‑TTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCagag                                                                                                                                                                      >  1:87463/1‑249 (MQ=255)                                                                                ttttGCCCGCGACG‑TTGCAGCATTAAAATCAGGAAAGTCGAGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCAAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGAt                                                                                                                                                                 >  1:502697/1‑251 (MQ=255)                                                                                    gCCCGCGACG‑TTGCAGCATTAAAATCAGGAAATCCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGa                                                                                                                                                             <  2:170575/251‑1 (MQ=255)                                                                                       cgcgACG‑TTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATc                                                                                                                                                          >  1:628662/1‑251 (MQ=255)                                                                                                     caTTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCACCCTGCAACAACGAAACACCACGAATGAGTTTAAAACCGATCGCCGCACCAATAACTTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACACAATAGAACTACACTACGG‑GCCGCGGATAGTGATCGCGACTCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGCTCCGAATCAGCATCGccctc                                                                                                                                              >  2:328204/1‑247 (MQ=255)                                                                                                            aTCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGTGAGGATCTGAATCAGCATCGCCATCAGGTTGGc                                                                                                                                      >  1:354831/1‑251 (MQ=255)                                                                                                                          gATCCCCGTCAACACCGCGCCCTGCAACAACGACACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATACCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGACTCTGCTCCGCCAGCTTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGCTCTGAATCAGCATCGCCATCAGGTTGCCCCACACAACAAccc                                                                                                                        >  2:43834/1‑251 (MQ=255)                                                                                                                                            gcCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAgg                                                                                                                                           >  1:359192/1‑214 (MQ=255)                                                                                                                                            gcCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGGCGGTTGCCATCGCAATAAGTTCTGCCTGAACCCAATAGATCTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAgg                                                                                                                                           <  2:359192/214‑1 (MQ=255)                                                                                                                                                      aacGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCAcc                                                                                            <  1:630628/251‑1 (MQ=255)                                                                                                                                                        cGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCACCCg                                                                                          >  1:400232/1‑251 (MQ=255)                                                                                                                                                         gAAACCCCAAGAATGAGTTTAAAACCTATCCCCGCACCAATAAATTCCGCCAGGGCGGTTGCCAGCGCAATAATTTCTGCCTGAACCCAATAGAACTACATTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCAGCAGGTGGTCCCAAACGCCAACCCGCAGTAGCTGATAGCCGAAGCTAGCACCCGc                                                                                         <  2:653302/251‑1 (MQ=255)                                                                                                                                                          aaaCACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCGGATAGCCGAAGCTAGCACCCGcc                                                                                        <  1:504373/251‑1 (MQ=255)                                                                                                                                                                 aaGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAatat                                                                                 >  1:88464/1‑251 (MQ=255)                                                                                                                                                                      tGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTc                                                                            >  1:320485/1‑251 (MQ=255)                                                                                                                                                                              aaaCCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGtt                                                                    <  1:302573/251‑1 (MQ=255)                                                                                                                                                                               aaCCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGACGCTAGCACCCGCCTGAATATTGGTCGCAAAGtt                                                                    >  1:628388/1‑250 (MQ=255)                                                                                                                                                                                   gATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCgg                                                               >  1:271387/1‑251 (MQ=255)                                                                                                                                                                                          gcACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCGat                                                        <  1:513674/251‑1 (MQ=255)                                                                                                                                                                                            aCCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCCCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGTATCGCCATCAGGTTGGCACAAACGACAACCCACAGTAGCTGATAGCGGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCGatat                                                      <  1:205706/251‑1 (MQ=255)                                                                                                                                                                                                  aaaTACCGCCAGGGCGGTTGCCACCGCAAAAATTTTTGCCTGAACCCAATAGAACCACACTTCGG‑GACGCGGAAAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTTAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCGATATAACCaa                                                <  2:230051/251‑1 (MQ=255)                                                                                                                                                                                                        cGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACTACACGACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCGATATAACCAATcgccg                                          <  1:33242/251‑1 (MQ=255)                                                                                                                                                                                                              gTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCACTAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCTGCTAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGAAAAGTTACCGGGATCGATATAACCAATCGCCGCAATg                                     >  2:332798/1‑250 (MQ=255)                                                                                                                                                                                                                                aatTTCTGCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCACCCGCCt                                                                                       <  1:619238/182‑1 (MQ=255)                                                                                                                                                                                                                                   ttcggCCTGAACCCAAAAGAACGACACTACGGTGA‑GCGGATAGGTATCGCGAAGCTGCTCCGCAATATTTTTATCGTTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACTACAACCCACAGTAGCTGATATCCGAAGCTAGCCCCCGCCTGAATATTGGTCGCAAAGGTACCGGGATCGATATAACCAATCGCCTCAATGAACGCAGGTACCATTAATGCg                <  2:340456/246‑1 (MQ=255)                                                                                                                                                                                                                                      tGCCTGAACCCAATCGAACTACACTACGG‑GACGCGGATAGTGATCCCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGcc            >  1:177971/1‑251 (MQ=255)                                                                                                                                                                                                                                       gCCTGAACCCAATAGAACTACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGcc            <  1:244494/250‑1 (MQ=255)                                                                                                                                                                                                                                                caaaaGAACTACACTACGG‑GGAGGGGATAGTGTACGCGATCCTGCCGCGCCAAATTTTTACCGGTGGCAATCCCGAGTTGGGCAGAGAGGATCTGAATCAGCATCGCCAGCAGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGTGCGATATATCCAATCGCCGCAATGAACGCAGGGCCCATTAATGCGAGCCTCATCTTgcg  <  2:177971/247‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   AGGCTGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACG‑TTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCATCGCAATAATTTCTGCCTGAACCCAATAGAACCACACTACGG‑GACGCGGATAGTGATCGCGAATCTGCTCCGCCAGATTTTTACCGGTGGCAATCCCTAGTTTGGCAGAGAGGATCTGAATCAGCATCGCCATCAGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATTGGTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCG  >  NC_000913/2512176‑2512664

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