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% 50.3 / NA 17intergenic (+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/8);  total (9/8)

GAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCA  >  NC_002947/352926‑353209
                                                                                                                                             |                                                                                                                                              
gAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                                          >  6:70685/1‑148 (MQ=255)
    ggCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCggg                                                                                                                                              <  8:24667/140‑1 (MQ=255)
             gcgcGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                             <  3:65795/148‑1 (MQ=255)
                               gCATCTTCCACTTTGGCGAGTACACAGTGCCGGATGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgt                                                                                                           <  2:65414/148‑1 (MQ=255)
                                cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgtg                                                                                                          >  6:64050/1‑148 (MQ=255)
                                    ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGTTTGCTGAGATCTCTGTggggg                                                                                                      >  2:61461/1‑146 (MQ=255)
                                       cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGa                                                                                                       <  4:65865/144‑1 (MQ=255)
                                        aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGAGCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCgg                                                                                                   <  1:66257/147‑1 (MQ=255)
                                                    acacAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                      >  3:163453/1‑148 (MQ=255)
                                                        aGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCGCCCCCCGGTGTGCGGAGATCTCTGTGGGAGCGGGCGTGCCCGCGCACGGc                                                                                  >  7:69575/1‑148 (MQ=255)
                                                                aGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCGCTCCCACAGGGGTGCTGAGATCGCGGTGGGCGCGGGCGGGCGCGCGCCGGGCCGCGCAGc                                                                          >  3:108770/1‑148 (MQ=37)
                                                                              aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCTCCCACAGGTTTTCTGAGCTCTTTGTGGGAGCCGGCGTGCCCGCGAAGGGCCGCGCAGCGGCCCCAATTGcg                                                             >  3:16983/1‑146 (MQ=38)
                                                                               tGGCTGCGCGCGGGATCGTCGTGCGCTGATCCATTGACCGGGGCGGCTTGGCAGCCCTTCGCGGCCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTc                                                           <  8:18319/148‑1 (MQ=38)
                                                                                                gtcgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCg                                          >  4:57732/1‑148 (MQ=255)
                                                                                                                              tttGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATccc              <  5:64050/146‑1 (MQ=255)
                                                                                                                               tggCAGCCCTGCGCGGGCCCGCCCGCGCCCACCGGTTGGCTGAGATCTCGGTGGGAGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa           <  1:61461/146‑1 (MQ=38)
                                                                                                                                 gCAGCCCTGCGCGGGCACGCCCGAGCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa           <  6:134418/146‑1 (MQ=255)
                                                                                                                                        ttCGCGGGCACCCCCGCACCCCCAGGGTTTCTGAGAACTCTGTGGGAGCGGGCGTGCCCGCCCAGGGCCCCCAAGCCGCCCCCATTGCTTAAACCGGGGTTAATGCCCACCACCGCGTGGTAAAGCACGGTAACCCCGCGTATCCcca  >  2:52086/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                              
GAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCA  >  NC_002947/352926‑353209

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: