Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,649 C→T 16.6% intergenic (+256/+459) 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,6490CT16.6% 88.9 / 10.2 47intergenic (+256/+459)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base C (16/22);  new base T (2/6);  total (18/29)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

ACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCAT  >  NC_000913/2304550‑2304671
                                                                                                   |                      
aCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTc                        >  2:1440254/1‑100 (MQ=33)
  gccgcATCCGAGATCCAACGCCCGTGCCGGTTGCCTGATGCGACGCGGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTgg                      <  1:1733043/100‑1 (MQ=11)
  gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCgg                      <  2:1289032/100‑1 (MQ=33)
  gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCgg                      <  2:2218949/100‑1 (MQ=33)
  gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCgg                      <  2:640519/100‑1 (MQ=33)
   ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGa                     <  1:1431865/100‑1 (MQ=17)
   ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                     >  2:9992/1‑100 (MQ=33)
   ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                     <  2:582721/100‑1 (MQ=32)
   ccgcAACCGAGATCCAACACCCGAGCCGGTTGCGCGATGCGACGCTGGCGCGTGTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCgg                      <  1:1943477/99‑1 (MQ=11)
    cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATTCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAt                    <  1:1742260/100‑1 (MQ=17)
    cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAt                    <  1:872024/100‑1 (MQ=33)
     gcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATa                   <  1:2800011/100‑1 (MQ=14)
      cATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa                  >  2:2292736/1‑100 (MQ=33)
      cATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATaa                  >  1:1579140/1‑100 (MQ=18)
        tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAgg                >  1:804828/1‑100 (MQ=33)
        ctcGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAgg                <  1:1337733/98‑1 (MQ=12)
         ccGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc               <  2:2525214/100‑1 (MQ=32)
          cGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCGTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGCTAAGGCg              >  1:1586071/1‑100 (MQ=16)
          cGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAggag              >  2:1337738/1‑98 (MQ=12)
          cGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCg              <  2:759196/100‑1 (MQ=33)
          cGACATCCAACGCCCGAGCCCGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGATCCACCCCGTAGGACGGATAAGGCg              >  2:2471095/1‑100 (MQ=255)
           gACATCCACCGCCCGCGGCGGGTGCCTGATGAGACGCTGGCGCGTCTTATCAGGTCTACACCGCTGCGGAGTGCTCCCCCCCGTAGGTCGGATAAGGCgg             <  1:459671/100‑1 (MQ=11)
           gACATCCAACGCCCGAGGCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACTGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg             <  1:1795450/100‑1 (MQ=12)
           gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCgg             <  1:2007617/100‑1 (MQ=14)
           gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg             <  2:538726/100‑1 (MQ=33)
           gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg             <  2:540153/100‑1 (MQ=33)
           gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg             <  1:1061379/100‑1 (MQ=32)
           gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg             <  1:102130/100‑1 (MQ=32)
           gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg             <  1:2090137/100‑1 (MQ=32)
           gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTAGGATAAGGCgg             <  1:2237985/100‑1 (MQ=18)
            aCATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGt            <  2:331460/100‑1 (MQ=33)
            aCATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGt            <  2:1212331/100‑1 (MQ=33)
             cATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGGGTCTGACCAGGGCCACACAGCGGAGAAGTGCGCCACCCCGTAGCTCGGAGAGGGc               >  1:23276/1‑96 (MQ=2)
             cATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCg              <  1:2169917/97‑1 (MQ=14)
             cATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGtt           <  2:2012053/100‑1 (MQ=33)
               tCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTAc         >  1:685743/1‑100 (MQ=18)
               tCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAc         >  1:849989/1‑100 (MQ=33)
               tCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAc         >  1:850108/1‑100 (MQ=33)
                ccAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTACTCCACCCCGTAGGTCGGATAAGGCGGTTAcg        <  1:551923/100‑1 (MQ=17)
                 cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg             <  1:210774/94‑1 (MQ=32)
                 cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgc       >  2:210775/1‑100 (MQ=33)
                 cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAcggcggaaccgc       >  1:2139101/1‑94 (MQ=11)
                  aaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgc       <  1:1742001/99‑1 (MQ=32)
                    cGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgccgc    >  2:527827/1‑100 (MQ=32)
                    cGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGGCGGATAAGGCGGCTAcgccgc    >  2:2223698/1‑100 (MQ=12)
                     gCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCa   >  2:925243/1‑100 (MQ=31)
                      cccGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATGAGGCGGTTAcgccgctt  >  2:2040707/1‑98 (MQ=11)
                      cccGAGCCGGTTGCCTGACGCGACGCTGGCGCGTCTGATCAGGCCGACACCGCTGTGAGGTGCTGCACCCCGTCGGTCGGTTAAGGCGGTTACGCCGCAt  >  1:1762718/1‑100 (MQ=11)
                                                                                                   |                      
ACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCAT  >  NC_000913/2304550‑2304671

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