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_002947	353,067	0	A	G	100.0%	86.7 / NA	26	intergenic (+34/+84)	hisF/cbcV	imidazole 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 (14/12);  total (14/12)

TCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGA  >  NC_002947/352922‑353200                                                                                                                                                  |                                                                                                                                      tCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCCTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGCCCGGGGCTGCTTTGCCGCCCTTCGCggg                                                                                                                                     >  8:96446/1‑148 (MQ=255)   tGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCa                                                                                                                                   >  3:90637/1‑148 (MQ=255)              gCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTccc                                                                                                                        >  3:87868/1‑148 (MQ=255)                          cTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTg                                                                                                            <  3:5421/148‑1 (MQ=255)                           tGGCGGCCAGCATCTTACACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                           <  4:4935/148‑1 (MQ=255)                                   aGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGAGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtg                                                                                                   >  1:20872/1‑148 (MQ=255)                                     cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgtg                                                                                                 <  5:68259/148‑1 (MQ=255)                                          tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCCCGCCCCCCCGGTTGCTGCGCTCTCAGTGGGAGc                                                                                            >  6:145783/1‑148 (MQ=255)                                           ccACTTTGGCGAGTACACAGCGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCCCCCCCGACCACACGGTTGCGGAGATATCTGTGGGAGcc                                                                                           >  2:65226/1‑147 (MQ=255)                                              cTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGAGCCATTGACCGGGGCTGCTTTGCAGCCCGTCGCGGGCACGCCCGCGCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                                        <  3:18992/148‑1 (MQ=255)                                                           acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACACCCGCTCCCACAGGTTTGCGGAGATGTCTGTTGGAGCGGGCGTGCCCGCGAAgg                                                                           >  6:71552/1‑148 (MQ=39)                                                             aGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGACCTTGCTGTGGGAGCGGGCGTGCCCGCGCACGGc                                                                         >  3:172122/1‑148 (MQ=38)                                                               tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcg                                                                       >  2:15196/1‑148 (MQ=255)                                                                        ccAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGCAGGGCCGCGAAGCGGc                                                              >  6:118806/1‑148 (MQ=255)                                                                        ccAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGc                                                              <  4:87868/148‑1 (MQ=255)                                                                         cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCCCCCGCGCCCCCCCGGGGGCGGCGCTCGCTGTGGGCGCGGGCGGGGCCGCGCACGGCCGCGCCGCgggc                                                             >  3:110005/1‑146 (MQ=255)                                                                         cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTCCCACAGGTTTGAGGAGATATCGGTGGGAGCGGGCGGGCCCGCGCAGGGCCGCGCAgccgcc                                                             >  4:99906/1‑148 (MQ=37)                                                                                                   tcgtcgTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTGCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGc                                   >  8:165354/1‑148 (MQ=255)                                                                                                   gcggcgctCGCTGATCCAGTGACCGGGGCTGCTTTGCAGCCCGTCGCGGGCCCGCCCGCTCCCACCGGTTTGCGGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGc                                   <  7:108468/141‑1 (MQ=38)                                                                                                                 gcctttGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGCGATCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTgg                     <  5:162647/144‑1 (MQ=255)                                                                                                                    aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTCCCCCAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGCAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGGTCTTGCCGCGCAGCGCGGGGTaa                  >  2:108983/1‑148 (MQ=38)                                                                                                                     ttGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGttt                                                                                                                >  6:93188/1‑53 (MQ=255)                                                                                                                     ttGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa                 <  3:142312/148‑1 (MQ=255)                                                                                                                      tGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTg                                                                                                               <  5:93188/53‑1 (MQ=255)                                                                                                                          ggggggggctTTGCAGCCCTGCGCGGGCACGCCCGCGCCCACCGGTTGGCGGAGATCGCTGGGGGCGGGGGCGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg            <  1:65226/141‑1 (MQ=255)                                                                                                                          cGGGGCGGCGTGGCAGCCCGGCGCGGGCACGCCCGCTCCCACCGGGTGGCGGCGAGCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg            <  7:96446/148‑1 (MQ=37)                                                                                                                                    ttccAGCCCTGCGCGGGCACGCCGGCTCCCACAGGTTTGCGGAGATCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAAGCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa  <  5:71552/145‑1 (MQ=38)                                                                                                                                                  |                                                                                                                                      TCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGA  >  NC_002947/352922‑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: ‑