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,3730GC95.1% 140.5 / ‑2.7 41intergenic (‑2/+4)caiE/caiDstimulator of CaiD and CaiB enzyme activities/carnitinyl‑CoA dehydratase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base C (20/19);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAGTGGGCCGATGTAGACACCGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCT  >  NC_000913/35250‑35503
                                                                                                                           |                                                                                                                                  
gCGAGTGGGCCGATGTAGACACCGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCttt                                                                                                                     <  1:717787/139‑1 (MQ=255)
  gAGTGGGCCGATGTAGACACCGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTcc                                                                                                                   >  2:757766/1‑139 (MQ=255)
      gggCCGATGTAGACACCGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCacac                                                                                                               >  2:548062/1‑139 (MQ=255)
                  acacCGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCtt                                                                                                   >  1:356909/1‑139 (MQ=255)
                      cGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCg                                                                                               >  1:443587/1‑139 (MQ=255)
                        gCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGc                                                                                             >  2:68379/1‑139 (MQ=255)
                         cTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCa                                                                                            >  2:259258/1‑139 (MQ=255)
                                            aaTCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCt                                                                             <  1:800999/135‑1 (MQ=255)
                                            aaTCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCt                                                                             >  2:800999/1‑135 (MQ=255)
                                                  gACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCAt                                                                   <  1:259258/139‑1 (MQ=255)
                                                       cACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATTGCATCTTcc                                                              >  2:599749/1‑139 (MQ=255)
                                                       cACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTcc                                                              >  2:612387/1‑139 (MQ=255)
                                                       cACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTcc                                                              >  2:300030/1‑139 (MQ=255)
                                                          tGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTc                                                                                                                         >  2:753656/1‑77 (MQ=255)
                                                          tGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTc                                                                                                                         <  1:753656/77‑1 (MQ=255)
                                                          tGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCg                                                                                               <  1:717390/103‑1 (MQ=255)
                                                          tGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCg                                                                                               >  2:717390/1‑103 (MQ=255)
                                                          tGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATcctttga                                                            >  1:438353/1‑132 (MQ=255)
                                                          tGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATcctttga                                                            <  2:438353/138‑7 (MQ=255)
                                                                 gACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGaa                                                    >  1:80596/1‑139 (MQ=255)
                                                                      aCGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGAt                                               >  2:258648/1‑139 (MQ=255)
                                                                          cGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGAt                                           <  2:799433/139‑1 (MQ=255)
                                                                          cGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGAt                                           <  1:150341/139‑1 (MQ=255)
                                                                                        ccGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACAcg                             <  1:599749/139‑1 (MQ=255)
                                                                                          ggAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACAcgcc                           <  2:551535/139‑1 (MQ=255)
                                                                                          ggAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACAcgcc                           >  2:136823/1‑139 (MQ=255)
                                                                                          ggAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACAcgcc                           <  1:134342/139‑1 (MQ=255)
                                                                                          ggAATTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACAcgcc                           >  1:404179/1‑139 (MQ=255)
                                                                                             aTTAACCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCt                        <  1:548062/139‑1 (MQ=255)
                                                                                                aaCCCCTCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCg                     >  1:397862/1‑139 (MQ=255)
                                                                                                      tCAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAAtata               >  1:728409/1‑139 (MQ=255)
                                                                                                       cAAAGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAg              <  1:136823/139‑1 (MQ=255)
                                                                                                          aGGCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGa           <  1:335036/139‑1 (MQ=255)
                                                                                                          aGGCGTAATAGCTCACGCTTGTTAACGGCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGa           <  1:731744/139‑1 (MQ=255)
                                                                                                            gCGTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATa         <  2:362183/139‑1 (MQ=255)
                                                                                                              gTAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACg       <  1:452666/139‑1 (MQ=255)
                                                                                                               tAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGaa                                                    <  2:628539/93‑1 (MQ=255)
                                                                                                               tAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGaa                                                    >  1:628539/1‑93 (MQ=255)
                                                                                                               tAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACAcg                             <  2:445284/116‑1 (MQ=255)
                                                                                                               tAATAGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACAcg                             >  1:445284/1‑116 (MQ=255)
                                                                                                                   aGCTCACGCTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGcttct  <  2:356909/139‑1 (MQ=255)
                                                                                                                           |                                                                                                                                  
GCGAGTGGGCCGATGTAGACACCGGCTCCCACAATCACATCGCCAATCAAGACGGCACTGGGATGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCT  >  NC_000913/35250‑35503

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