Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,738,324	C→T	G16S (GGC→AGC)	sdhA ←	succinate dehydrogenase flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,738,324	0	C	T	100.0%	68.9 / NA	21	G16S (GGC→AGC)	sdhA	succinate dehydrogenase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/14);  total (7/14)

GGTCGGCCGAAGCGATGGCGCAGGTGATGCCGCCCTGGGCGGATACGGTGTGCGAACGGGTCGGGAAGACCTTGGTGACTACGGCAGTCTTGTGGCCGCCTTGAGCCAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCCACCGATGATGATGGCGTCGAAGGAAATGGTTGGAATGCTAGCCATGGATCAGATACCCCAGAGAATCTGCACACCCCAGACGAAGTAAGCGAACATCGCAACGCCGCATACCGCCTGGAACAGGAAACG  >  NC_002947/4738183‑4738453                                                                                                                                              |                                                                                                                                  ggTCGGCCGAAGCGATGGCGCAGGTGATGCCGCCCTGGGCGGATACGGTGTGCGAACGGGTCGGGAAGACCTTGGTGACTACGGCAGTCTTGTGGCCGCCTTGAGCCAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCTACCgat                                                                                                                             <  4:139139/148‑1 (MQ=255)      gCCGAAGCGATGGCGCAGGTGATGCCGCCCTGGGCGGATACGGGGTGCGAACGGGTCGGGAAGACCTTGGTGACTACGGCAGTCTTGTGGCCGCCTTGAGCCAGTTGCAGCGCGGCGCGCATGCCGGCGCCGCCGCTACCgatgatga                                                                                                                        <  6:40894/148‑1 (MQ=255)        cGAAGCGATGGCGCAGGTGATGCCGCCCTGGGCGGATACGGTGTGCGAACGGGTCGGGAAGACCTTGGTGACTACGGCAGTCTTGTGGCCGCCTTGAGCCAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCTACCGatgatgatg                                                                                                                      >  1:10989/1‑148 (MQ=255)        cGAAGCGATGGCGCAGGTGATGCCGCCCTGGGCGGATACGGTGTGCGAACGGGTCGGGAAGACCTTGGTGACTACGGCAGTCTTGTGGCCGCCTTGAGCCAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCTACCGatgatgatg                                                                                                                      <  3:124702/148‑1 (MQ=255)                    gcAGGTGATGCCGCCCTGGGCGGATACGGTGTGCGAACGGGTCGGGAAGACCTTGGTGACTACGGCAGTCTTGTGGCCGCCTTGAGCCAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCTACCGATGATGATGGCGTCGAAGGaa                                                                                                          >  5:94483/1‑148 (MQ=255)                      aGGTGATGCCGCCCTGGGCGGATACGGTGTGCGAACGGGTCGGGAAGACCTTGGTGACTACGGCAGTCTTGTGGCCGCCTTGAGCCAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCTACCGATGATGATGGCGTCGAAGGAAAt                                                                                                        <  3:37254/148‑1 (MQ=255)                        gTGATGCCGCCCTGGGCGGATACGGTGTGCGAACGGGTCGGGAAGACCTTGGTGACTACGGCAGTCTTGTGGCCGCCTTGAGCCAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCTACCGATGATGATGGCGTCGAAGGAAATgg                                                                                                      <  2:107186/148‑1 (MQ=255)                                               ggTGTGCGAACGGGTCGGGAAGACCTTGGTGACTACGGCAGTCTTGTGGCCGCCTTGAGCCAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCTACCGATGATGATGGCGTCGAAGGAAATGGTTGGAATGCTAGCCATGGa                                                                                   >  2:51734/1‑144 (MQ=255)                                                                ggAAGACCTTGGTGACTACGGCAGTCTTGTGGCCGCCTTGAGCCAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCTACCGATGATGATGGCGTCGAAGGAAATGGTTGGAATGCTAGCCATGGATCAGATACCCCAGAGAATCTg                                                              <  6:128525/148‑1 (MQ=255)                                                                     aCCTTGGTGACTACGGCAGTCTTGTGGCCGCCTTGAGCCAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCTACCGATGATGATGGCGTCGAAGGAAATGGTTGGAATGCTAGCCATGGATCAGATACCCCAGAGAATCTGCacac                                                         >  6:84496/1‑148 (MQ=255)                                                                       cTTGGTGACTACGGCAGTCTTGTGGCCGCCTTGAGCCAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCTACCGATGATGATGGCGTCGAAGGAAATGGTTGGAATGCTAGCCATGGATCAGATACCCCAGAGAATCTGCACAcc                                                        >  5:66407/1‑147 (MQ=255)                                                                               cTACGGCAGTCTTGTGGCCGCCTTGAGCCAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCTACCGATGATGATGGCGTCGAAGGAAATGGTTGGAATGCTAGCCATGGATCAGATACCCCAGAGAATCTGCACACCCCAGACGaa                                               <  3:20008/148‑1 (MQ=255)                                                                                          ttGTGGCCGCCTTGAGCCAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCTACCGATGATGATGGCGTCGAAGGAAATGGTTGGAATGCTAGCCATGGATCAGATACCCCAGAGAATCTGCACACCCCAGACGAAGTAAGCGAACa                                    >  3:84559/1‑148 (MQ=255)                                                                                             tgGCCGCCTTGAGCCAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCTACCGATGATGATGGCGTCGAAGGAAATGGTTGGAATGCTAGCCATGGATCAGATACCCCAGAGAATCTGCACACCCCAGACGAAGTAAGCGAACATc                                  <  6:33139/147‑1 (MQ=255)                                                                                                          ccAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCTACCGATGATGATGGCGTCGAAGGAAATGGTTGGAATGCTAGCCa                                                                                       >  7:83914/1‑81 (MQ=255)                                                                                                           cAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCTACCGATGATGATGGCGTCGAAGGAAATGGTTGGAATGCTAGCCAt                                                                                      <  8:83914/81‑1 (MQ=255)                                                                                                               tGCAGCGCTGCGCGCATGCCGGCGCCACCGCTACCGATGATGATGGCGTCGAAGGAAATGGTTGGAATGCTAGCCATGGATCAGATACCCCAGAGAATCTGCACACCCCAGACGAAGTAAGCGAACATCGCAACGCCGCATACCGCCt               <  7:103129/148‑1 (MQ=255)                                                                                                                 cAGCGCTGCGCGCATGCCGGCGCCGCCGCTACCGATGATGATGGCGTCGAAGGAAATGGTTGGAATGCTAGCCATGGATCAGATACCCCAGAGAATCTGCACACCCCAGACGAAGTAAGCGAACATCGCAACGCCGCATACCGCCTg              <  3:128965/147‑1 (MQ=255)                                                                                                                     gcTGCGCGCATGCCGGCGCCGCCGCTACCGATGATGATGGCGTCGAAGGAAATGGTTGGAATGCTAGCCATGGATCAGATACCCCAGAGAATCTGCACACCCCAGACGAAGTAAGCGAACATCGCAACGCCGCATACCGCCTGGAACa         <  1:53139/148‑1 (MQ=255)                                                                                                                          gcgcATGCCGGCGCCGCCGCTACCGATGATGATGGCGGCGAAGGAAATGGTTGGAATGCTAGCCATGGATCAGATACCCCAGAGAATCTGCACACCCCAGACGAAGTAAGCGAACATCGCAACGCCGCATACCGCCTGGAACAGGaaa    <  5:40645/148‑1 (MQ=255)                                                                                                                            gcATGCCGGCGCCGCCGCTACCGATGATGATGGCGTCGAAGGAAATGGTTGGAATGCTAGCCATGGATCAGATACCCCAGAGAATCTGCACACCCCAGACGAAGTAAGCGAACATCGCAACGCCGCATACCGCCTGGAACAGGAAACg  <  5:25398/148‑1 (MQ=255)                                                                                                                                              |                                                                                                                                  GGTCGGCCGAAGCGATGGCGCAGGTGATGCCGCCCTGGGCGGATACGGTGTGCGAACGGGTCGGGAAGACCTTGGTGACTACGGCAGTCTTGTGGCCGCCTTGAGCCAGTTGCAGCGCTGCGCGCATGCCGGCGCCGCCGCCACCGATGATGATGGCGTCGAAGGAAATGGTTGGAATGCTAGCCATGGATCAGATACCCCAGAGAATCTGCACACCCCAGACGAAGTAAGCGAACATCGCAACGCCGCATACCGCCTGGAACAGGAAACG  >  NC_002947/4738183‑4738453

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