Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 37,524 T→G S101R (AGC→CGC)  caiC ← crotonobetaine/carnitine‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927337,5240TG100.0% 94.2 / NA 26S101R (AGC→CGC) caiCcrotonobetaine/carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (15/11);  total (15/11)

GGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCA  >  NZ_CP009273/37390‑37655
                                                                                                                                      |                                                                                                                                   
gggAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTcc                                                                                                                                 >  1:142941/1‑139 (MQ=255)
      cgCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGca                                                                                                                            >  1:210524/2‑138 (MQ=255)
            aTCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAgg                                                                                                                     <  2:262621/139‑1 (MQ=255)
                 tCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAggc                                                                                                                    <  1:371701/135‑1 (MQ=255)
                                                      ctgAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGc                                                                           <  2:377772/139‑1 (MQ=255)
                                                                 gATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGc                                                                 >  1:1872/1‑138 (MQ=255)
                                                                        aGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGAt                                                          >  1:166811/1‑138 (MQ=255)
                                                                         ggATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCgg                                                                                                       <  1:367073/92‑1 (MQ=255)
                                                                         ggATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCgg                                                                                                       >  2:367073/1‑92 (MQ=255)
                                                                          gATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATaaa                                                       <  2:417530/139‑1 (MQ=255)
                                                                               aaTTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCg                                                  >  2:188637/1‑139 (MQ=255)
                                                                                    cgcACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAg                                             >  1:376633/1‑139 (MQ=255)
                                                                                    cgcACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAg                                             >  2:371676/1‑139 (MQ=255)
                                                                                    cgcACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAg                                             >  1:383821/1‑139 (MQ=255)
                                                                                    cgcACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAg                                             >  1:166544/1‑139 (MQ=255)
                                                                                           gTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGACCCAGCAAAGGATAATTTCCGGCCGGTTGtggg                                      >  1:86472/1‑136 (MQ=255)
                                                                                              aCCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCgaga                                    <  2:420350/138‑1 (MQ=255)
                                                                                                     ggCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAg                                             <  2:99705/122‑1 (MQ=255)
                                                                                                     ggCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAg                                             >  1:99705/1‑122 (MQ=255)
                                                                                                        aCGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAAc                         >  2:415390/1‑139 (MQ=255)
                                                                                                               gCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCg                  >  2:354311/1‑139 (MQ=255)
                                                                                                                    tttGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCttt             >  2:154351/1‑139 (MQ=255)
                                                                                                                               cACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCa  <  2:166811/139‑1 (MQ=255)
                                                                                                                               cACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCa  <  2:166544/139‑1 (MQ=255)
                                                                                                                               cACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCa  <  1:212213/139‑1 (MQ=255)
                                                                                                                               cACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCa  <  1:188637/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                   
GGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCA  >  NZ_CP009273/37390‑37655

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