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%	128.2 / NA	42	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 (22/20);  total (22/20)

TTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCAGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCC  >  NC_002947/352921‑353206                                                                                                                                                    |                                                                                                                                              ttCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGATTCGCGCGGGATCGTCGTTCGCTGATCCATTGAACGGGGCTGCTTTGCA‑GCCCTTCGCggg                                                                                                                                             >  6:354920/1‑149 (MQ=255)       aaGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTTCA‑GCCCTTCGCGGGccagcc                                                                                                                                       >  8:298311/1‑144 (MQ=255)            cctGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTc                                                                                                                                <  5:311535/147‑1 (MQ=255)             aTGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTcc                                                                                                                               >  8:165446/1‑149 (MQ=255)                 cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCca                                                                                                                             >  7:113361/1‑147 (MQ=255)                         tGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCt                                                                                                                   <  7:167791/149‑1 (MQ=255)                         tGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCt                                                                                                                   <  7:165446/149‑1 (MQ=255)                           cTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTga                                                                                                                 <  8:113361/149‑1 (MQ=255)                                ggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGAtct                                                                                                            >  2:242022/1‑149 (MQ=255)                                    aGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgt                                                                                                        <  1:359015/149‑1 (MQ=255)                                             cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTACCACAGGTTTGCTGAGATCTCTGTGGGAGCggg                                                                                               >  2:68737/1‑149 (MQ=255)                                              aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCACTACCACATGTTTGATGAGATAGATGTGGGAGCGGGc                                                                                              >  5:10579/1‑149 (MQ=38)                                                  tGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGt                                                                                            <  2:40728/147‑1 (MQ=255)                                                   ggCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGc                                                                                          <  4:261836/148‑1 (MQ=255)                                                           cacaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAgg                                                                                 >  6:70316/1‑149 (MQ=255)                                                            acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCCCAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                >  5:122098/1‑149 (MQ=255)                                                                tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATACATTGACCGGCGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCCCACCCACAGGTTTGGGGAGATCTCTTTGGGACCGGGCGTGGCCCCCAAgggaccc                                                                            >  3:273615/1‑145 (MQ=37)                                                                 gCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcgcg                                                                           <  5:70316/149‑1 (MQ=255)                                                                  ccGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcgcg                                                                           >  1:291812/1‑148 (MQ=255)                                                                     gAGGCCAAGGCCTACCTGGCGTCGCGCGGGATCGGCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGc                                                                       <  8:266549/149‑1 (MQ=255)                                                                        gCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCGGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGc                                                                    <  8:23728/149‑1 (MQ=255)                                                                          cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGccc                                                                  <  2:246116/149‑1 (MQ=255)                                                                          cAAGGCCAACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTATGCA‑GCCCTTCGCGGGC‑‑ACCCCCCCTGCCCCACGTTTGCTGAGATGTCAGAGGGAGCGGGCGTGCCCGCGAAGGGCCCCCAAGCGGcac                                                                  >  8:70194/1‑147 (MQ=255)                                                                                  aCATGGCTTCGCGCGGGATCGTCGTACGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCt                                                          >  6:360021/1‑149 (MQ=255)                                                                                     tctcTTCGGGCGGTATGGTGGTTCGCTTATCCATTGACGGGGGCTGCTTTGAC‑GCCTTTCGCGGGC‑‑ACGCCGGCTCCCACTGGTTTGGTTATTTCTCTGGGGTTGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTc                                                        <  2:106604/145‑1 (MQ=19)                                                                                       gCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGGCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCaa                                                      >  1:88556/1‑149 (MQ=255)                                                                                         ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGCGGACCCCCCCAACCACAGGTTTGGTGAGATGTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAACCGGCCCCAATTTCTTCaa                                                      >  1:106604/1‑148 (MQ=37)                                                                                            gcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTg                                                 <  2:291812/148‑1 (MQ=255)                                                                                              gcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCGGAGATATCTGTGGGAGCGGGCGTGCCCGCGCAGGGCCGCGCAGCGGCCCCAATTTCTTGAACCGTGGtt                                              >  3:322366/1‑149 (MQ=37)                                                                                                    tcgtcgTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGGTGAGATCTCTGTGGGAGCGGGCGTCGCCGCGAAGGGCCGCGAAGCGGCCCCAAATGCTTCAACCGTGGTTCTTGcc                                        >  7:190883/1‑149 (MQ=38)                                                                                                     cgtcgtTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACTGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGcc                                        <  3:65961/148‑1 (MQ=255)                                                                                                       tcgtTCGCTTATCCATTGACCGGGGCTGCTTTGCA‑GCCCGTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGag                                     <  7:298311/149‑1 (MQ=255)                                                                                                       tcgtTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGag                                     <  1:242022/149‑1 (MQ=255)                                                                                                            cGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAgcg                                >  1:59679/1‑149 (MQ=255)                                                                                                                     aTTGACCGGGGCTGCTTTGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa                       <  6:122098/149‑1 (MQ=255)                                                                                                                          ccGTGGCGGGTTGGCT‑GCCGTGCGCGGGC‑‑AGGCCCGCTCCCACAGGTTTGCTGAGAGCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                  <  6:10579/149‑1 (MQ=37)                                                                                                                                    tttGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcccc          >  5:319543/1‑147 (MQ=255)                                                                                                                                    tttGCA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa        <  2:88556/149‑1 (MQ=255)                                                                                                                                    tttGCA‑GCCCTTCGCGGGC‑‑ACCCCCGCGCCCCCAGGTTTGCTGAGATATAGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCGTCAACCGTGGTTGTTGGCGAGCAGCGCGGGGTAAAGCTCGGGATaccc          >  4:202845/1‑147 (MQ=25)                                                                                                                                        cA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTAt    <  2:59679/149‑1 (MQ=255)                                                                                                                                        cA‑GCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGGGGCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCCACCGTGGGTCTTGCCGAGCAGCGCGTGGTCAAGCTCGGGATCCCCGAGTAt    >  8:324987/1‑149 (MQ=38)                                                                                                                                           gCCCTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATcc  >  7:102035/1‑149 (MQ=255)                                                                                                                                              cTTCGCGGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                  >  6:63376/1‑130 (MQ=255)                                                                                                                                                    |                                                                                                                                              TTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCA‑GCCCTTCGCAGGC‑‑ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCC  >  NC_002947/352921‑353206

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