Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,492,438 C→T 7.5% Q621* (CAG→TAG)  PP_2185 → formate dehydrogenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,492,4380CT 7.5% 153.5 / 2.3 53Q621* (CAG→TAG) PP_2185formate dehydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (23/26);  new base T (2/2);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

GTAGGGGTGAACCCGCTGCGCGGGCAGAACAACGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCAGAGCGACCCCAACACCCAGCACGTCACCGCGGCCTTGCTGGCCATGGAGTG  >  NC_002947/2492294‑2492583
                                                                                                                                                |                                                                                                                                                 
gtaggggtGAACCCGCTGCGCGGGCAGAACAACGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAgat                                                                                                                                               >  5:241347/1‑147 (MQ=255)
   ggggtGAACCCGCTGCGCGGGCAGAACAACGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGa                                                                                                                                            <  3:231308/149‑1 (MQ=255)
     ggtGAACCCGCTGCGCGGGCAGAACAACGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGAcc                                                                                                                                          >  4:37435/1‑149 (MQ=255)
      gtGAACCCGCTGCGCGGGCAGAACAACGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGAccc                                                                                                                                         <  4:253405/149‑1 (MQ=255)
        gAACCCGCTGCGCGGGCAGAACAACGTGCAGGGCTCGTGCGTCATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCgg                                                                                                                                       <  5:231130/149‑1 (MQ=255)
        gAACCCGCTGCGCGGGCAGAACAACGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCgg                                                                                                                                       <  8:179518/149‑1 (MQ=255)
              gCTGCGCGGGCAGAACAACGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTgc                                                                                                                                 >  6:175500/1‑149 (MQ=255)
                   gcgGGCAGAACAACGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCa                                                                                                                              >  8:25449/1‑147 (MQ=255)
                   gcgGGCAGAACAACGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCa                                                                                                                              <  7:25449/147‑1 (MQ=255)
                       gCAGAACAACGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTgc                                                                                                                                 <  5:262754/140‑1 (MQ=255)
                       gCAGAACAACGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTgc                                                                                                                                 >  6:262754/1‑140 (MQ=255)
                        cAGAACAACGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCaa                                                                                                                       >  7:63028/1‑149 (MQ=255)
                         aGAACAACGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAAc                                                                                                                      <  6:250321/149‑1 (MQ=255)
                             caacGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTc                                                                                                                           <  3:85272/140‑1 (MQ=255)
                             caacGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTc                                                                                                                           >  4:85272/1‑140 (MQ=255)
                             caaaGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGt                                                                                                                  >  8:147076/1‑149 (MQ=255)
                                cGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGTAGCCCGACCCGGGTCTGCGCATTCCCa                                                                                                                        >  7:135815/1‑140 (MQ=255)
                                cGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGTAGCCCGACCCGGGTCTGCGCATTCCCa                                                                                                                        <  8:135815/140‑1 (MQ=255)
                                cGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGTAGCCCGA‑CCGGGTCTGCGCATTCCCa                                                                                                                        >  3:96359/1‑139 (MQ=255)
                                cGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGTAGCCCGA‑CCGGGTCTGCGCATTCCCa                                                                                                                        <  4:96359/139‑1 (MQ=255)
                                   gCAGGGCTCGTGCGACATGGGCTCGTACCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAgg                                                                                                            >  1:145304/1‑149 (MQ=255)
                                      gggCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAggcgg                                                                                                         <  7:195087/149‑1 (MQ=255)
                                                         tCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCaa                                                                                      <  4:22315/149‑1 (MQ=255)
                                                           gTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAgg                                                                                    >  1:34950/1‑149 (MQ=255)
                                                           gTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAg                                                                                     <  3:37435/148‑1 (MQ=255)
                                                           gTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGAGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAgg                                                                                    >  1:85313/1‑149 (MQ=255)
                                                            ttCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGgc                                                                                   >  8:107618/1‑149 (MQ=255)
                                                                           cccGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAg                                                                      <  5:175500/147‑1 (MQ=255)
                                                                            ccGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGt                                                                   <  8:63028/149‑1 (MQ=255)
                                                                             cGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGt                                                                   >  3:224147/1‑148 (MQ=255)
                                                                                      gCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCg                                                         <  3:250905/149‑1 (MQ=255)
                                                                                        aTGTCTCCAACGAGGGCGGGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGcc                                                       <  6:211437/149‑1 (MQ=255)
                                                                                           tctcCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGc                                                                                           <  4:44376/110‑1 (MQ=255)
                                                                                           tctcCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGc                                                                                           >  3:44376/1‑110 (MQ=255)
                                                                                           tctcCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCag                                                    <  7:147076/149‑1 (MQ=255)
                                                                                             tcCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCagag                                                  >  3:31427/1‑149 (MQ=255)
                                                                                               cAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCaga                                                   <  8:48141/146‑1 (MQ=255)
                                                                                               cAACGAGGGCGAGCCCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCCGGTCTGCGCATTCCCCACATGTTCGAGGCGGCGCTGGATCGCAGCTTCAAGGCGCTGTACTGCCACGGTGAGGACATCGCCCAGAGCg                                                >  7:226302/1‑149 (MQ=255)
                                                                                                   gAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCAGAGCGAccc                                            >  6:40405/1‑149 (MQ=255)
                                                                                                    aGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCAGAGCGAcccc                                           <  6:59418/149‑1 (MQ=255)
                                                                                                      ggCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCAGAGCGACCCCaa                                         >  3:241612/1‑149 (MQ=255)
                                                                                                                      tCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCAGAGCGACCCCAACACCCAGc                                 <  8:172625/141‑1 (MQ=255)
                                                                                                                      tCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCAGAGCGACCCCAACACCCAGc                                 >  7:172625/1‑141 (MQ=255)
                                                                                                                        gAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCAGAGCGACCCCAACACCCAGCACGTCACcg                        <  2:145304/148‑1 (MQ=255)
                                                                                                                          aCAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCAGAGCGACCCCAACACCCAGCACGTCACCGCgg                     <  2:34950/149‑1 (MQ=255)
                                                                                                                              gCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCAGAGCGACCCCAACACCCAGCACGTCACCGCGGCCtt                 >  5:39843/1‑149 (MQ=255)
                                                                                                                                  ggggCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCAGAGCGAACCCAACACCCAGCACGTCACCGCGGCCTTGCTg             >  2:251683/1‑149 (MQ=255)
                                                                                                                                      cgagACCTTGCTGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCAGAGCGACCCCAACACCCAGCACGTCACCGCGGCCTTGCTGGCCa         <  2:85313/146‑1 (MQ=255)
                                                                                                                                        tGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCAGAGCGACCCCAACACCCAGCACGTCACCGCGGCCTTGCTGGCCATg       <  1:91449/149‑1 (MQ=255)
                                                                                                                                            cTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCAGAGCGACCCCAACACCCAGCACGTCACCGCGGCCTTGCTGGCCATGGAGt   <  4:224147/149‑1 (MQ=255)
                                                                                                                                             ttGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACt                                                                           <  4:149107/76‑1 (MQ=255)
                                                                                                                                             ttGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACt                                                                           >  3:149107/1‑76 (MQ=255)
                                                                                                                                             ttGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCAGAGCGACCCCAACACCCAGCACGTCACCGCGGCCTTGCTGGCCATGGAGTg  >  4:217737/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GTAGGGGTGAACCCGCTGCGCGGGCAGAACAACGTGCAGGGCTCGTGCGACATGGGCTCGTTCCCCCACGAGCTGCCCGGCTACCGGCATGTCTCCAACGAGGGCGTGCGCGCCGAGTTCGAACAGGCGTGGGGCGTGACCTTGCAGCCCGACCCGGGTCTGCGCATTCCCAACATGTTCGAGGCGGCGCTGGATGGCAGCTTCAAGGCGCTGTACTGCCAGGGTGAGGACATCGCCCAGAGCGACCCCAACACCCAGCACGTCACCGCGGCCTTGCTGGCCATGGAGTG  >  NC_002947/2492294‑2492583

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