Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 37,396 T→C A143A (GCA→GCG caiC ← crotonobetaine/carnitine‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927337,3960TC100.0% 100.5 / NA 29A143A (GCA→GCGcaiCcrotonobetaine/carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/16);  total (13/16)

GGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCG  >  NZ_CP009273/37282‑37531
                                                                                                                  |                                                                                                                                       
ggTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGc                                                                                                                 <  1:543294/139‑1 (MQ=255)
                 ccgTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCAGGTTGTTGGTTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATc                                                                                                >  1:228644/1‑139 (MQ=255)
                   gTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCtt                                                                                              >  1:82493/1‑139 (MQ=255)
                                 tAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGt                                                                                <  1:491091/139‑1 (MQ=255)
                                        ggTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                                                         <  1:13242/139‑1 (MQ=255)
                                        ggTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                                                         <  2:438324/139‑1 (MQ=255)
                                        ggTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                                                         <  1:630935/139‑1 (MQ=255)
                                             aTAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGa                                                                    <  2:47174/139‑1 (MQ=255)
                                                  acaAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGaa                                                               <  2:182957/139‑1 (MQ=255)
                                                    aaGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCa                                                                                                  <  2:573842/102‑1 (MQ=255)
                                                    aaGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCa                                                                                                  >  1:573842/1‑102 (MQ=255)
                                                    aaGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAAtt                                                             >  2:608879/1‑139 (MQ=255)
                                                      ggTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTgc                                                           >  2:591811/1‑139 (MQ=255)
                                                      ggTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTgc                                                           >  1:576653/1‑139 (MQ=255)
                                                       gTGGCAGGTTGTTGATGTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTgcg                                                          <  2:509575/139‑1 (MQ=255)
                                                           cAGGTTGTTGATTTTTCAGTTGCGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACt                                                      >  1:98797/1‑139 (MQ=255)
                                                           cAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACt                                                      >  1:391619/1‑139 (MQ=255)
                                                           cAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACt                                                      >  1:412629/1‑139 (MQ=255)
                                                                  ttgATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCAcc                                               <  2:313865/139‑1 (MQ=255)
                                                                  ttgATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCAc                                                >  2:75575/1‑138 (MQ=255)
                                                                         ttCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAgg                                         >  1:341876/1‑138 (MQ=255)
                                                                                 aGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACact                                                                       <  1:586140/100‑3 (MQ=255)
                                                                                 aGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACact                                                                       >  2:586140/1‑98 (MQ=255)
                                                                                           cTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCa                      <  1:362180/139‑1 (MQ=255)
                                                                                           cTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCa                      <  1:504137/139‑1 (MQ=255)
                                                                                                gCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATc                 <  2:576653/139‑1 (MQ=255)
                                                                                                   atcatcAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCAc              >  1:104720/1‑139 (MQ=255)
                                                                                                             ggAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCt    <  1:16179/139‑1 (MQ=255)
                                                                                                               aaGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCg  <  2:104720/139‑1 (MQ=255)
                                                                                                                  |                                                                                                                                       
GGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCG  >  NZ_CP009273/37282‑37531

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