Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,318,784 G→T 9.5% Q342H (CAG→CAT)  PP_2039 → acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,318,7840GT 9.5% 142.7 / 4.7 49Q342H (CAG→CAT) PP_2039acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (25/19);  new base T (3/2);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCG  >  NC_002947/2318641‑2318929
                                                                                                                                               |                                                                                                                                                 
gCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATCACCCCGGCGGCAATATCCGAACGCAGGCCCCGATGTGCAAGCTGGTCGTTACGGAGACGGCCTTCCGCATCGCCGATAGCGCCGTGCAGATaca                                                                                                                                              >  1:366057/1‑149 (MQ=255)
     cGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACAcggcg                                                                                                                                         <  5:230196/149‑1 (MQ=255)
           gCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCg                                                                                                                                   >  2:46188/1‑149 (MQ=255)
                ctgcACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGgct                                                                                                                               <  4:91639/148‑1 (MQ=255)
                     cGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCACACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGTCGGGCTGCTGCg                                                                                                                         >  7:120472/1‑149 (MQ=255)
                      gcagcaTGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGt                                                                                                                        <  7:301488/149‑1 (MQ=255)
                        agcaTGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTgg                                                                                                                      <  4:211324/149‑1 (MQ=255)
                              gTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAAccc                                                                                                                >  2:257871/1‑149 (MQ=255)
                                  aCGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGgtg                                                                                                            >  5:53837/1‑149 (MQ=255)
                                   cGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGgtg                                                                                                            <  1:46188/148‑1 (MQ=255)
                                   cGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGgtg                                                                                                            <  4:145376/148‑1 (MQ=255)
                                     ccACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGgtggag                                                                                                         >  6:299032/1‑149 (MQ=255)
                                      cACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGgtggagt                                                                                                        >  8:7159/1‑149 (MQ=255)
                                      cACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGgtggagt                                                                                                        <  5:299032/149‑1 (MQ=255)
                                       aCGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCtgctgc                                                                                                                          <  7:280788/130‑1 (MQ=255)
                                       aCGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCtgctgc                                                                                                                          >  8:280788/1‑130 (MQ=255)
                                         ggCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGgtggagtgga                                                                                                     <  3:351404/149‑1 (MQ=255)
                                               ggggAATGACGCCGTCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTTTTCGTGTCGGAGCCGGCCTTCCGCAGCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCGGCGTGGCAACCCGGTGGAGTGGATATTcc                                                                                               <  7:301403/146‑1 (MQ=255)
                                               ggTGAATGACGCCGGCGCCAATATCCGCACCCAGGCCCCGATGTGCAAGCGGTTCGTTTCGGCTACGGCCTTCCGCATCGCGGTTGGCGCCTTGCATTTACACGGGGGTGCTGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTcc                                                                                               <  1:257871/149‑1 (MQ=255)
                                                 tGAATGACGCCAGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGt                                                                                             >  4:147216/1‑149 (MQ=255)
                                                       aCGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGaca                                                                                       >  7:193441/1‑149 (MQ=255)
                                                               ggcAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTatgat                                                                               >  5:31781/1‑149 (MQ=255)
                                                                         gCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATtctc                                                                     <  6:53837/149‑1 (MQ=255)
                                                                            cccAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCAcc                                                                  <  7:7159/149‑1 (MQ=255)
                                                                              cAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCg                                                                 >  7:325486/1‑148 (MQ=255)
                                                                              cAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCg                                                                 >  5:256740/1‑148 (MQ=255)
                                                                                           gCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATc                                                   >  2:84811/1‑149 (MQ=255)
                                                                                            cAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATc                                                   >  4:299638/1‑148 (MQ=255)
                                                                                               gCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAgc                                               >  3:4655/1‑149 (MQ=255)
                                                                                                  gTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAgcgc                                             >  2:300541/1‑148 (MQ=255)
                                                                                                          cgttgtcGGCGTTCCGCTTCGCGGATGGCGCCGTGCAGATACACTGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGTAGTGGTTATTCTGTGCGTGACGTATGTTGCGGATTCTCACCTGTACCAGCGAGATCCAGCGCATCACGAt                                     <  2:366057/142‑1 (MQ=255)
                                                                                                            gagcCGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAAACCGGGGGAGTGCAGATTCCGTGCGACACGTATGATGCGCATTCTCACCGGACCCACCGAGATCCAGCGCAACCCCAACGc                                  >  7:36300/1‑149 (MQ=255)
                                                                                                               aCGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGcc                                 >  7:85894/1‑147 (MQ=255)
                                                                                                                   ccTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGc                           >  6:18033/1‑149 (MQ=255)
                                                                                                                    cTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCg                          >  6:23423/1‑149 (MQ=255)
                                                                                                                           aTCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCAt                   >  6:44180/1‑149 (MQ=255)
                                                                                                                               ccGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCt               >  6:86544/1‑149 (MQ=255)
                                                                                                                                 gATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGa             <  6:340651/149‑1 (MQ=255)
                                                                                                                                  gcaccgcCGTGTATCTGCAAGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAg            >  8:175470/5‑149 (MQ=255)
                                                                                                                                  aTCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCTCAACACCATCGCCAAGGTCGTGCTCATGCCTGAg            <  8:124947/149‑1 (MQ=255)
                                                                                                                                   tCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGc           <  3:147216/149‑1 (MQ=255)
                                                                                                                                     ccgcCGTGTATCTGCACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGc                 >  1:23386/2‑141 (MQ=255)
                                                                                                                                     ccgcCGTGTATCTGCACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGc                 <  2:23386/140‑1 (MQ=255)
                                                                                                                                      cgcCGTGTATCTGCACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAAt        >  1:37551/1‑149 (MQ=255)
                                                                                                                                      cgcCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAAt        >  3:231621/1‑149 (MQ=255)
                                                                                                                                       gcCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGc                 >  5:347534/1‑139 (MQ=255)
                                                                                                                                       gcCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGc                 <  6:347534/139‑1 (MQ=255)
                                                                                                                                             cAGTTACACGGCGGTGCCGGGCGGCGGCGGGGCTACCCGGTGGTGTGGCTATTCCGTGCGCCCCGTTTGATGCGCATTCTCTCCGGAACCAGCGAGATCCAGCGCAACACCATCTCCAAGGGCGTGCTCATGCCTGAGCAATGAcgcg  <  8:36300/148‑1 (MQ=255)
                                                                                                                                             cAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGAcgcg  <  8:85894/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCG  >  NC_002947/2318641‑2318929

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