Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,214,360	G→A	A438A (GCG→GCA)	rpoD →	RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,214,360	0	G	A	100.0%	76.3 / NA	23	A438A (GCG→GCA)	rpoD	RNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/13);  total (10/13)

ACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACT  >  NC_000913/3214233‑3214494                                                                                                                                |                                                                                                                                       aCCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGctct                                                                                                                            >  1:248348/1‑140 (MQ=255) aCCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGctct                                                                                                                            >  1:271263/1‑140 (MQ=255) aCCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGctct                                                                                                                            >  2:400353/1‑140 (MQ=255)          tGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGCTCTATCGCGGAt                                                                                                                   <  1:400353/140‑1 (MQ=255)            cAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGCTCTATCGCGGATCa                                                                                                                 >  1:340806/1‑140 (MQ=255)                  cTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTAATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGCTCTATCGCGGATCAGgcgcg                                                                                                           <  2:356620/140‑1 (MQ=255)                    tGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGCTCTATCGCGGATCAGGCGCGCa                                                                                                         >  2:13702/1‑140 (MQ=255)                     gACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCATCCGCTCTATCGCGGATCAGGCGCGCAc                                                                                                        <  1:55423/140‑1 (MQ=255)                     gACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGCTCTATCGCGGATCAGGCGCGCAc                                                                                                        <  2:272288/140‑1 (MQ=255)                      aCCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTTGTGGATCCGTCAGGCAATCACCCGCTCTATCGCGGATCAGGCGCGCAcc                                                                                                       <  2:340806/140‑1 (MQ=255)                           aTTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCg                                                                                                  >  1:22365/1‑140 (MQ=255)                            ttCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGCTCTATCGCGGATCAGGCGCGCAc                                                                                                        <  1:33955/133‑1 (MQ=255)                            ttCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGCTCTATCGCGGATCAGGCGCGCAc                                                                                                        >  2:33955/1‑133 (MQ=255)                                 gAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTcc                                                                                            >  2:146711/1‑140 (MQ=255)                                  aaGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCg                                                                                           <  2:70654/140‑1 (MQ=255)                                    ggCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGt                                                                                         <  2:271263/140‑1 (MQ=255)                                    ggCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGt                                                                                         <  2:22365/140‑1 (MQ=255)                                                      aaaGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTgaga                                                                           <  2:364062/136‑1 (MQ=255)                                                      aaaGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTgaga                                                                           >  1:364062/1‑136 (MQ=255)                                                               gATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCt                                                              <  1:146711/140‑1 (MQ=255)                                                                                    ggTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCAATCACCCGCTCGATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGAt                                         <  2:154710/140‑1 (MQ=255)                                                                                                                    gATCCGTCAGGCAATCCCCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCgg         <  2:248348/140‑1 (MQ=255)                                                                                                                           cAGGCAATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCACCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACt  >  2:375411/1‑140 (MQ=255)                                                                                                                                |                                                                                                                                       ACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACT  >  NC_000913/3214233‑3214494

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