Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,305,086 C→G 14.0% intergenic (+693/+22) 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,305,0860CG14.0% 174.2 / 19.0 57intergenic (+693/+22)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base C (28/21);  new base G (4/4);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAAC  >  NC_000913/2304939‑2305231
                                                                                                                                                   |                                                                                                                                                 
gcgcGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATAcc                                                                                                                                                  >  4:75359/1‑149 (MQ=255)
      cTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCg                                                                                                                                            >  4:355034/1‑149 (MQ=255)
         aTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGAt                                                                                                                                         <  2:566771/149‑1 (MQ=255)
         aTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGAt                                                                                                                                         <  1:86329/149‑1 (MQ=255)
         aTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAt                                                                                                                                         >  1:502245/1‑149 (MQ=255)
         aTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAt                                                                                                                                         <  2:385880/149‑1 (MQ=255)
          tCATGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAtg                                                                                                                                        >  3:120745/1‑149 (MQ=255)
                 tACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTa                                                                                                                                 >  3:545538/1‑149 (MQ=255)
                    acCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCa                                                                                                                              >  3:226167/1‑149 (MQ=255)
                     ccGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCAt                                                                                                                             <  1:57894/149‑1 (MQ=255)
                      cGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCAtt                                                                                                                            <  4:316122/149‑1 (MQ=255)
                         tgtgAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACa                                                                                                                         >  2:454836/1‑149 (MQ=255)
                             aaGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGc                                                                                                                     >  2:588325/1‑149 (MQ=255)
                                  cTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAatat                                                                                                                >  4:224914/1‑149 (MQ=255)
                                    ccACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGAt                                                                                                                                         >  1:28499/1‑122 (MQ=255)
                                    ccACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGAt                                                                                                                                         <  2:28499/122‑1 (MQ=255)
                                      aCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGc                                                                                                            <  1:561079/149‑1 (MQ=255)
                                         ccGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCAc                                                                                                         >  4:325131/1‑149 (MQ=255)
                                              ggTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCActtctt                                                                                                    <  2:85321/149‑1 (MQ=255)
                                                ttGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCtttt                                                                                                  <  2:194217/149‑1 (MQ=255)
                                                ttGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCtttt                                                                                                  >  4:670827/1‑149 (MQ=255)
                                                ttGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCtttt                                                                                                  >  3:159526/1‑149 (MQ=255)
                                                            gTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGc                                                                                                            >  1:97784/1‑127 (MQ=255)
                                                            gTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGc                                                                                                            <  2:97784/127‑1 (MQ=255)
                                                            gTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGc                                                                                      <  3:670827/149‑1 (MQ=255)
                                                               aCGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGt                                                                                   <  1:498611/149‑1 (MQ=255)
                                                                           aCATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGt                                                                       >  2:567069/1‑149 (MQ=255)
                                                                             aTCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACgg                                                                                           >  2:295839/1‑127 (MQ=255)
                                                                             aTCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACgg                                                                                           <  1:295839/127‑1 (MQ=255)
                                                                             aTCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGttt                                                                     >  4:587238/1‑149 (MQ=255)
                                                                                tAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCgg                                                                  <  1:369594/149‑1 (MQ=255)
                                                                                 aaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGc                                                                 >  1:365019/1‑149 (MQ=255)
                                                                                    gCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTa                                                              <  1:516400/149‑1 (MQ=255)
                                                                                    gCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTa                                                              <  2:502245/149‑1 (MQ=255)
                                                                                    gCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTa                                                              <  3:75359/149‑1 (MQ=255)
                                                                                      ccGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACt                                                            >  3:10589/1‑149 (MQ=255)
                                                                                                        gCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCtt                                          <  4:590299/149‑1 (MQ=255)
                                                                                                          gACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGt                                        >  3:29774/1‑149 (MQ=255)
                                                                                                            cGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGAtgtg                                                                                                                                      <  4:531087/53‑1 (MQ=17)
                                                                                                            cGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGAtgtg                                                                                                                                      >  3:531087/1‑53 (MQ=17)
                                                                                                            cGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCg                                      <  1:567069/149‑1 (MQ=255)
                                                                                                              cTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTa                                    <  1:261044/149‑1 (MQ=255)
                                                                                                                ggCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGt                                  <  4:420998/149‑1 (MQ=255)
                                                                                                                     gTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAAc                             <  3:333231/149‑1 (MQ=255)
                                                                                                                     gTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAAc                             >  3:390053/1‑149 (MQ=255)
                                                                                                                      tCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAAcc                            >  1:671376/1‑149 (MQ=255)
                                                                                                                          aTCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAg                        >  1:418785/1‑149 (MQ=255)
                                                                                                                              gCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTTTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACt                    >  1:86221/1‑149 (MQ=255)
                                                                                                                              gCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACt                    >  1:358696/1‑149 (MQ=255)
                                                                                                                                      cacaTCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGgt            >  3:435888/1‑149 (MQ=255)
                                                                                                                                         aTCCGGCATGGTATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTa         >  2:106994/1‑149 (MQ=255)
                                                                                                                                         aTCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTa         >  4:576550/1‑149 (MQ=255)
                                                                                                                                         aTCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGATTCGGCGTACTATCTGCTGCCTCCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTa         >  3:276277/1‑149 (MQ=255)
                                                                                                                                          tCCGGCATGGTATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTAc        >  3:384408/1‑149 (MQ=255)
                                                                                                                                           ccGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTTCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACt       >  2:581960/1‑149 (MQ=255)
                                                                                                                                             ggCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGc     >  2:462965/1‑149 (MQ=255)
                                                                                                                                                aTACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAAc  <  1:452248/149‑1 (MQ=255)
                                                                                                                                                aTACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAAc  <  3:355034/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAAC  >  NC_000913/2304939‑2305231

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