Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,649 C→T 15.0% 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,6490CT15.0% 144.3 / 15.1 66intergenic (+256/+459)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base T (8/2);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGG  >  NC_000913/2304551‑2304713
                                                                                                  |                                                                
cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCg                                                                 <  2:1688553/100‑1 (MQ=33)
cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCg                                                                 >  1:572924/1‑100 (MQ=33)
cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCg                                                                 <  2:360058/100‑1 (MQ=33)
  ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTGCACCCCGTAGGTCGGa                                                               >  1:330356/1‑100 (MQ=17)
  ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                                                               <  2:2230982/100‑1 (MQ=33)
  ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                                                               <  2:1658027/100‑1 (MQ=33)
  ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                                                               <  2:1173809/100‑1 (MQ=32)
  ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                                                               <  2:98777/100‑1 (MQ=33)
    gcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATa                                                             <  2:510048/100‑1 (MQ=33)
     cATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa                                                            <  2:572922/100‑1 (MQ=32)
     cATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa                                                            <  2:666271/100‑1 (MQ=33)
        ccGACATCCAACGCCCGAGCCGGTTGCCTGATGTGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc                                                         >  1:2129979/1‑100 (MQ=18)
        ccGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc                                                         >  1:1578872/1‑100 (MQ=33)
        ccGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc                                                         <  2:1945028/100‑1 (MQ=32)
        ccGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc                                                         <  2:1496731/100‑1 (MQ=33)
        ccGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc                                                         <  1:186681/100‑1 (MQ=32)
         cGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCg                                                        >  2:414136/1‑100 (MQ=33)
          gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg                                                       <  1:1815511/100‑1 (MQ=33)
          gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg                                                       <  1:1812919/100‑1 (MQ=32)
          gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg                                                       <  2:869633/100‑1 (MQ=33)
          gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg                                                       <  2:987025/100‑1 (MQ=33)
           ggagcacttcacagcggtgtaggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGt                                                      >  1:71211‑M2/48‑100 (MQ=255)
           aCATCCAACGCCCGAGCCGGTTGCCTGATGTGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGt                                                      <  2:2129962/100‑1 (MQ=18)
            cacccAACGCCCGTGCCGGTTGCCTGATGGGACGCTGGCGCGGCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGt                                                      <  2:1738352/96‑1 (MQ=11)
             aTCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGGTa                                                    >  2:1113461/1‑100 (MQ=16)
               ccAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCgggcacg                                                  >  1:2144860/1‑95 (MQ=255)
                cAACGCCTGAGACGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGAGCGGCACCCCGTAGGTCGGATAAGGCGGTTAcgc                                                 >  2:1250240/1‑100 (MQ=11)
                cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTCATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGtt                                                     <  2:1292132/96‑1 (MQ=17)
                 cttcacagcggtgtaggcctgataagacgcgccagcgtcgcATCAGGCCTCCACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCAcgcc                                                >  1:1556033‑M2/42‑100 (MQ=255)
                 aaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgcc                                                <  1:1285337/100‑1 (MQ=33)
                   cGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgccgc                                              <  2:2115634/100‑1 (MQ=32)
                   cGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgccg                                               <  2:1020759/99‑1 (MQ=32)
                           gtgtaggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACa                                      >  1:148602‑M2/32‑100 (MQ=255)
                            tgtaggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACAt                                     >  2:377794‑M2/31‑100 (MQ=255)
                               aggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTa                                  >  1:1094091‑M2/28‑100 (MQ=255)
                               aggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTa                                  >  2:1730485‑M2/28‑100 (MQ=255)
                                ggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCCCCCCGTAGGTTGGATACGGCGGTTCCGCCGCATCCGACATc                                    >  2:2010198‑M2/27‑97 (MQ=255)
                                ggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGAcatccaa                                 >  2:738179‑M2/27‑97 (MQ=255)
                                ggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTaa                                 >  2:589211‑M2/27‑100 (MQ=255)
                                 gcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAAc                                >  2:1743434‑M2/26‑100 (MQ=255)
                                 gcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATACAACATATaa                                 >  2:2050697‑M2/26‑99 (MQ=255)
                                 gcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCAGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACCACTAAc                                >  1:1292139‑M2/26‑100 (MQ=255)
                                  cctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACg                               >  1:1658328‑M2/25‑100 (MQ=255)
                                  cctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACg                               <  2:148601‑M2/76‑1 (MQ=255)
                                    tgataagacgcgcccgcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGTGGTTACGCCGCATCCGACATCTAACGcc                             >  1:1160788‑M2/23‑100 (MQ=255)
                                     gataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGccc                            >  2:1228580‑M2/22‑100 (MQ=255)
                                        aagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACAtccaacgtccgag                         <  1:377796‑M2/82‑7 (MQ=255)
                                              cgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGtt                   >  2:1711498‑M2/13‑100 (MQ=255)
                                                      tcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGAATAAGGCGGTTACGCCGCATCCTAAATCTAACGCCCAAGCCGGTTGCCTGATg           >  2:1900013‑M2/5‑100 (MQ=255)
                                                       ctcATCCGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg           <  2:925589‑M2/96‑1 (MQ=255)
                                                          ATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGAc       >  2:1540477‑M2/1‑100 (MQ=255)
                                                           TCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACg      <  2:1228324‑M2/100‑1 (MQ=255)
                                                           TCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACg      >  1:573716‑M2/1‑100 (MQ=255)
                                                            CAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGc     >  1:616401‑M2/1‑100 (MQ=255)
                                                            CAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGc     >  1:1401314‑M2/1‑100 (MQ=255)
                                                            CAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGc     >  2:2334580‑M2/1‑100 (MQ=255)
                                                            CAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGc     >  1:1973422‑M2/1‑100 (MQ=255)
                                                            CAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATACGACGc     >  2:2007157‑M2/1‑100 (MQ=255)
                                                             AGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCt    >  1:1892827‑M2/1‑100 (MQ=255)
                                                             AGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCAAGCCGGTTGCCTGATGCGACGCt    >  2:396602‑M2/1‑100 (MQ=255)
                                                              GGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGc     >  1:1077054‑M2/1‑98 (MQ=255)
                                                              GGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTg   <  1:1536938‑M2/100‑1 (MQ=255)
                                                              GGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTg   >  1:1749701‑M2/1‑100 (MQ=255)
                                                              GGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTg   >  2:1105557‑M2/1‑100 (MQ=255)
                                                               GCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCt    <  2:921887‑M2/98‑1 (MQ=255)
                                                                CCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGCCATCAAACGCCCAAGCCGGTTGCCTGATGCCCCGCTgg  >  1:1096181‑M2/1‑99 (MQ=255)
                                                                                                  |                                                                
CGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGG  >  NC_000913/2304551‑2304713

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