Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,282,627	C→T	Q922* (CAG→TAG)	narG →	nitrate reductase A subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,282,627	0	C	T	85.7%	112.8 / ‑0.4	50	Q922* (CAG→TAG)	narG	nitrate reductase A subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (6/1);  new base T (12/30);  total (18/32)
Fisher's exact test for biased strand distribution p-value = 7.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCC‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCA  >  NC_000913/1282386‑1282869                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        tCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTTGCTGGGACAACC‑AGAGCgggg                                                                                                                                                                                                                                               >  1:582635/1‑248 (MQ=255)   tGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGGTTGCCTTGAACACAC‑AGAGCGAGATg                                                                                                                                                                                                                                             >  1:440268/1‑251 (MQ=255)       gcaaGAAACCGACATCGTCACGCTGCCTATCCAGCATGTCTCTGCCGCTGAACGGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGACCACCT‑AGAGCGAGATGGATc                                                                                                                                                                                                                                         <  2:119506/249‑1 (MQ=255)                cGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAAGGATTGACTGGAAAAACT‑AGAGCGAGAtcggattgttgcggg                                                                                                                                                                                                                                >  1:248158/1‑237 (MQ=255)                 gACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGGTTGCCTGGAACACCC‑AAAGCGAGGAGGAACTGCTG‑Cgga                                                                                                                                                                                                                               >  1:319263/1‑248 (MQ=255)                 gACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGGTTGCCTGGAACAACT‑AGAGCGGGATGGGTCTGCTG‑CGTa                                                                                                                                                                                                                               >  1:192138/1‑250 (MQ=255)                         cACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATGATCCGGCGACTTACGAACGCTTTACCTCTATCTGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACt                                                                                                                                                                                                                      <  1:579755/251‑1 (MQ=255)                                       aGCAGGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTAGCGGCCCGCTGATGGAGAAAAGCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGATCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGaa                                                                                                                                                                                                        <  2:68447/251‑1 (MQ=255)                                                  gCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGc                                                                                                                                                                                                <  2:191157/248‑1 (MQ=255)                                                  gCCGCTGAACTGGCGCAGCCGCTGGATGAGAAAGACTGGAAAAAAGGCGAGTGCGGCCTGATCCCAGGTAAAACCGCGCCACACATTAGGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCg                                                                                                                                                                                               <  2:314686/249‑1 (MQ=255)                                                       tGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAAGGAGTGCCTGGGACAACC‑AGAGCGCGAGTGAACTGCTG‑CGTAAGCCCAACTACACCCAAACGGAAGGGCCGGCgcaaagcc                                                                                                                                                                                        >  1:60468/1‑244 (MQ=255)                                                             ggCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGcc                                                                                                                                                                                    <  1:272667/249‑1 (MQ=255)                                                                 cAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCt                                                                                                                                                                             <  2:112873/251‑1 (MQ=255)                                                                  aGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTg                                                                                                                                                                            <  2:499638/251‑1 (MQ=255)                                                                     cGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AACGCGAGGTGGGTCTGCTG‑CGTAAGCCCCACTACACCAAAACGGGAGGACCGGAGGAAGGCCCGCCG‑ctgctgcac                                                                                                                                                                         >  2:156050/1‑248 (MQ=255)                                                                                 aaGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAAGGATTGGCCGGAAAAACT‑AGATCGAGATGGCTATGGTT‑GGGAAGCGCAACTACCCCAAAACGGAAGGAGCGGCGGAAGGGCAGCAG‑ATACTGCAACAC‑GCAAtata                                                                                                                                                             >  1:541808/1‑247 (MQ=255)                                                                                 aaGACTGGAAAAAAGGCGAGTGCGACCTGACCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGAt                                                                                                                                                            <  1:229728/251‑1 (MQ=255)                                                                                    aCTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGGAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCg                                                                                                                                                         <  2:403790/251‑1 (MQ=255)                                                                                              aGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGAGGTTGCCTGGAAAACCT‑AGAGCGAGATGGATCTGCTG‑CGGAAGGTCCACTACACCCAAGCGGAAGGTGCGGCGCGAGGGCAGCCGGAGGCT‑‑AAAACCCCAATTAATGCGGCCGAgatg                                                                                                                                                >  2:399175/1‑250 (MQ=255)                                                                                                 cgaatGCGACCTGATCCCAGTTAAAACCGCGCCACACATTAGGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTGTTCCTCTATCGGCACGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACCCCAAAGCGGAAGGTCCGGCGAAAGGCCAGACG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATcc                                                                                                                                            <  1:9944/247‑1 (MQ=255)                                                                                                 cGAGTGCGACCTGATCCCAGTTAAAACCGCGCCACACATGTGGGTCGTAGAGCGCGATTATCCGGCGACTTACGAAGGCTTTACCTCTATCGTCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGTGATTGCCTGGAACACCT‑AGAGCGAGATGTATCGGCTG‑CGTAAGCTCATCTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATcc                                                                                                                                            <  2:255532/251‑1 (MQ=255)                                                                                                   aGTGCGACCTGATCCCAGGTAAACCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAAGGATTGCCCGGCACACCC‑AGGGGGGGAAGGGGCTGCTGGCG‑AAACTCCACGACCCCAAAGAGGAAGGGCAGGCGAAAAGCCAGCCG‑GAGCTG‑GACACCGCGATTGATGGGGCAGGgatggaaccg                                                                                                                                          >  1:154133/1‑245 (MQ=255)                                                                                                       cGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAAGGCTTGCCCCGAACACCT‑ATAGCGAGATGGATCTGCTG‑CGTTAACTCCACTCCACCCAAGCGGGAGGTGCGGCGAAAGGCCCGCCG‑CTGCTG‑AAAACCCCCATTGATGCGGCcgagctaatcccgcaac                                                                                                                                      >  2:286655/1‑238 (MQ=255)                                                                                                                   agggaaaaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCTACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCTAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGTCGAAAGGCCAGCCG‑ATGCGG‑AAGACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGaaa                                                                                                                         <  1:236611/247‑1 (MQ=255)                                                                                                                    ggTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAAc                                                                                                                        <  2:319263/251‑1 (MQ=255)                                                                                                                        aaaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGACAATCGGTAATGGCGGCAAAGGGATTGCCTGGAACACCT‑AGAGCGCGATGGATCTGCTG‑CGTAAGCTCCACTAAACCAAAACGGAAGGTCCGGCGAAAGGCCAGACG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACAATGGCACCGGGAaacaac                                                                                                                    >  1:559897/1‑251 (MQ=255)                                                                                                                        aaaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCGCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAAc                                                                                                                    <  1:433822/251‑1 (MQ=255)                                                                                                                           ccGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGt                                                                                                                 <  1:423807/251‑1 (MQ=255)                                                                                                                                   cacaTTATGGGCGTAGGGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCG‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACGCCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGTCAGGTAGc                                                                                                         <  2:482034/251‑1 (MQ=255)                                                                                                                                          tGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAAGGACTGGCCGGGACACCT‑AGCGAGAGGTGGATCTGCCG‑CGTAAGGTCCACTCCACCAAAACGGACGGGCCCGCGGAAGGCCAGACG‑ATGCTG‑AACACCCCAATTGATGGGGCAGAGAAGATCCCGCACAC‑GGCACCGGGAACCACCGGTCAGGTAGCCgcgga                                                                                                   >  1:516853/1‑246 (MQ=255)                                                                                                                                                 agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTggg                                                                                           <  2:584183/251‑1 (MQ=255)                                                                                                                                                         aTTATCCGGCGACTTACGAACGCTTTACCTCTAGCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATGGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGGCCGGCGAAAGGCCAGCCG‑ATGCTG‑AGCACCGCAGTTGATTCGGCAGATATGATCCTG‑TCAGTGGCTCCGGAAGCCATCGGTCAGGTACCCGTTAAAGCCTGGGCTTCCCTg                                                                                   <  2:87751/251‑1 (MQ=255)                                                                                                                                                                         cGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTc                                                                   <  1:456014/251‑1 (MQ=255)                                                                                                                                                                                 ttACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGGTTGCCTGGAACACCC‑AGAGCGAGAAGGATCTGCTG‑CGTAAGATCCACTACACCAAAACGGAAGGGCCGGCGACAGGACCGCCC‑ATGCTG‑AAAACCGCAAATGATGCGGGAGAGATGATCCTG‑ACACTGGCACCGGAAAACAACGGGCAGGGAGCCGGGAAAACCCGGGCTGCCCTGAGCGGATTTACCGGGCGTGACCac                                                           >  1:240135/1‑250 (MQ=255)                                                                                                                                                                                       ctATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGTCAGGGAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTGACCGGTCGTGACCATACGc                                                       <  2:240135/249‑1 (MQ=255)                                                                                                                                                                                       ctATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCGG‑CGTAAGCTCAACTACACCAAAGCGTAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGACTGCCCTGAGCGAATTTACCGGTCGTGACCATACGc                                                       <  2:124198/249‑1 (MQ=255)                                                                                                                                                                                           cGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTgg                                                 <  2:192138/251‑1 (MQ=255)                                                                                                                                                                                            ggCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAgaga                                                                                                                                                  <  2:130461/154‑1 (MQ=255)                                                                                                                                                                                            ggCCCGCTGATGGAGAAAATCGGTAATGGCGGAAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGGAAGCTCAACTACACCAAAACGGAAGGTCCGGCGAAAAGCCAGCCG‑ATGCTG‑GACACCGCAATTGATGCGGCAgaga                                                                                                                                                  >  1:130461/1‑154 (MQ=255)                                                                                                                                                                                             gCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGgcg                                               <  2:408092/251‑1 (MQ=255)                                                                                                                                                                                              cccgcgtaTGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACCCCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCT‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGgcgc                                              <  2:502880/244‑1 (MQ=255)                                                                                                                                                                                                       ggAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACGCCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACGCTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAaag                                     <  1:240505/251‑1 (MQ=255)                                                                                                                                                                                                        aagaAAATCGTGAATGGGCGTAAAGGGATTGCC‑CGAAAACCCTAGAACGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCTGCAGAGATGAGCCTG‑ACACTGGCAACGGAAACCAACGTTCAGGTAGCTGATAAATCCTGGGCTGCCCGGAGCAAATTTGCCGGTCGTGCCCATACGCATCTGGCCCTTaagaaaga                                    <  2:292572/250‑1 (MQ=255)                                                                                                                                                                                                              aTCGGTAATGGCGGTAAAGGGAGTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGGAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACga                              <  2:153727/251‑1 (MQ=255)                                                                                                                                                                                                                            tAAAGGGATTGCCTGGCACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGACGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTcc                >  2:351035/1‑251 (MQ=255)                                                                                                                                                                                                                                   aTTGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCGTCCGAGatat         <  2:362373/251‑1 (MQ=255)                                                                                                                                                                                                                                     tGCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACCCCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATAtt        <  1:119076/250‑1 (MQ=255)                                                                                                                                                                                                                                      gCCTGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGATAAGATCCGCTTCCGCGATATTCa      <  1:169920/251‑1 (MQ=255)                                                                                                                                                                                                                                         tGGAACACCT‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCCAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGc   >  2:406022/1‑251 (MQ=255)                                                                                                                                                                                                                                          ggAACACTT‑AGAGCGAGATGGATCTGCTG‑CGTAAACTCAACTACACCAAAGCGGCAGGTCCGGCGACAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGCTGATCCTG‑AAACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCAGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGcc  >  2:235794/1‑250 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        TCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCC‑AGAGCGAGATGGATCTGCTG‑CGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCG‑ATGCTG‑AACACCGCAATTGATGCGGCAGAGATGATCCTG‑ACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCA  >  NC_000913/1282386‑1282869

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