Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,318,096 C→A 8.0% P99T (CCG→ACG)  PP_2039 → acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,318,0960CA 8.0% 142.9 / 2.7 50P99T (CCG→ACG) PP_2039acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (22/24);  new base A (2/2);  total (24/26)
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.06e-01

TGCGAGCAGGGCTTCTACAACATCATGCTGCCCGAGGCGATAGGTGGCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTGTGCCGTGCCGAAAAGGCCGTGTGCTTCGCCCTGACTGAAAGCGAAGCCGGCTCTGATGCAAC  >  NC_002947/2317952‑2318226
                                                                                                                                                |                                                                                                                                  
tGCGAGCAGGGCTTCTACAACATCATGCTGCCCGAGGCGATAGGTGGCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGccc                                                                                                                               <  2:1580715/150‑1 (MQ=255)
         ggCTTCTACAACATCATGCTGCCCGAGGCGATAGGTGGCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTTCGGGCCGCCCCGTATCGGc                                                                                                                      >  1:660473/1‑150 (MQ=255)
                 caacaTCATGCTGCCCGAGGCGATAGGTGGCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCg                                                                                                                        >  1:768890/1‑140 (MQ=255)
                 caacaTCATGCTGCCCGAGGCGATAGGTGGCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCg                                                                                                                        <  2:768890/140‑1 (MQ=255)
                    catcatGCTGCCCGAGGCGATAGGTGGCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCAAGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCaa                                                                                                           <  2:1521967/150‑1 (MQ=255)
                    catcatGCTGCCCGAGGCGATAGGTGGCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCAAGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCaa                                                                                                           >  1:1521967/1‑150 (MQ=255)
                     atcatGCTGCCCGAGGCGATAGGTGGCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAg                                                                                                          >  1:97326/1‑150 (MQ=255)
                        atGCTGCCCGAGGCGATAGGTGGCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTg                                                                                                       >  1:393090/1‑150 (MQ=255)
                          gctgcCCGAGGCGATAGGTGGCGCCGGGTTGATCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGCCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCGGCTCAACGTGGc                                                                                                     >  2:556784/1‑150 (MQ=255)
                            tgcCCGAGGCGATAGGTGGCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGa                                                                                                   <  1:360389/150‑1 (MQ=255)
                               ccGAGGCGATAGGTGGCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGCGCTTTCCGGGCCGTCCCGTGTCGGCCTTCTGGTCAGGTTGGCGAGcc                                                                                                <  2:1252793/150‑1 (MQ=255)
                               ccGAGGCGATAGGTGGCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGcc                                                                                                >  1:647629/1‑150 (MQ=255)
                                 gAGGCGATAGGTGGCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGc                                                                                                 >  1:308993/1‑147 (MQ=255)
                                 gAGGCGATAGGTGGCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGc                                                                                                 <  2:308993/147‑1 (MQ=255)
                                   ggCGATAGGTGGCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGAc                                                                                            >  2:775449/1‑150 (MQ=255)
                                             ggCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAg                                                                                  <  1:1598165/150‑1 (MQ=255)
                                             ggCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAg                                                                                  <  1:775449/150‑1 (MQ=255)
                                             ggCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAg                                                                                  <  2:393090/150‑1 (MQ=255)
                                                gcCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGACGCCCCGTATCGGCCATCt                                                                                                                 <  1:248383/116‑1 (MQ=255)
                                                gcCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGACGCCCCGTATCGGCCATCt                                                                                                                 >  2:248383/1‑116 (MQ=255)
                                                 ccGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCGTCTGTTCAAGGTGGCGAGCCCGGCCCAGAGCGAGGCGt                                                                              >  2:1434364/1‑150 (MQ=255)
                                                       tGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAACGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGc                                                                        <  2:641643/150‑1 (MQ=255)
                                                            gtgtCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCg                                                                   <  2:97326/150‑1 (MQ=255)
                                                              gtCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGgt                                                                 <  1:141769/150‑1 (MQ=255)
                                                              gtCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGgt                                                                 <  1:985640/150‑1 (MQ=255)
                                                              gtCAGACCTGTGTGCGGTCAAGGAAGCCACGCTGCTCAGCTGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTCGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGgt                                                                 <  2:950402/150‑1 (MQ=255)
                                                               tCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGgtg                                                                >  2:1017077/1‑150 (MQ=255)
                                                                 aGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAg                                                                                                          <  1:667076/106‑1 (MQ=255)
                                                                 aGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAg                                                                                                          >  2:667076/1‑106 (MQ=255)
                                                                     cTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTgtgccgt                                                          >  1:1290507/1‑150 (MQ=255)
                                                                     cTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGCCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTgtgccgt                                                          >  2:1186835/1‑150 (MQ=255)
                                                                        tgtgCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGACGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTgtgccgtgcc                                                       >  1:503655/1‑150 (MQ=255)
                                                                          tgCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTGTGCCGTGCCGa                                                     >  2:641431/1‑150 (MQ=255)
                                                                           gCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGAt                                                                                      >  1:676271/1‑116 (MQ=255)
                                                                           gCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGAt                                                                                      <  2:676271/116‑1 (MQ=255)
                                                                             ggTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTg                                                                         >  1:649094/1‑127 (MQ=255)
                                                                             ggTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTg                                                                         <  2:649094/127‑1 (MQ=255)
                                                                               tCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTGTGCCGTGCCGAAAAgg                                                >  2:402410/1‑150 (MQ=255)
                                                                                 aaGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTGTGCCGTGCCGAAAAGGcc                                              >  1:623510/1‑150 (MQ=255)
                                                                                   gccaGCCACGGTGCTCACCGGCCCGATGCTGGCGCGGCGTATGCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTTTTCAAGGTGGCGAGCCCGACCCCGATCGCGGCGTTCCTGCAGCCGGTGTGCCGTGCCGAAAAGGCCgt                                            <  2:660473/147‑1 (MQ=255)
                                                                                             gTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGACGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTGTGCCGTGCCGAAAAGGCCGTGTGCTTCGcc                                  <  2:503655/150‑1 (MQ=255)
                                                                                               gCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTGTGCCGTGCCGAAAAGGCCGTGTGCTTCGCCct                                >  1:466319/1‑150 (MQ=255)
                                                                                                 tCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTGTGCCGTGCCGAAAAGGCCGTGTGCTTCGCCctga                              <  1:1351724/150‑1 (MQ=255)
                                                                                                 tCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTGTGCCGTGCCGAAAAGGCCGTGTGCTTCGCCctga                              >  2:509931/1‑150 (MQ=255)
                                                                                                 tCACCGGCGCGATGCAGGCGCCGCATAGTCTCGGCGAGCTTTCCGGGCAGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTGTGCCGTGCCGAAAAGGCCGTGTGCTTCGCCctga                              <  2:536650/150‑1 (MQ=255)
                                                                                                     cGGCTCGATGCTGGCGCCGACCATTCTTGGCCAGATTTCCGGGCCGCCCCGGATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTGTGCCGTGCCGAAAAGGCCGTGTGCTTCGCCctgactga                          >  2:1252440/1‑150 (MQ=255)
                                                                                                         tCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTGTGCCGTGCCGAAAAGGCCGTGTGCTTCGCCCTGACTGAAAGc                      <  2:466319/150‑1 (MQ=255)
                                                                                                             tGCTGGCGCCGCATAGTCTCGGCGCGCGTTCCGGGCCGCCCCGTATCGGCCATCTGTTGACGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTGTGCCGTGCCGAAAAGGCCGTGTGCTTCGCCCTGACTGAAAGCGAAg                  <  1:556784/150‑1 (MQ=255)
                                                                                                              gCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTGTGCCGTGCCGAAAAGGCCGTGTGCTTCGCCCTGACTGAAAGCGAAGc                 <  1:1170716/150‑1 (MQ=255)
                                                                                                                             tctcGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTGTGCCGTGCCGAAAAGGCCGTTTGCTTCGCCCTGACTGAAAGCGAAGCCGGCTCTGATGCAAc  <  2:1290507/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                  
TGCGAGCAGGGCTTCTACAACATCATGCTGCCCGAGGCGATAGGTGGCGCCGGGTTGAGCGTGTCAGACCTGTGTGCGGTCAAGGAAGCCACGGTGCTCACCGGCTCGATGCTGGCGCCGCATATTCTCGGCGAGCTTTCCGGGCCGCCCCGTATCGGCCATCTGTTCAAGGTGGCGAGCCCGACCCAGATCGAGGCGTTCCTGCAGCCGGTGTGCCGTGCCGAAAAGGCCGTGTGCTTCGCCCTGACTGAAAGCGAAGCCGGCTCTGATGCAAC  >  NC_002947/2317952‑2318226

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