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%	107.8 / NA	33	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 (17/16);  total (17/16)

GGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGT  >  NC_002947/352925‑353202                                                                                                                                                |                                                                                                                                        ggAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCa                                                                                                                                     >  7:37350/1‑147 (MQ=255)  gAAGGGCAAGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCAc                                                                                                                                    >  4:46099/1‑147 (MQ=255)   aaGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACg                                                                                                                                   >  5:251043/1‑147 (MQ=255)   aaGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACg                                                                                                                                   >  4:255408/1‑147 (MQ=255)          tGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCCCGCCCGCTc                                                                                                                            >  5:151804/1‑147 (MQ=255)            ccAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTccc                                                                                                                          >  7:11830/1‑147 (MQ=255)            ccAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTccc                                                                                                                          >  6:203825/1‑147 (MQ=255)                cgcgGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAg                                                                                                                      <  6:251043/147‑1 (MQ=255)                    gTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGttt                                                                                                                  <  1:115855/147‑1 (MQ=255)                    gTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCCGCTTCGCGGCCC‑TTCGCGGGCACGCCCGCTCCCACAgagat                                                                                                                  >  6:264603/1‑143 (MQ=38)                        tGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGTTTGc                                                                                                                >  5:55960/1‑145 (MQ=255)                         ggcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGttt                                                                                                                  <  7:64913/142‑1 (MQ=255)                         gccgcccAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATAGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCACCCC‑TTCTCGGGCACGCCACCTCCCACACGTTTGcag                                                                                                              >  8:64913/6‑144 (MQ=255)                             gcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAtc                                                                                                         <  4:261263/147‑1 (MQ=255)                                                 cGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGcc                                                                                     <  4:133824/147‑1 (MQ=255)                                                      acacAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGa                                                                                <  1:214872/147‑1 (MQ=255)                                                         caGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGTTTGGGGAGATCTGTGTGGGAGCGGGGGTGCCCGCGAAgg                                                                             >  2:186123/1‑147 (MQ=38)                                                                     ccAAGGCCTACATGGCTTCGCGCGGGAGCGGCGTTCGCTGATCCATTCCCCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCCCAGGTTTGCTGCGATCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCgg                                                                 <  6:151804/147‑1 (MQ=37)                                                                                                tcgtcgTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCttt                                      >  2:28587/1‑146 (MQ=255)                                                                                                    cgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGa                                  >  1:38191/1‑147 (MQ=255)                                                                                                     gtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGag                                 <  3:190249/147‑1 (MQ=255)                                                                                                       tCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGTTTTCTGAGATATCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagca                               >  1:4083/1‑147 (MQ=255)                                                                                                       tCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagca                               <  5:203825/147‑1 (MQ=255)                                                                                                         gCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagcag                              <  3:46099/146‑1 (MQ=255)                                                                                                               ccATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                        >  7:70305/1‑146 (MQ=255)                                                                                                                 aTTGACCGGGGCTGCTTTGCAGCCCGTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGCAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTgg                       >  1:3188/1‑146 (MQ=255)                                                                                                                       cGGGGCTGCTTTGCAGCGC‑TTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCGCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTc               <  5:212150/147‑1 (MQ=255)                                                                                                                       cGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTc               <  8:70305/147‑1 (MQ=255)                                                                                                                       cGGGGCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTc               <  2:3188/147‑1 (MQ=255)                                                                                                                         gggCTGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACCGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCgg             <  5:10765/147‑1 (MQ=255)                                                                                                                         gggCGGCTTTGCAGCCC‑TTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGTTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCgg             <  3:255408/147‑1 (MQ=255)                                                                                                                               ctTTGCAGCCC‑TTCGCGGGCACGCCCCCTCCCACAGGGTTGCGGAGATATCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTTCTTAAACCGTGGTTATTTCCGAGCAGCGCCTGGTAAAGCTAGGtttgcc       >  4:235695/1‑142 (MQ=25)                                                                                                                                     aGCCC‑TGCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCTCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGt  <  3:93649/146‑1 (MQ=255)                                                                                                                                                |                                                                                                                                        GGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCC‑TTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGT  >  NC_002947/352925‑353202

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: ‑