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%	178.6 / NA	53	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 (28/25);  total (28/25)

GAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑AGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATC  >  NC_002947/352926‑353205                                                                                                                                               |                                                                                                                                             gAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGc                                                                                                                                       <  2:337075/149‑1 (MQ=255)    gggCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACcccc                                                                                                                                     >  6:224207/1‑148 (MQ=255)         tGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCGTGCCCGcgaag                                                                                                                               <  1:246099/149‑5 (MQ=255)           ccAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGGAGCCCTTCGC‑GGGCACGCCCGCTCCcac                                                                                                                             <  1:230701/149‑1 (MQ=255)                      cTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑t                                                                                                                   >  1:69220/1‑147 (MQ=255)                        ggcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGc                                                                                                                     >  5:133435/1‑144 (MQ=255)                          cggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAGGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑Tgag                                                                                                                <  3:241064/146‑1 (MQ=255)                              cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCca                                                                                                                              <  8:157006/129‑1 (MQ=255)                              cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCca                                                                                                                              >  7:157006/1‑129 (MQ=255)                                 cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTgg                                                                                                      >  7:342313/1‑149 (MQ=255)                                       ccACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGc                                                                                                  <  8:342313/147‑1 (MQ=255)                                        cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑GGAGATGTCTGTGGGAGCgg                                                                                                >  4:74586/1‑148 (MQ=255)                                              ggCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTg                                                                                           <  6:133435/147‑1 (MQ=255)                                                 gAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGGGATCTCTGTGGGAGCGGGCGTGCCcgc                                                                                      >  7:175224/1‑149 (MQ=255)                                                  aGTACCCAGTGCCGGAGCCCAAGGCCTACATGGCTTCGCGCGGGATCGCCGTTCGCTGACCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCcgcg                                                                                     <  1:367247/149‑1 (MQ=38)                                                  aGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCGCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCcgc                                                                                      <  6:103163/148‑1 (MQ=255)                                                       acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCGACAGGGTTGC‑GGAGATCTGTGTGGGAGCGGGCGTGCCCGCGAAgg                                                                                 >  4:95643/1‑148 (MQ=39)                                                       acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCACCCACAGGTTTGC‑AGAGATATCTGTGGGAGCGGGCGTGCCCGCGCAggg                                                                                >  8:23640/1‑149 (MQ=39)                                                         aGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCcac                                                                                                                             >  5:74255/1‑103 (MQ=255)                                                          gTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCcaca                                                                                                                            <  6:74255/103‑1 (MQ=255)                                                            gCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCCCTACCACACGGTTGC‑GGAGATCACTGTGGGAGCGGGCGGGCCCGCCAAGGGccc                                                                             >  6:3412/1‑146 (MQ=38)                                                                   gCCAAGGCCTACATGGCTCCGCGCGGGATCGGCGGTCGCTGATCCATTGCCCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACCGGTTGGC‑TGAGACCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCgg                                                                     <  3:74586/148‑1 (MQ=37)                                                                     cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCCCGCCCCCACCCACACGGGTGG‑GGAGATCGATGTGGGAGCGGGGGTGCCCGCCCAGGGCCCCGAAGCgggc                                                                   >  5:43872/1‑146 (MQ=255)                                                                     cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcg                                                                             >  8:291448/1‑138 (MQ=255)                                                                       aGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCACCCACAGGGTTGG‑AGAGGACGCAGTTGGAGCGGGGGTGGCCCCGACGGGCCGCGCACCGG‑cccc                                                                >  1:218590/1‑148 (MQ=25)                                                                         gCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCGTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAAt                                                             <  2:237920/149‑1 (MQ=255)                                                                              cttGGCTGCGCGGGGGTTCGTGGTTCGTTTATCCATTGCCGGGGGCGGCTTGGCAGCCGTTCGC‑GGGCCCGCCCGCTCCCCCCGGTTTGG‑TGAGTTCTCTGTGGTTGGGGGGGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCtt                                                        <  8:175224/147‑1 (MQ=255)                                                                                tGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTc                                                       <  4:259644/148‑1 (MQ=255)                                                                                   cTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGGACCCCCCCACCCACAGGGTTTG‑TGAGATATCTGTGGGGGCGGGGGGGGCCGCCCACGGCCCCGCAGCGGCCCCCAATTgttgaac                                                     >  6:260321/1‑142 (MQ=255)                                                                                        cgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTg                                                <  5:318578/147‑1 (MQ=255)                                                                                        cgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGt                                              >  7:262290/1‑149 (MQ=255)                                                                                               tcgtcgTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGc                                        >  5:168756/1‑148 (MQ=255)                                                                                                  tcgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCg                                      <  4:160445/147‑1 (MQ=255)                                                                                                   cgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTACCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGag                                    >  3:267118/1‑148 (MQ=255)                                                                                                     ttCGCTTATCCTTTGCCCGGGGGTCCTTGGCTGCCAGGCGC‑GGGCCCCCCGGCTCCCACTGGTTGGC‑TCTCATGTCTGGGGGAGGGGGCTGGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGagcag                                 <  2:218590/149‑1 (MQ=255)                                                                                                     ttCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGagca                                  <  3:182127/148‑1 (MQ=255)                                                                                                                aTTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑GCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTa                        >  7:105765/1‑147 (MQ=255)                                                                                                                aTTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTACCGAGCAGCGCGTGGTaaa                      <  7:286077/149‑1 (MQ=255)                                                                                                                aTTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTACCGAGCAGCGCGTGGTaaa                      <  7:280625/149‑1 (MQ=255)                                                                                                                aTTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTACCGAGCAGCGCGTGGTa                        >  8:280625/1‑147 (MQ=255)                                                                                                                aTTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGGGGTTCTTACCGAGCAGCGCGTGGTa                        >  8:286077/1‑147 (MQ=255)                                                                                                                aTTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCGATGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCCATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTa                        >  1:327231/1‑147 (MQ=255)                                                                                                                    aCCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACCCCCCCTCCCACACCTTTTC‑AGACAACACAGAGGGAGCGGGCGGGCCCGCCAAGCGCCGCGCAGAGG‑CCCCAAATGCTTGAAACCTGGGTCGTGCCGCGCAGCGCGTCGTATAGCt                   >  4:364765/1‑148 (MQ=12)                                                                                                                    aCCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACCCCCCCGCCCACACGGTTGC‑TGGGATATAGGTGGGAGCGGGCGGGGCCGCGAACGGCCGGGCAGCGG‑CCCCCAATGCGTGAACCCGGGGTATTTGCGAGCACCGCGGGGGTAAGCt                   >  4:142343/1‑148 (MQ=255)                                                                                                                     ccGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGGTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGGCCGCGAAGGGCCGCGCAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGGGCAGCGCGTGGTAAAGCTCg                 >  1:93514/1‑149 (MQ=38)                                                                                                                      cGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                 <  6:168756/148‑1 (MQ=255)                                                                                                                      cGGGGCGGCTTTGC‑CCCCTTGGCGGGGCACCCCCGCTCCC‑CCCGTTTGCGTGCTAACTCTTGGGGGGGGGGGGTCCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                 <  3:213967/148‑1 (MQ=255)                                                                                                                       ggggCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGGTTGG‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCggg               >  1:327368/1‑148 (MQ=255)                                                                                                                                 tGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa       <  7:23640/147‑1 (MQ=255)                                                                                                                                   cAGCCCTTCGC‑GGCCCCGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGt     <  5:260321/147‑1 (MQ=255)                                                                                                                                    aGCCCTTCGC‑GGGCCCGCCCCCACCCACAGGGTTGG‑TGAGATCTAGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGCAGCGG‑CCCCAATTGCTTCAACCGGGGTTGTTGCCGCACAGCGCGGGGGAAAGGACGGTGTGCCCGAGGAt   >  4:228461/1‑148 (MQ=19)                                                                                                                                    aGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATc  <  2:69220/149‑1 (MQ=255)                                                                                                                                          tCGC‑GGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTc                  >  4:106028/1‑127 (MQ=255)                                                                                                                                           cgc‑gGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                 <  3:106028/127‑1 (MQ=255)                                                                                                                                               |                                                                                                                                             GAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑AGGCACGCCCGCTCCCACAGGTTTGC‑TGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATC  >  NC_002947/352926‑353205

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