Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 40,278 C→T V37I (GTC→ATC)  caiA ← crotonobetaine reductase subunit II, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091340,2780CT98.0% 159.0 / ‑4.6 51V37I (GTC→ATC) caiAcrotonobetaine reductase subunit II, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (23/27);  total (23/28)
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.97e-01

CCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGG  >  NC_000913/40155‑40411
                                                                                                                           |                                                                                                                                     
ccAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCa                                                                                                                        <  1:97675/139‑1 (MQ=255)
         tccaCACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCaaa                                                                                                               >  1:223520/1‑139 (MQ=255)
         tccaCACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCaaa                                                                                                               >  2:665686/1‑139 (MQ=255)
            acacGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATa                                                                                                            <  1:828907/139‑1 (MQ=255)
                   ggcgAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTcc                                                                                                     >  1:78137/1‑139 (MQ=255)
                              cAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCt                                                                                          <  2:223520/139‑1 (MQ=255)
                                   cccGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCa                                                                                     <  2:165096/139‑1 (MQ=255)
                                     cgcgTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTc                                                                                                                     <  1:366522/105‑1 (MQ=255)
                                     cgcgTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTc                                                                                                                     >  2:366522/1‑105 (MQ=255)
                                     cgcgTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCaaa                                                                                                               <  1:414797/111‑1 (MQ=255)
                                     cgcgTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCaaa                                                                                                               >  2:414797/1‑111 (MQ=255)
                                       cgTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAg                                                                                 >  1:234398/1‑139 (MQ=255)
                                       cgTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAg                                                                                 >  2:603043/1‑139 (MQ=255)
                                        gTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGt                                                                                >  2:86328/1‑139 (MQ=255)
                                                   cGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATAcc                                                                     >  2:559104/1‑139 (MQ=255)
                                                   cGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATAcc                                                                     >  2:157848/1‑139 (MQ=255)
                                                           tCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACaa                                                             <  2:234398/139‑1 (MQ=255)
                                                             aGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCa                                                                                                                        >  1:480356/1‑78 (MQ=255)
                                                             aGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCa                                                                                                                        <  2:480356/78‑1 (MQ=255)
                                                             aGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAAc                                                           <  1:405588/139‑1 (MQ=255)
                                                             aGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAAc                                                           <  1:86328/139‑1 (MQ=255)
                                                                    agcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGt                                                                                                                                  <  2:897618/61‑1 (MQ=255)
                                                                    agcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGt                                                                                                                                  >  1:897618/1‑61 (MQ=255)
                                                                    agcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCa                                                                                                                        <  2:756627/71‑1 (MQ=255)
                                                                    agcagGCTGTCGATACCCACATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCa                                                                                                                        >  1:756627/1‑71 (MQ=255)
                                                                     gcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTg                                                   <  2:1039491/139‑1 (MQ=255)
                                                                      caggcTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTc                                                                                                                     <  2:305374/68‑1 (MQ=255)
                                                                      cagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTc                                                                                                                     >  1:305374/1‑72 (MQ=255)
                                                                                    ccATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCg                                                                                            <  2:487276/83‑1 (MQ=255)
                                                                                    ccATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCg                                                                                            >  1:487276/1‑83 (MQ=255)
                                                                                       tatCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTGGCTGGCCATCAGTTCGCGGATACCGGCGACaaa                                                            >  2:1020857/1‑112 (MQ=255)
                                                                                       tatCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAg                                                                                 <  1:34837/91‑1 (MQ=255)
                                                                                       tatCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAg                                                                                 >  2:34837/1‑91 (MQ=255)
                                                                                       tatCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACaaa                                                            <  1:1020857/112‑1 (MQ=255)
                                                                                       tatCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATttaaattaaa                                 >  2:586672/1‑139 (MQ=255)
                                                                                             ccAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATAc                                                                      <  1:696168/96‑1 (MQ=255)
                                                                                             ccAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATAc                                                                      >  2:696168/1‑96 (MQ=255)
                                                                                              cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCAt                                                                                    >  2:326876/1‑81 (MQ=255)
                                                                                              cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCAt                                                                                    <  1:326876/81‑1 (MQ=255)
                                                                                              cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATc                              <  2:296543/135‑1 (MQ=255)
                                                                                              cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATc                              >  1:296543/1‑135 (MQ=255)
                                                                                                   cTTTGACAAAACGTTCCGGGTAGATGTTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCa                     <  2:689088/139‑1 (MQ=255)
                                                                                                   cTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCa                     <  1:157848/139‑1 (MQ=255)
                                                                                                    tttGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCaa                    <  1:603043/139‑1 (MQ=255)
                                                                                                    tttGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCaa                    <  2:78137/139‑1 (MQ=255)
                                                                                                          aaaaCGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCtt              <  1:665686/139‑1 (MQ=255)
                                                                                                           aaaCGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCCCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTg             >  1:79192/1‑139 (MQ=255)
                                                                                                                  ccGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATAtttt      >  1:182775/1‑139 (MQ=255)
                                                                                                                    gggTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTgg    <  1:170925/139‑1 (MQ=255)
                                                                                                                      gTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTgggg  <  1:948578/139‑1 (MQ=255)
                                                                                                                      gTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTgggg  <  1:831751/139‑1 (MQ=255)
                                                                                                                           |                                                                                                                                     
CCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGG  >  NC_000913/40155‑40411

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