Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,304,522 2 bp→TC intergenic (+129/+585) 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_0009132,304,5220GT86.7% 19.5 / 0.9 15intergenic (+129/+586)eco/mqoserine protease inhibitor ecotin/malate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (10/3);  total (11/4)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01
*NC_0009132,304,5230GC86.7% 23.0 / 1.5 15intergenic (+130/+585)eco/mqoserine protease inhibitor ecotin/malate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base C (10/3);  total (11/4)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

CTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT  >  NC_000913/2304283‑2304757
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                          
ctgcGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGt                                                                                                                                                                                                                                    >  1:132921/1‑249 (MQ=255)
               aagttgcagcgcgttgtgtATACGCCAGAAAAGGTAGATGGGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACGGCCGTGAAGTGCGGCACCCCGTAGGTCAGAAAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGt                                                                                                                                                                                                                                <  2:46837/223‑1 (MQ=18)
               aaGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGt                                                                                                                                                                                                                                >  1:46837/1‑238 (MQ=255)
                  cTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAACCCGTAGGTCGGATAAGGCGGttac                                                                                                                                                                                                                >  1:351462/1‑248 (MQ=255)
                           gtggtgTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCCGGCGCGTCTTCTCCGGCCTACACCGCTGTGAAGTGCTCCACCCCGTATGTcggctaaagcgggtaatccgcactc                                                                                                                                                                                                       >  2:388491/1‑227 (MQ=21)
                            tggtgTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAGATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACCCCGCTGTGAAGTTCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATccc                                                                                                                                                                                                      >  1:257663/1‑250 (MQ=18)
                                     cGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGGCGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGCTACGCCGCATCCGCCATCTAAc                                                                                                                                                                                             >  1:156732/1‑251 (MQ=17)
                                             aaTGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTCCACCGCCGGGAAGTGCTCCTCCCCGGAGGTCGGATAAAGCGGTTAcgccccctacgccaactaacacccccgc                                                                                                                                                                                     >  2:114653/1‑227 (MQ=17)
                                                gTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACCTCCAACGcccaaaccgg                                                                                                                                                                                  >  1:111549/1‑244 (MQ=17)
                                                          aGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCCCCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                                                                        <  1:256188/251‑1 (MQ=14)
                                                                        aaGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCCCATCCGACATCTAACGCCCCAACCGGTAGACTGAAGCGGCGCTGGcgcgg                                                                                                                                                          >  1:550739/1‑250 (MQ=14)
                                                                                    aaaaTTGACAACGCGGTGGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCt                                                                                                                                              <  1:344042/251‑1 (MQ=11)
                                                                                     aaaTTGACAACGCGGTAGTTCGCTAAACTGCCGTGAGGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGTCGGTTGCCTGATGCGACGCTAGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAAGCCCCATCCGACATCTAACGCCCCAGCCCGTTGCCTGATGCGACGCTGGCGCGTCTTATccagcct                                                                                                                                              >  1:46186/1‑244 (MQ=11)
                                                                                                                                                                                                                      agttcgctaaactgccgtgAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTGGTGCGGCGCTGGCGCGTCTTATCAGGCCCACAACGCTGTGAAGTGc            >  1:322187/17‑251 (MQ=11)
                                                                                                                                                                                                                                actgccgtgAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGt  <  1:82327/245‑1 (MQ=255)
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                          
CTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT  >  NC_000913/2304283‑2304757

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