Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,688 A→G 15.1% intergenic (+295/+420) 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,6880AG15.1% 64.5 / 8.5 40intergenic (+295/+420)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base A (12/16);  new base G (3/2);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

CGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAT  >  NC_000913/2304545‑2304836
                                                                                                                                               |                                                                                                                                                    
cGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCgg                                                                                                                                                <  1:1106966/150‑1 (MQ=34)
cGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCgg                                                                                                                                                <  1:1106960/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:1713149/150‑1 (MQ=32)
     aCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGcc                                                                                                                                           <  1:937659/150‑1 (MQ=32)
          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)
                                                                                                                                         aCGCCCAAGCCGGTGGCCTGCTGCGACGCTGGCGCGTCGTCTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGAGAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt       <  2:2669950/150‑1 (MQ=11)
                                                                                                                                           gCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGGAACGCAGGAGCGTcg     >  2:1159066/1‑149 (MQ=255)
                                                                                                                                             ccAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTa   >  2:338031/1‑150 (MQ=37)
                                                                                                                                              cAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt  <  2:1001024/150‑1 (MQ=37)
                                                                                                                                               |                                                                                                                                                    
CGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAT  >  NC_000913/2304545‑2304836

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: