Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,680 T→C 21.8% intergenic (+287/+428) 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_0009132,304,6800TC21.8% 63.8 / 13.1 38intergenic (+287/+428)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base T (9/19);  new base C (3/5);  total (12/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGC  >  NC_000913/2304531‑2304829
                                                                                                                                                     |                                                                                                                                                     
tAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCt                                                                                                                                                       <  1:2257461/150‑1 (MQ=37)
  ggTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTaa                                                                                                                                                     <  1:1098865/150‑1 (MQ=40)
  ggTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTaa                                                                                                                                                     <  1:1098801/150‑1 (MQ=40)
   gTCAGACAAGGCGGTCACGCCGAATCCGACATACAAAGCCCGAGCCGGTTGCCTTATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTGCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAAc                                                                                                                                                    <  2:107267/150‑1 (MQ=14)
    tCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACg                                                                                                                                                   >  2:484117/1‑150 (MQ=39)
    tCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACg                                                                                                                                                   <  1:2063769/150‑1 (MQ=39)
              cGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCgg                                                                                                                                         <  1:1106960/150‑1 (MQ=34)
              cGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCgg                                                                                                                                         <  1:1106966/150‑1 (MQ=34)
               gggcACGCCGCCTCCTCCCTCCCCCGCCCGTGCCGGTGCCCTGCGGCGCCGCGGGCGCGTCTTATCTGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGt                                                                                                                                        <  1:2719394/147‑1 (MQ=11)
                gTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGTTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCCGtt                                                                                                                                       >  1:306906/1‑150 (MQ=255)
                 tCACGCCGAATCCGACATCAAACGACCGAGCCGGTTGCCTGTTGCGACGCTGGAGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTg                                                                                                                                      <  2:484723/150‑1 (MQ=11)
                   aCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGcc                                                                                                                                    <  1:937659/150‑1 (MQ=32)
                   aCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGcc                                                                                                                                    <  1:1713149/150‑1 (MQ=32)
                   aCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTAGGATAAGGCGGTTACGCCGCATCCAACATCt                                                                                                                                                       >  1:1590847/1‑131 (MQ=17)
                        gccaccGACAACCACCGCCAGACCCTGTTGCCTGCGGCGACGCTGGCCCGAGCTATCAGGCCTACACCGCTGTGAAGTGCTGCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                               <  2:1852310/146‑1 (MQ=11)
                              gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGc                                                                                                                         >  2:2321274/1‑150 (MQ=32)
                              gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCGACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGc                                                                                                                         >  1:1822230/1‑150 (MQ=17)
                               aCATCCAACGCCCGTGCCGGTTGCCTGTTTCTACGCTGGCGCGTCTTATCATGCCTACACCGCTGTGAAGTGCTGCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCt                                                                                                                        <  2:2495133/150‑1 (MQ=11)
                               aCATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCt                                                                                                                        <  2:2495132/150‑1 (MQ=17)
                                    cctcGCCCGAGCCGGTGGCCCGCTGCGTCGCGGGCGCGTCTGCGCCGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCCTCCGACATCGACCGCCCCAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                                                   <  2:2486754/147‑1 (MQ=11)
                                    cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                                                   <  1:1411968/150‑1 (MQ=32)
                                         cccGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGACtt                                                                                                              >  2:1643524/1‑150 (MQ=11)
                                                                                                                 tAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATcc                                      >  1:2301744/1‑150 (MQ=32)
                                                                                                                       ggATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCTAGCCGGTTGCCTGATGCGACGCTGGAGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTc                                <  2:1187345/150‑1 (MQ=25)
                                                                                                                       ggATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTc                                <  2:1187774/150‑1 (MQ=35)
                                                                                                                        gacaaGGCGGTCACGCCGCATCCGACATCCAACGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTcc                               <  2:2617470/147‑1 (MQ=0)
                                                                                                                              gCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCg                         <  2:1822230/150‑1 (MQ=37)
                                                                                                                              gCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGGCCGAGCCg                         >  2:284539/1‑150 (MQ=21)
                                                                                                                              gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCg                         >  2:1215081/1‑150 (MQ=12)
                                                                                                                              gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCg                         >  1:474554/1‑150 (MQ=12)
                                                                                                                                   tACGCCGTTTCCTACATCGAACGCCCTAGCCGTTTGCGTGATTCTAATCTGGCTCGTCTTATGTTGCCTAAACCGCTGTGAAGTGCTGCACCCCGTAGGTATGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGc                    <  1:1153141/150‑1 (MQ=11)
                                                                                                                                   tACGACGCATCCTAGATCCTAAGCCCTAGCCGGTGGCGTTATGCGAATCTGGCTCGTCTTTTCTGGCCTAAACCGCTGTGAAGTGCTGCACCCCGTAGGTTTGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGc                    <  1:2114956/150‑1 (MQ=11)
                                                                                                                                    aCGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTTTGATGTGCTCCACCCCTTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGcc                   <  1:484117/150‑1 (MQ=255)
                                                                                                                                    aCGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGcc                   <  2:2313797/150‑1 (MQ=37)
                                                                                                                                     cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCt                  <  1:507964/150‑1 (MQ=21)
                                                                                                                                       ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGtt                      >  2:1520022/1‑144 (MQ=21)
                                                                                                                                       ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGtt                      <  1:1520022/144‑1 (MQ=21)
                                                                                                                                         gcATCCGACATCCACCGCCCCAGCCGGTTGCCGGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTGGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATg              <  2:261003/150‑1 (MQ=11)
                                                                                                                                               gACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCGAACGGCCGAGCCGGTTGCCTGAGGCGgcgc        >  2:2276645/1‑150 (MQ=255)
                                                                                                                                                   tCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTGGGATAAGGCGGTTACGCCGCATCCGACATCCAACGGCCGAGCCGGTTGCCTGAGGCGACGCTGgc    >  2:1530065/1‑150 (MQ=255)
                                                                                                                                                     tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGGCCGAGCCGGTTGCCTGATGCGACGCTGgcgc  >  1:1140336/1‑150 (MQ=21)
                                                                                                                                                     tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCATACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTGACGCCGCATCCGACATCCATCGTCCGAGCCGGTTGCCTGATGCGACGCTGgcgc  >  1:687479/1‑150 (MQ=255)
                                                                                                                                                     tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAGCAGGCCTACACCGCTGGGAAGAGCGCCACCCCGTAGGTGGGATCAGGCGGTTACGCCGCATCCGACATCCAAAGTCCGAGCCTGATGCCAGAGGCGACGCTGgccc  >  2:892380/1‑148 (MQ=11)
                                                                                                                                                     |                                                                                                                                                     
TAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGC  >  NC_000913/2304531‑2304829

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