Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,304,549	C→T	intergenic (+156/+559)	eco → / ← mqo	serine protease inhibitor ecotin/malate:quinone oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,304,549	0	C	T	92.9%	20.4 / ‑3.2	14	intergenic (+156/+559)	eco/mqo	serine protease inhibitor ecotin/malate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (5/8);  total (5/9)
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.60e-01

TATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGT‑CTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAC  >  NC_000913/2304262‑2304847                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                            tATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCTGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGACCTACACCGCTGCGACGTGCTCCACCCCGTAGGTTGGATACGGCGGTTACGCCGCATcc                                                                                                                                                                                                                                                                                                  >  1:843664/1‑299 (MQ=25) tATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATcc                                                                                                                                                                                                                                                                                                  >  1:843524/1‑299 (MQ=255)                        gcGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGtt                                                                                                                                                                                                                                                                           >  1:31633/1‑299 (MQ=255)                                             gATCGTGGTGTATACGCCAGACAAGGTAGATGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                                                                                                                                                     <  2:712976/300‑1 (MQ=17)                                                     tgtATACGCCAGACAATGTAGCTGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCTTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGTCCTACACCGCTTTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGT‑CTTATCAg                                                                                                                                                                                                                                            <  2:407899/300‑1 (MQ=17)                                                      gtATACGCCAGACAATGTAGATGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGT‑CTTATCa                                                                                                                                                                                                                                             >  1:108577/1‑298 (MQ=255)                                                           cGCCAGACAATGTAGATGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGT‑CTTATCAGGCCt                                                                                                                                                                                                                                        >  1:270754/1‑298 (MQ=255)                                                               aGACGATGTAGATGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACACCGCGGTAGTTCGTTAAACTGCCGTGAAGTGTGGCAGCTAGTGGGTTAGATAAGGCGGTCATGTGGCATCTGATATCCAACGCCCTAGCTGGTTGCCTGATGCTACGCTGGCGTGTCTTATCGGGCCTACACCGCTGTGAAGTGCGTCACCCCGTAGGTCGGGTAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGTGT‑TTTATCAggtttacatt                                                                                                                                                                                                                                   <  2:237580/299‑9 (MQ=12)                                                                  aatggTAGATGTGAAGNTACCGCGTCTGGAAGGCGNGAAGAGAGAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTNCTTATCAGGCCTacac                                                                                                                                                                                                                                    <  1:700276/295‑1 (MQ=14)                                                                   acgctAGATGTGAAGNTACCGCGTCTGGAAGGCGNGAAGAGAGAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTNCTTATCAGGCCTacac                                                                                                                                                                                                                                    <  1:700277/294‑1 (MQ=14)                                                                           atgttaaggaccgcggtgggAAGGGG‑GAAGGTA‑AAATTGACAGCTCGTTAGTTCGCTAAACTGCCGTTAAATGCGGCACCCCTTAGGTCAGACACGTCGGTCACGCCGCATCCGACATCCAACGCCCGAGCTGGTTGCCGGATTCGACACTGTCGCTTCTGATCATGCCTACACCGCTTTTAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCAAGCCGGTTGCCTGAGGCGACGCTGGCGCGT‑CTTATCAGGCCTACACCGCTGTGAAGTGc                                                                                                                                                                                                                       <  2:463229/282‑1 (MQ=2)                                                                                         gTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGGAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGT‑CTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGt                                                                                                                                                                                                         <  1:415940/300‑1 (MQ=11)                                                                                             ggAAGGCG‑GAAGAGA‑AAATCGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGTTCAGACAAGGCTGTCACGCCGCATCCGACAGCCAACGCCCGAACCGGTTGCCTGATGTGACACTGGCGCGTCGGATCATGCCTACACTGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGT‑CTTATCAg                                                                                                                                                                                                                                            <  2:409433/261‑1 (MQ=11)                                                                                                                                                                                                                                                                                                  tCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGT‑CTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTacac  <  1:204731/300‑1 (MQ=39)                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                            TATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGT‑CTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAC  >  NC_000913/2304262‑2304847

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