Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,649 C→T 15.5% intergenic (+256/+459) 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,6490CT15.5% 43.8 / 4.1 28intergenic (+256/+459)eco/mqoserine protease inhibitor ecotin/malate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (11/10);  new base T (2/2);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

TACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGT  >  NC_000913/2304504‑2304798
                                                                                                                                                 |                                                                                                                                                     
tACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCTACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAt                                                                                                                                                   <  2:31748/150‑1 (MQ=255)
 acacCGCCGTGAAGTGCGGCACCCCGTAGGTCAGACCAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGTTGCGACGCTTTCTCGTCTTATGAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATa                                                                                                                                                  <  2:1396694/150‑1 (MQ=14)
       cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg                                                                                                                                            >  1:1227406/1‑150 (MQ=40)
        tgtgAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc                                                                                                                                              <  2:608985/147‑1 (MQ=40)
        tgtgAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc                                                                                                                                              >  1:608985/1‑147 (MQ=40)
         gtgAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATTCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGtt                                                                                                                                          <  1:98421/150‑1 (MQ=255)
                   gCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATc                                                                                                                                >  1:850915/1‑150 (MQ=255)
                     aCCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTGCACCCCGTAGGTGGGAAGCGGCGGTAACGCCGGATCCg                                                                                                                              >  2:31010/1‑150 (MQ=11)
                       cccGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGAc                                                                                                                            >  1:263111/1‑150 (MQ=255)
                         cGTAGGGCCCACAAGGCGGTCACGCCGCATCGTACATCCATCGCTCGAGGAGGTTGCCTGTTTCTAATTTGGCGCGGCGTTTGTTGGCTACACTGCTGTGAAGTGCTGCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACAt                                                                                                                          <  2:1115017/150‑1 (MQ=2)
                            cgcccAGACAAGGCGGTCACGCCGCATCCGACATCCACCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCAGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTa                                                                                                                       <  1:512015/146‑1 (MQ=255)
                               tCAGCCAACTCGGTCACTCCTCATCCTACATCCAACGCCCGAGCCGGTTGCCTGATTCTACTCTGCCTCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACg                                                                                                                    <  2:751460/150‑1 (MQ=11)
                                cAGACATGGCTGTCACGCCGCATCCGACATCCAACGCCCGAGCCGTTTGCGTTATTGTAATCTGGCTCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGc                                                                                                                   <  2:1020409/150‑1 (MQ=11)
                                  gACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATGCGACATCTAACGccc                                                                                                                 >  2:1431231/1‑150 (MQ=37)
                                  gACAAGGCGGTCACGCCGCATCCGACATCCAACGACCGAGCCGGTTGACATCTGAAAAGCTTCCGCTTCTTAAAAAGACTACACAGCTGTGAAGATCTCAACACCGAAGATCGGATAAGAAAGTTACGCCGCATCCGACATCGAACGgcc                                                                                                                 >  2:81687/1‑150 (MQ=2)
                                      aGGCGGTCACGCCGCATCCGACATCCAACTCCCAAGCCCGTAGACTGATGCAACGCTGTCGAGTCTTATAAGGCAAACACCGCTGAGAAGTGATAAACAAAGTAGGTAGGATAAGGCGGTTACGCCGCATCCGACATCGAGCGCGCAAGc                                                                                                             >  2:668478/1‑150 (MQ=2)
                                           gTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCgg                                                                                                                                            <  1:412748/114‑1 (MQ=25)
                                           gTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCgg                                                                                                                                            >  2:412748/1‑114 (MQ=25)
                                             cACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGTCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                      >  2:418355/1‑150 (MQ=17)
                                                gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTg                                                                                                   >  1:608564/1‑150 (MQ=32)
                                                   gcATCCGACATCCACCGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAta                                                                                                <  1:1279957/150‑2 (MQ=17)
                                                               cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTAcgc                                                                                                                                      <  1:984706/100‑1 (MQ=14)
                                                               cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTAcgc                                                                                                                                      >  2:984706/1‑100 (MQ=14)
                                                                   gCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCt                                                                                <  2:263111/150‑1 (MQ=31)
                                                                                                                                      aCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCg             >  1:204090/1‑150 (MQ=11)
                                                                                                                                      aCCCCGTAGGTAGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCg             >  1:1195859/1‑150 (MQ=11)
                                                                                                                                       ccccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa            >  1:422708/1‑150 (MQ=31)
                                                                                                                                               gTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAAc    >  1:163753/1‑150 (MQ=34)
                                                                                                                                                 cGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGt  >  2:574644/1‑150 (MQ=35)
                                                                                                                                                 |                                                                                                                                                     
TACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGT  >  NC_000913/2304504‑2304798

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