Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,304,535	A→G	17.6%	intergenic (+144/+571)	eco → / ← mqo	ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,304,535	0	A	G	17.6%	62.3 / 4.1	40	intergenic (+144/+571)	eco/mqo	ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base A (19/12);  new base G (4/4);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTC  >  NC_000913/2304310‑2304747                                                                                                                                                                                                                                    |                                                                                                                                                                                                                     ggttgataCGCCAGACAATGTAGATGGGAAGTACCGCGTTTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTGCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTC‑GGATAAGGCGGtta                                                                                                                                                                                                         <  1:314410/234‑3 (MQ=18)       taCGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCCACGCCCGCGCCGGTTGCCCGATGCGACGCTGGCGCGTCTTATCAGGCCTCCACCGCCGGGAAGTGCGGCACCCCGTAGGCC‑ataaccgccggtaacccccc                                                                                                                                                                                                   >  2:842298/1‑220 (MQ=21)       taCGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCCCCCCGTAGGTC‑AGCAAAGGCGGGCAccccc                                                                                                                                                                                                    >  1:1325652/1‑237 (MQ=255)       taCGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCAcgccg                                                                                                                                                                                                    >  1:494364/1‑238 (MQ=255)                   tGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGCAATCCa                                                                                                                                                                                       >  1:852788/1‑239 (MQ=255)                    gTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGGCGCTGGCGCTGCTTTTCCGGCCTCCACCACGGGGAGGTGCGCCACCCCGCAGGGC‑GGACAAGggggacccccccccaccaccaacca                                                                                                                                                                                       >  2:1172977/1‑217 (MQ=18)                           tgAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACCTCCAACGCCCg                                                                                                                                                                                >  1:993372/1‑238 (MQ=255)                             aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑GGATAAGGCGGTTACCCCCCACCCgccccttaacgcccaca                                                                                                                                                                              >  1:1045504/1‑222 (MQ=17)                             aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGCCGGTCACGCCGCATCCGACATCCAACCCCCGAg                                                                                                                                                                              >  1:554121/1‑238 (MQ=255)                             aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGCCTTATCCGGCCTACACCGCTGTGAAGGGCGGACCCCCGCAGGGC‑GGACAAGGCGGTCACGCCCCACCCGACCCCAAACCCCCGacc                                                                                                                                                                             >  2:1361748/1‑237 (MQ=18)                                     cgTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACCCCGCATCCGCCATCCAACGCCCGACCCGGTTg                                                                                                                                                                       >  1:1152441/1‑237 (MQ=255)                                     cgTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGCC‑AGACAAGGCGGTCACGCCGCATCCGCCATCCAACGcccagaccgattg                                                                                                                                                                       >  1:1167289/1‑227 (MQ=255)                                                                cagcggggtcttcgcaaaaaggcccGTGAAGGGCGCCCCCCGGTGGGGCAGAGAAGGCGGTCCCCCCGCAG‑CCGCCCTCCACCGCCCGGCGGGGTGGGCTGAGGGGACGCGGCCGCGTCTTAGCAGGCCTACACCGCTGGGAAGTGCGGCCCCCCGTAGGTC‑AGAAAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAg                                                                                                                                          <  2:1152441/216‑1 (MQ=9)                                                                  aCGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAACCCGTTGCCTGTTGCGACGCTGGCGCGTCCTATCCgg                                                                                                                                         >  1:157834/1‑238 (MQ=255)                                                                             cGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGCCATCCAACGCCCGAGCCGGTTGCCCGCTGCGACGCGGGGCCGCCTTTTCAGGCCTACCCCCCtg                                                                                                                             >  1:809569/1‑239 (MQ=14)                                                                                  aaCTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCCTATCAAGCCTACACCCCTGGGAAg                                                                                                                        >  1:462131/1‑239 (MQ=255)                                                                                  aaCTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTCCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCCCTGTGAAg                                                                                                                        >  1:932975/1‑239 (MQ=255)                                                                                                 gccgccctttggggtcggggaaGGCGGTC‑CGCCCGATGCCGGCATCGAGCGCCCGGGCCGGTGGGCGGATGGGCAGCTGGCGCGGCGTATCCTGCCTACACCGGCGTGAAGTGCGGCACCCCGTAGGTC‑AGGCAAGGCGGCCACGCCGCAGCCGACATCCATCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAg                                                                                                         <  2:360376/219‑1 (MQ=9)                                                                                                     cccccGTAGGCCAGACGAGGCGGCCACCCCGCAG‑CCTACATCCAACGCCCGAGCCGGTGGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTc                                                                                                      <  1:108331/237‑1 (MQ=255)                                                                                                       cccgggtgtcgggcAAGGCGGTCACGCCGCAC‑CCGATATCCAACCCCCGGGTCGGTTGCGTGATGCGACGCGGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                                                                                                   <  1:624018/226‑1 (MQ=14)                                                                                                                   aCAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGGAGTGCTCCACCCCGCAGGTCGGATAAGGCGGttg                                                                                        >  1:454780/1‑237 (MQ=255)                                                                                                                                                     ccggggCCGGTTGCCTGATGCGCCGCTGGCGCGTCTTATCAGGCCTACCCCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCgg                                                         <  1:285334/231‑1 (MQ=255)                                                                                                                                                             ggTTGCCTGATGCGACTCGGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                               <  1:1037113/238‑1 (MQ=255)                                                                                                                                                                     tgtggtcaggcggcggcttttATCATGCCCACCCCGCTGGGAAGTGCGCCCCCCCGGAGGTC‑GGAGAGGGCGGTACGGCCGCACCCGACAGCAAACGCCCGAGCCGGTGGGCTGATGCGACGCGGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTgg                                      <  2:448045/220‑1 (MQ=11)                                                                                                                                                                               tttgcgtgttttagtggCCTACACCCGTGTGAAGGGCGCCCCCCCG‑GGGTCGGGACAAGGCGGTGACGCCGCACCCGACATCCAACGCCCGGGCCGGTTGCCGGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                            <  2:570936/224‑1 (MQ=11)                                                                                                                                                                               tGGCGGGTCTTATCAGGCCTCAACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTGCGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                            <  1:1112678/239‑1 (MQ=255)                                                                                                                                                                                  gggatcttcacaggccgctacCGTTGGGAAGTGCGCCACACCTGAGGTC‑CGACAGGGCGGACACCCCGCACCCGCCACCCACCGCCGGAGCCGGTTGCCGGATGCGACGCGGGCGCGTCTGATCAGGCCTACACCGCTGGGAAGTGCTCCCCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAg                         <  2:454780/220‑1 (MQ=11)                                                                                                                                                                                         tagcAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                  <  1:244733/236‑1 (MQ=255)                                                                                                                                                                                               gCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCCCCTCCGACACTCAACGCCCCAACCGGTTGCCCGATGCTGCGCTCGCGCGTCttttactgcctccacccccg            >  1:448045/1‑221 (MQ=14)                                                                                                                                                                                               gCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACCCCGCAACCGGCATCTAACCCCCAACCCGGTGTCCTGATTGGACGCTGGCGGGCCCTATCAAGCCTACAccc               >  1:1348268/1‑235 (MQ=14)                                                                                                                                                                                                   accccGCTGGGAAGTGCGGCCCCCCGGAGGGC‑AGACAGGGCGGTCACGCCGCAGCCGACATCCAACGCCCGAGCCGGTTGCCTGATGGGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGaa        <  1:842298/236‑1 (MQ=14)                                                                                                                                                                                                        tcctggaagggcgcacccccGTAGGTC‑GGCCAGGGCGGCCCCCGCGCATCCGCCAGCCAACCGCCGGGCCTGTGGCCTGATGCAACGCTGGCGCGTCTTAACAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCt   <  2:157834/223‑1 (MQ=11)                                                                                                                                                                                                         cTGTGAAGTGCGGCACCCCGTAGGTC‑GGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAACCCGTTGCCTGATGCGACGGTGGCGCGTCTTATCCAGCCTACACCCCTGTGAAGTggc   >  1:1161178/1‑236 (MQ=255)                                                                                                                                                                                                         cTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                             >  2:1232601/1‑100 (MQ=32)                                                                                                                                                                                                         cTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                             <  1:1232601/100‑1 (MQ=32)                                                                                                                                                                                                         cTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAACCCGTTGCCTGATGCGCGCCTGGCCCGGCTTTTAAAGCCCACAcccccggtgacagccc  >  1:591420/1‑225 (MQ=14)                                                                                                                                                                                                         cTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGAGTCTTATCa                                                                                                                                           >  1:777447/1‑102 (MQ=14)                                                                                                                                                                                                          tgtgAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                             >  2:521725/1‑99 (MQ=31)                                                                                                                                                                                                          tgtgAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                             <  1:521725/99‑1 (MQ=31)                                                                                                                                                                                                                   cGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCg                                                                             <  1:1219068/154‑1 (MQ=40)                                                                                                                                                                                                                   cGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCg                                                                             >  2:1219068/1‑154 (MQ=40)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                     GGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAT‑CCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTC‑AGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTC  >  NC_000913/2304310‑2304747

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