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% 43.3 / NA 15intergenic (+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 (7/8);  total (7/8)

CAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCC  >  NC_002947/352937‑353208
                                                                                                                                  |                                                                                                                                             
cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCca                                                                                                                               <  7:126194/147‑1 (MQ=255)
  gcgcGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                            <  3:38029/148‑1 (MQ=255)
              cggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTTGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTgaga                                                                                                                 <  5:49288/147‑1 (MQ=255)
                  cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATctct                                                                                                            >  6:80595/1‑148 (MQ=255)
                    gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGAGCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgt                                                                                                          <  1:131109/148‑1 (MQ=255)
                                                   ggAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcgcgca                                                                           >  7:127238/1‑146 (MQ=255)
                                                   ggAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGCTGAGATATCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcgcgca                                                                           >  6:139470/1‑146 (MQ=255)
                                                                  cATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTCCCACACGTGTGGTGAGAAGTCTGTGGGAGCGGGCCTGCCCGCGAACGGCCCAGAAGCGGCCCCAATTGc                                                             >  8:83439/1‑147 (MQ=37)
                                                                                          tCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCCCCCCCACCCACAGGTTTGCTGAGATGTGTGTGGGGGCGGGCGTGCCCCCGCAGGGCCGCGCAGCGGCCCCAAATGGGTAAACCCTGGGTGTTGCCGagcag                                    >  2:56940/1‑148 (MQ=255)
                                                                                                        aCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCGCACCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGCAGGGCCCCGAAGCGGCCCCCATTTCTTCACCCCTGGTTATTGGCGAGCAGCGCGTGGTAAAg                        >  2:55754/1‑146 (MQ=37)
                                                                                                          cGGGGCTGCTTTGCCGCCCTGCGCGGGCCCGCCCGCGCCCACCGGTTTGCTGAGAGCGCTGGGGGAGCGGGCGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                    <  1:56940/148‑1 (MQ=38)
                                                                                                            gggCTGCTTTGCAGCCCTTGGCGGGCACGCCCGCTCCCACAGGTTTGCTTAGAGCTCTGGGGGAGCGGGCTTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGt                  <  2:52073/148‑1 (MQ=38)
                                                                                                               ctgctTTGCAGCCCGGCGCGGGCACGCCCGCTCCCACAGGTTTGCTGATATCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATc               <  4:12395/148‑1 (MQ=39)
                                                                                                                tgctTTGCAGCCCTTCGCGGGCCCGCCCGCTCCCCCAGGGTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTGCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcc              <  3:22482/148‑1 (MQ=38)
                                                                                                                            cTTCGCGGGCACGCCCCCTCCCACAGGTTTGCGGAGATATCTGTGGGAGCGGGCGTGCCCGCGCAGGGCCCCGCAGCGGCCCCAATTGCGTCAACCGTGGTTGTTGGCGAGCCGCGCGTGGGAAAGCTCGGTGTACCCCCGTATcccc  >  2:56470/1‑148 (MQ=25)
                                                                                                                                  |                                                                                                                                             
CAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCC  >  NC_002947/352937‑353208

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