Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 35,373 G→C intergenic (‑2/+4) caiE ← / ← caiD stimulator of CaiD and CaiB enzyme activities/carnitinyl‑CoA dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091335,3730GC98.3% 210.5 / ‑6.6 60intergenic (‑2/+4)caiE/caiDstimulator of CaiD and CaiB enzyme activities/carnitinyl‑CoA dehydratase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (28/31);  total (28/32)
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.99e-01

GCAGTGAGGCGAGTGGGCCGATGTAGACACCGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGC  >  NC_000913/35242‑35499
                                                                                                                                   |                                                                                                                              
gCAGTGAGGCGAGTGGGCCGATGTAGACACCGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTaa                                                                                                                         >  2:428896/1‑139 (MQ=255)
   gTGAGGCGAGTGGGCCGATGTAGACACCGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTa                                                                                                                          >  1:72668/1‑135 (MQ=255)
   gTGAGGCGAGTGGGCCGATGTAGACACCGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTa                                                                                                                          <  2:72668/135‑1 (MQ=255)
       ggCGAGTGGGCCGATGTAGACACCGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCtt                                                                                                                  <  1:428896/139‑1 (MQ=255)
           aGTGGGCCGATGTAGACACCGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCca                                                                                                              >  2:1429811/1‑139 (MQ=255)
           aGTGGGCCGATGTAGACACCGGCTCCCACAATCACAGCGCCAATCAAGACGGCACTGGGATGGACAAACGCAGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCca                                                                                                              >  1:923806/1‑139 (MQ=255)
             tGGGCCGATGTAGACACCGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCcaca                                                                                                            <  1:124510/139‑1 (MQ=255)
                gCCGATGTAGACACCGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCg                                                                                                         >  1:835660/1‑139 (MQ=255)
                              cGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTcc                                                                                                               >  2:771622/1‑119 (MQ=255)
                              cGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTcc                                                                                                               <  1:771622/119‑1 (MQ=255)
                              cGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCg                                                                                           >  2:479005/1‑139 (MQ=255)
                              cGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCg                                                                                           >  2:1606617/1‑139 (MQ=255)
                              cGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCg                                                                                           >  1:1321361/1‑139 (MQ=255)
                                gCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGc                                                                                         >  1:842784/1‑139 (MQ=255)
                                  tCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCaa                                                                                       <  1:63498/139‑1 (MQ=255)
                                    ccACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAAc                                                                                     >  2:721187/1‑139 (MQ=255)
                                    ccACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAAc                                                                                     <  2:644097/139‑1 (MQ=255)
                                            acaTCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCgg                                                                             <  1:1606617/139‑1 (MQ=255)
                                             caTCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGc                                                                            <  1:918656/139‑1 (MQ=255)
                                              aTCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGcc                                                                           <  2:923806/139‑1 (MQ=255)
                                                   cAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCa                                                                      <  1:616234/139‑1 (MQ=255)
                                                          gACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAAc                                                                                                                        <  1:223288/82‑1 (MQ=255)
                                                          gACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAAc                                                                                                                        >  2:223288/1‑82 (MQ=255)
                                                            cGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTc                                                                                                                >  1:243436/1‑88 (MQ=255)
                                                            cGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTc                                                                                                                <  2:243436/88‑1 (MQ=255)
                                                               cACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCAt                                                               >  2:1106993/1‑134 (MQ=255)
                                                               cACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCAt                                                               <  1:1106993/134‑1 (MQ=255)
                                                                    ggATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCAt                                                               >  2:53616/1‑129 (MQ=255)
                                                                    ggATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCAt                                                               <  1:53616/129‑1 (MQ=255)
                                                                       tGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAg                                                  <  1:721187/139‑1 (MQ=255)
                                                                         gACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGaa                                                >  1:829386/1‑139 (MQ=255)
                                                                              aCGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTtc                                                                                              <  2:50675/88‑1 (MQ=255)
                                                                              aCGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTtc                                                                                              >  1:50675/1‑88 (MQ=255)
                                                                              aCGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGAt                                           >  2:1239016/1‑139 (MQ=255)
                                                                                 ccGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGa                                        <  1:1551094/139‑1 (MQ=255)
                                                                                  cGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGAt                                       <  2:1263579/139‑1 (MQ=255)
                                                                                   gTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATa                                      <  2:829386/139‑1 (MQ=255)
                                                                                   gTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATa                                      >  1:771496/1‑139 (MQ=255)
                                                                                   gTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATa                                      <  1:479005/139‑1 (MQ=255)
                                                                                    tCGGGTGAACCACCGTAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAg                                     <  1:78634/139‑1 (MQ=255)
                                                                                             ccaccGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAac                            <  2:842784/139‑1 (MQ=255)
                                                                                              caccGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAaca                           <  2:1321361/139‑1 (MQ=255)
                                                                                                  ggAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAg                                                                                >  1:686436/1‑82 (MQ=255)
                                                                                                  ggAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAg                                                                                <  2:686436/82‑1 (MQ=255)
                                                                                                  ggAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACAcgcc                       <  2:1245268/139‑1 (MQ=255)
                                                                                                      ttAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTATCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTg                   >  1:759179/1‑139 (MQ=255)
                                                                                                       tAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCAt                                                               <  1:1613609/94‑1 (MQ=255)
                                                                                                       tAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCAt                                                               >  2:1613609/1‑94 (MQ=255)
                                                                                                        aaCCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCg                 >  1:588164/1‑139 (MQ=255)
                                                                                                             cTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGa                                        <  2:818755/111‑1 (MQ=255)
                                                                                                             cTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGa                                        >  1:818755/1‑111 (MQ=255)
                                                                                                                  aGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGa       <  2:771496/139‑1 (MQ=255)
                                                                                                                     cGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACAcg                         <  1:53181/118‑1 (MQ=255)
                                                                                                                     cGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACAcg                         >  2:53181/1‑118 (MQ=255)
                                                                                                                     cGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATAc    <  2:588164/139‑1 (MQ=255)
                                                                                                                       tAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAAcc                                              >  2:1002368/1‑95 (MQ=255)
                                                                                                                       tAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAAcc                                              <  1:1002368/95‑1 (MQ=255)
                                                                                                                       tAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGc  >  1:1444394/1‑139 (MQ=255)
                                                                                                                               cACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCAt                                                               <  1:97466/70‑1 (MQ=255)
                                                                                                                               cACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCAt                                                               >  2:97466/1‑70 (MQ=255)
                                                                                                                                   |                                                                                                                              
GCAGTGAGGCGAGTGGGCCGATGTAGACACCGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGC  >  NC_000913/35242‑35499

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