Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 353,067 A→G intergenic (+34/+84) hisF → / ← cbcV imidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947353,0670AG100.0% 67.0 / NA 22intergenic (+34/+84)hisF/cbcVimidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/13);  total (9/13)

TGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGA  >  NC_002947/352924‑353200
                                                                                                                                               |                                                                                                                                     
tGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGc                                                                                                                                    <  1:170858/147‑1 (MQ=255)
 ggAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCAc                                                                                                                                  >  3:20472/1‑148 (MQ=255)
   aaGGGCATGCCAGCGCGGGGCTGGCGGCCAGCATCTTCCACTTTGGCGAGGACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGAGCGGCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCAc                                                                                                                                  <  3:43487/146‑1 (MQ=255)
              aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGGCCGCCCCCGCCccca                                                                                                                     >  6:15510/1‑148 (MQ=255)
                             gcgaGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAtct                                                                                                      <  4:145448/145‑1 (MQ=255)
                              ccAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAtctc                                                                                                     <  2:143008/148‑1 (MQ=255)
                               cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTACCACAGGGGTGCTGAGctctgg                                                                                                    >  4:142709/1‑146 (MQ=38)
                                 gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgt                                                                                                  <  3:57215/148‑1 (MQ=255)
                                  cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgtg                                                                                                 <  2:3443/148‑1 (MQ=255)
                                                     tACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGc                                                                                                              >  3:125210/1‑116 (MQ=255)
                                                      acacAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                             <  4:125210/116‑1 (MQ=255)
                                                       cacaGGCCCGGCGGCCATGGCCTAAAGGGCTTCGGGCGGGTTCGGCTTGCGCTCATCCATTCACCGGGGCGGCTTGGCAGCCGTTCGCGGGAACGCCCGCGCCCACAGGTTTGCTGATATCTCTGGGGGAGCGGGCGTGCCCGCGAAg                                                                            <  6:137789/148‑1 (MQ=19)
                                                          aGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGc                                                                         >  4:161317/1‑148 (MQ=255)
                                                                ggAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGaa                                                                   <  2:95919/148‑1 (MQ=255)
                                                                      cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGGGGAGATCTCTGTGGGCGCGGGCGTGCCCGCGACGGGCCGCGCAGCGGcc                                                             >  7:23067/1‑148 (MQ=38)
                                                                                     ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTACCCCAGGTTTGCTGAGATCTCTGTGGGAGCCGGCGTGCCCGCGAAGGGCCCCGAAGCGGCCCCAATTGCTTCAAcc                                              >  3:46988/1‑148 (MQ=39)
                                                                                                      tgctgtGATCCTTGGACGGGGGGTGCTTTCCTGCCCGGGGCGGGCCGGCCGGCACCCACGGGTTTGTGTTGATCGATGTGGGCGGGGGGTGGCCCGCGAATGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagca                             <  7:135471/143‑1 (MQ=255)
                                                                                                                 aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTCCCCCAGGTTTGCTGCGATCTGTGTGGGCGCGGGCGTGCCCGCGCAGGGCCGCGCAGCGGCCCCAATTGCTTCCACCGTGGTTCTTGCCGAGCAGCGCGGGGt                    >  5:54056/1‑146 (MQ=37)
                                                                                                                 aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGGGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCACCCGTGGTTGTTGCCGAGCAGCGCGGGGt                    >  2:106162/1‑146 (MQ=38)
                                                                                                                       cGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGTTTTGCTGAGATGTCTGTGGTAGCGGGCGTGCCCGCGAAGGGGCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGACGAGCAGCGCGTGGTAAAGCTCg            <  2:107668/148‑1 (MQ=38)
                                                                                                                                   gCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa  <  8:86537/146‑1 (MQ=255)
                                                                                                                                   gCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa  <  4:163801/146‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                     
TGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGA  >  NC_002947/352924‑353200

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