Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,304,796	T→C	15.1%	intergenic (+405/+310)	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,796	0	T	C	15.1%	125.4 / 7.3	56	intergenic (+405/+310)	eco/mqo	ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base T (29/17);  new base C (4/6);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCAT  >  NC_000913/2304509‑2305082                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACCTCCAACGt                                                                                                                                                                                                                                                                                                 >  1:553817/1‑288 (MQ=255) cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGt                                                                                                                                                                                                                                                                                                 >  1:95107/1‑288 (MQ=40)    tgAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGccca                                                                                                                                                                                                                                                                                              >  1:682678/1‑287 (MQ=18)         tGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCCCTGTGAAGTGCTCCACCCCGTAGGTCGGATCAGGCGGTTACGCCCCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTACCAGCCCTCCCCCCCTGTGCACTGCTCCACCCCGCAGGCTGCATCAGGCGGTTACCCCCCCCCCGACCTCCTACGCCCTCGCCgg                                                                                                                                                                                                                                                                                        >  2:25498/1‑289 (MQ=14)           cGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCCACACCGCTGTGACGTGCCCCCCCCCGTCGGTCGGCTAAGCCGGTTACGCCGCATCCGACATCCCACCTCCGCGCCGCtt                                                                                                                                                                                                                                                                                      >  2:672190/1‑289 (MQ=11)           cGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGtt                                                                                                                                                                                                                                                                                      >  1:14771/1‑289 (MQ=40)           cGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCCCCACCCCGTAGGTTGGGTAGGCCGGTCACGCCGCATCCCACAACCAACCCCCGAGCCAGtt                                                                                                                                                                                                                                                                                      >  2:345629/1‑289 (MQ=11)                   ggTAGGTCAGAACAGACGGTCAGGCCGCATCCGACGCCCGACGCCCGAGCAGGTTGCCGGATGGGAGGCTGGCGCGTCTTAGCAGGCCTAAACCGCTGTGAAGGGCTCAACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATg                                                                                                                                                                                                                                                                              <  2:407024/288‑1 (MQ=11)                             aCAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTTgc                                                                                                                                                                                                                                                                    >  1:141519/1‑289 (MQ=37)                             aCAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgc                                                                                                                                                                                                                                                                    >  1:549258/1‑289 (MQ=37)                              cAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCGAGCCGGTTGcccgcatccgacatcca                                                                                                                                                                                                                                                                    <  1:313840/288‑16 (MQ=9)                                 ggCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGATGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGcgcg                                                                                                                                                                                                                                                                 >  1:652929/1‑288 (MQ=25)                                               gtccataagccgacggccagacgggTTGCCTGATGCGACGCTGGCGCGTCTTATGAGGCGTACACCGCGGTGAAGTGCGCCACCCCGAAGGTCGGGTAAGGCGGTTAGGCCGCATCCGACATCTAACGCCCAAGCCGGT‑GGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTa                                                                                                                                                                                                                                                  <  2:358348/265‑1 (MQ=11)                                                    aCATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCTCTGTGAAGTCCCCCACCCCGTAGGTCGGATACGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGACGCT‑TTCCTGATGCGACGCTGGCTCGTCTTCTCAGGCCACCCCCCCCGTGAAGTCCTCCACCCCGTAGGTTGGTTTAGGCGGTTCCGCCACATCCCAACTTCAACCTCCCAGCCGGTTGCCTGTTGGCACGCTGGCTCGTCTTacacggacaccacct                                                                                                                                                                                                                                              >  2:614696/1‑274 (MQ=11)                                                    aCATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCCGGCGCGTCTTATCAGGCCTACACCGCTGCGAAGCGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACAACTAACGTCCAAGCCGGTTGCCCGTTCCGACGGCGGCCGGTCTTACCAC‑GCCTACAccc                                                                                                                                                                                                                                             >  2:98420/1‑287 (MQ=11)                                                                   gCCGGTTCCCAGAGGGGAGGCCGGGGCGTCTTACCAGGCATACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCGAGCCGGT‑TGCCTGCTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCa                                                                                                                                                                                                                              <  2:95107/288‑1 (MQ=11)                                                                        ttGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGACGTGCCCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGTCTCCGACATCCCACGCAAGAGCCTGTCGCCTGATGCGCCGCTGTCGCGCCTTATCAGCCCCTAC‑CCGCTGTCACTTGCccccccccc                                                                                                                                                                                                                         >  2:186807/1‑287 (MQ=11)                                                                                  acaggtgtgcgagcctgagcAGGCCGACACCGCTGTGAAGTGCTACACCCGTTAGGTCGGATAAGGCTGTGACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAt                                                                                                                                                                                                               <  2:553817/269‑1 (MQ=11)                                                                                   gacccgtgcccggcTTCTCGGGCCTTCACCGCTGTGAAGTGGTCACCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑GGCCTGATGCGACGCTGGGGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa                                                                                                                                                                                                             <  2:307827/276‑1 (MQ=11)                                                                                     cGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTaca                                                                                                                                                                                                                                                <  2:193511/252‑1 (MQ=37)                                                                                     cGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTaca                                                                                                                                                                                                                                                >  1:193511/1‑252 (MQ=37)                                                                                     cGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTaca                                                                                                                                                                                                                                                >  1:531148/1‑252 (MQ=37)                                                                                     cGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTaca                                                                                                                                                                                                                                                <  2:531148/252‑1 (MQ=37)                                                                                                             cacCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCCTATCCG‑GCCTCCACCGCTCTGCAGTTCTCCACCCCGTAGGGCGGATACGTCGGTTACGCTGCATCCCACACCTa                                                                                                                                                                                   >  2:548394/1‑289 (MQ=11)                                                                                                                  cTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCg                                                                                                                                                                                                                                                                                              >  2:230935/1‑178 (MQ=35)                                                                                                                  cTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCg                                                                                                                                                                                                                                                                                              <  1:230935/178‑1 (MQ=35)                                                                                                                  cTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGc                                                                                                                                                                               >  1:407033/1‑288 (MQ=255)                                                                                                                                 ccccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                           >  1:225395/1‑198 (MQ=34)                                                                                                                                 ccccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTAGCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                           <  2:225395/198‑1 (MQ=21)                                                                                                                                              aTAAGGCGGTTACGCCGCATCCGACATCTAACGCCAGCGCCGGTGTGAC‑GAAGGGACGCTGGCGCGTCTGATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                                                   <  2:407033/288‑1 (MQ=11)                                                                                                                                                       tgaccgcgcagccgtagtccatcGCCCGACCCGGT‑GGCCGGATGCGCGACTGGCGCGGCTGATGAGGCCTAGACCGCTGTGAAGTTCTCCACCCCGTAGGTTTGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCGGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                         <  2:652929/267‑1 (MQ=11)                                                                                                                                                        gaCGCCGCCTCCGGCATCTATCGCCAGAGCCGGT‑TGCCTGATGCGACGCTGGGGCGTCTTATCAGGCCTCCACCGCTGTGAAGAGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATc                                                                                                                                        <  2:549258/288‑1 (MQ=255)                                                                                                                                                             cgcATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                                                            <  1:539480/264‑1 (MQ=255)                                                                                                                                                             cgcATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGGTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                                                            >  2:539480/1‑264 (MQ=255)                                                                                                                                                                         cTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                           <  1:104562/158‑1 (MQ=35)                                                                                                                                                                         cTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACCCCGCTGCGAACTGCTCCCCCCCGTAGGTCGGCTAAGGCTGTTACGCCGCATCCCACCTCTAACCCCCAAGCCGGTTGCCTGATGCGCCGCTGGCGCGTCTTATCCGGCCTACCCCGCttg                                                                                                                        >  2:438640/1‑286 (MQ=14)                                                                                                                                                                           aaCGCCCGAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCCTATCAG‑CCCTACACCGCTGTGAAGTGCTCAACCCCGTAGGTCGGACAAGGCGGTCACGCCCCATCCGACCTCTAACGCCCAAGCCGGTTGCCTGATCCGCCGTCTGCGCGTCCTATCCGGCCCACCCCGCttcgt                                                                                                                      >  2:302424/1‑284 (MQ=11)                                                                                                                                                                                ccgaGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTaca                                                                                                                                                                                                                                                >  1:363913/4‑161 (MQ=14)                                                                                                                                                                                   aGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTaca                                                                                                                                                                                                                                                >  2:104562/1‑158 (MQ=38)                                                                                                                                                                                         t‑tGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                           <  2:260318/142‑1 (MQ=37)                                                                                                                                                                                         t‑tGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                           >  1:260318/1‑142 (MQ=37)                                                                                                                                                                                                                               cacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                           <  2:341369/105‑1 (MQ=40)                                                                                                                                                                                                                               cacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                           >  1:341369/1‑105 (MQ=40)                                                                                                                                                                                                                                    cTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAATGCGGTTACGCCGCATCCGACATCTACCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCTCTACCCCGCTGTGACTTGCTCAACCCCGGCGGTTCGATAACGCGTTTCCCCCCCATGCAACATCCAACGc                                                              >  2:122743/1‑288 (MQ=11)                                                                                                                                                                                                                                    cTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCACGCCGGTTGCCTGATGCGACGCTTCCGCGTCTTATCCGCCCTACCCCGCTCTTCAGTTCCCCACCCCGCAGGTTGGATAAGTCCGTTACGCCGCCTCCAACTTCTAACGc                                                              >  2:284770/1‑288 (MQ=14)                                                                                                                                                                                                                                    cTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                               <  2:160898/223‑1 (MQ=40)                                                                                                                                                                                                                                    cTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                               >  1:160898/1‑223 (MQ=40)                                                                                                                                                                                                                                                      cGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTCCACCGCTCTGCAGTGCCCCACCCCGTAGGTTGGATTATGCGGTTACGCCGCCTCCAACCTCTAACCCCCGCCCCGGTTGCCTGAt                                            >  2:106989/1‑288 (MQ=14)                                                                                                                                                                                                                                                        tAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAAAGCCCGAGCCGGTTGCCTGATGCg                                         >  1:446023/1‑289 (MQ=39)                                                                                                                                                                                                                                                              ggATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTg                                   >  1:307827/1‑289 (MQ=37)                                                                                                                                                                                                                                                                aTAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                         <  1:333807/185‑1 (MQ=34)                                                                                                                                                                                                                                                                aTAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTaca                                                                                                                                                                                                                                                <  1:201519/82‑1 (MQ=17)                                                                                                                                                                                                                                                                  aggggggggaCGGCGGATCCGACAGACACGGCCCGAGGCGTTTGAGTGCGGGGAGGCTGGCGGGTGTTGTCAG‑GGCGAAAAGGCTGTGAAGGGCACCACCCCTTAGGTCGGATAAGGCGGTTACGCCTCATCCGACATCCAACGCCCAAGCCGGTTGCCTGATGGGCCGCTGGGGCGTCTTATCAGGCCTAAACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGgcgc                               <  2:14771/280‑1 (MQ=11)                                                                                                                                                                                                                                                                     gCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCt                                                                                                                                                                                                                                 <  2:350029/92‑1 (MQ=37)                                                                                                                                                                                                                                                                          tACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATc                       >  1:274103/1‑289 (MQ=34)                                                                                                                                                                                                                                                                               cgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACCTCTAACGCCCGAGCCGGTTTCCTTATGCCACGCTCGCCCGTCTTCTCAGCAt                  >  2:25197/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                   tCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACg               >  1:438242/1‑288 (MQ=37)                                                                                                                                                                                                                                                                                    ccGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGcc             <  1:345629/289‑1 (MQ=38)                                                                                                                                                                                                                                                                                               cGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGAAGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACCCCCGAGCCGGTTGCCTGATGCGACGATCGCGCGCCTTATCAGCCTTCGCCCCCTCCGGCAt  >  1:410919/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                CTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGT‑TGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCAT  >  NC_000913/2304509‑2305082

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