Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 36,231 T→C intergenic (‑69/+40) caiD ← / ← caiC carnitinyl‑CoA dehydratase/putative crotonobetaine/carnitine‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091336,2310TC95.2% 210.4 / ‑1.6 63intergenic (‑69/+40)caiD/caiCcarnitinyl‑CoA dehydratase/putative crotonobetaine/carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (28/32);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCT  >  NC_000913/36139‑36327
                                                                                            |                                                                                                
ggTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGatta                                                                                           >  1:968940/1‑100 (MQ=255)
ggTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGatta                                                                                           >  1:1146477/1‑100 (MQ=255)
   cAGATGTAATGATTCACTCATTTCTTTTCTCCATTTTTGCTTTTTCATGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGAttatta                                                                                        <  2:1059206/100‑1 (MQ=255)
   cAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGAttatta                                                                                        <  1:236340/100‑1 (MQ=255)
    aGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAg                                                                                       <  2:78195/100‑1 (MQ=255)
      aTGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCt                                                                                     >  2:728797/1‑100 (MQ=255)
      aTGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCt                                                                                     <  1:170122/100‑1 (MQ=255)
       tGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTg                                                                                    >  2:334354/1‑100 (MQ=255)
        gTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTgg                                                                                   >  2:1250596/1‑100 (MQ=255)
         tAATGATTCACTCATTCCTTTTCTCCATTTTGGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGc                                                                                  <  2:114704/100‑1 (MQ=255)
           aTGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAg                                                                                <  2:877997/100‑1 (MQ=255)
              aTTCACTCATTCCTTTTCTCCATTTTTTCTTTTTCTTTGACGACAAAAAAAAAAAAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGt                                                                             <  2:1251588/100‑1 (MQ=255)
              aTTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGt                                                                             <  2:1237754/100‑1 (MQ=255)
               ttCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTa                                                                            >  2:901625/1‑100 (MQ=255)
                   cTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTc                                                                        >  1:978946/1‑100 (MQ=255)
                   cTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTc                                                                        >  2:395111/1‑100 (MQ=255)
                     cATTCCTTTTCTCCATTTTTGCTTTTTCGGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTccc                                                                      <  2:231591/100‑1 (MQ=255)
                     cATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTccc                                                                      <  2:47699/100‑1 (MQ=255)
                         ccTTTTCTCCATTTTTTTTTTTTTTTTTACGACAACATCCCTGCAAAAAATATATTTTGTTTTAGAGTGTGAATATTAGCTTTCAGGGTATTTctcttc                                                                   >  2:995495/1‑97 (MQ=255)
                         ccTTTTCTCCATTTTTTTTTTTTTTTTGACGACAACATCCCTGCAAAAAATGTATATTGTTTTAGAGCGTGATTATTAGCTGTCAGGGTAGTTCCct                                                                     >  2:434684/1‑97 (MQ=255)
                           tttaCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTATAGAGCGTGATTATAAGCTGGCAGGGTAGTTCCCTGCTGt                                                                >  1:97834/1‑100 (MQ=255)
                             ttCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGAACTACCctgccagct                                                              <  1:825214/100‑6 (MQ=255)
                              tctcCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTc                                                             >  2:1416973/1‑100 (MQ=255)
                               ctcCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCa                                                            <  2:1515845/100‑1 (MQ=255)
                               ctcCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCa                                                            >  1:1037115/1‑100 (MQ=255)
                                      tttGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCatttattt                                                     <  1:933913/100‑1 (MQ=255)
                                       ttGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTc                                                    <  2:1037390/100‑1 (MQ=255)
                                        tGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCa                                                   >  2:227783/1‑100 (MQ=255)
                                          cTTTTTCAGGGACGAGAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGa                                                 <  1:1550328/100‑1 (MQ=255)
                                          cTTTTTCAGGGACGACAACATCCCTGCAAAAACTGCATATTGTTTTAGAGCGTGTTTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAg                                                  >  2:65931/1‑99 (MQ=255)
                                          cTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGa                                                 <  2:513721/100‑1 (MQ=255)
                                           tttttCCGGGACAACAACATAAAAGCAAAAAATGTATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAg                                                  <  2:869130/98‑1 (MQ=255)
                                             tttCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttc                                              <  1:895812/100‑1 (MQ=255)
                                               tCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttctt                                            <  1:395112/100‑1 (MQ=255)
                                                 aGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttctttc                                          >  2:1488261/1‑100 (MQ=255)
                                                 aGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttctttc                                          >  1:1141054/1‑100 (MQ=255)
                                                  gggACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttctttct                                         >  2:452631/1‑100 (MQ=255)
                                                  gggACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttctttct                                         <  1:233076/100‑1 (MQ=255)
                                                  gggACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttctttct                                         <  2:1049697/100‑1 (MQ=255)
                                                   ggACGACAACCTCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTa                                        <  2:463877/100‑1 (MQ=255)
                                                   ggACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTa                                        <  2:1146470/100‑1 (MQ=255)
                                                   ggACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTa                                        >  2:952843/1‑100 (MQ=255)
                                                            caTCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTc                               >  2:463923/1‑100 (MQ=255)
                                                            caTCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTc                               >  2:461928/1‑100 (MQ=255)
                                                            caTCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTc                               >  1:630543/1‑100 (MQ=255)
                                                              tCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTccc                             >  2:577273/1‑100 (MQ=255)
                                                              tCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTccc                             <  2:792850/100‑1 (MQ=255)
                                                                  tGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAg                         <  1:1241215/100‑1 (MQ=255)
                                                                   gCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGc                        >  1:1233789/1‑100 (MQ=255)
                                                                    cAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCa                       <  1:461931/100‑1 (MQ=255)
                                                                        aaaTGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTa                   <  1:1262714/100‑1 (MQ=255)
                                                                            gCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCTCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTg               >  2:378750/1‑100 (MQ=255)
                                                                             cATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTgg              >  1:135173/1‑100 (MQ=255)
                                                                               tatTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCa            <  1:728800/100‑1 (MQ=255)
                                                                                  tGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGAt         >  1:1120298/1‑100 (MQ=255)
                                                                                   gTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATc        <  2:851658/100‑1 (MQ=255)
                                                                                    ttttAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCt       >  1:593738/1‑100 (MQ=255)
                                                                                       tAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTCTTTTCCCCGAGCAATTACGTGGCAGATCtttt    >  1:1225167/1‑100 (MQ=255)
                                                                                       tAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCtttt    <  1:707492/100‑1 (MQ=255)
                                                                                        agagCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTc   <  1:334355/100‑1 (MQ=255)
                                                                                        agagCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTc   <  1:452633/100‑1 (MQ=255)
                                                                                        agagCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTc   <  1:1004590/100‑1 (MQ=255)
                                                                                         gagCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCt  >  2:1357247/1‑100 (MQ=255)
                                                                                            |                                                                                                
GGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCT  >  NC_000913/36139‑36327

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