Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,600,220	G→T	D188E (GAC→GAA)	rpoH ←	RNA polymerase, sigma 32 (sigma H) factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,600,220	0	G	T	100.0%	115.8 / NA	33	D188E (GAC→GAA)	rpoH	RNA polymerase, sigma 32 (sigma H) factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/18);  total (15/18)

GTCCAGACCCTGCATCGCGTCGGTCAGACGGTTTGCCGCCTGCTCTTCCCAGTTATCATCTTCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCGTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGA  >  NC_000913/3600085‑3600349                                                                                                                                        |                                                                                                                                  gTCCAGACCCTGCATCGCGTCGGTCAGACGGTTTGCCGCCTGCTCTTCCCAGTTATCATCTTCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAAtcttcgt                                                                                                                               <  1:571819/140‑1 (MQ=255)             cATCGCGTCGGTCAGACGGTTTGCCGCCTGCTCTTCCCAGTTATCATCTTCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGt                                                                                                                   <  1:94612/140‑1 (MQ=255)                cgcgTCGGTCAGACGGTTTGCCGCCTGCTCTTCCCAGTTATCATCTTCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCaa                                                                                                                <  1:1007419/140‑1 (MQ=255)                 gcgTCGGTCAGACGGTTTGCCGCCTGCTCTTCCCAGTTATCATCTTCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCaaa                                                                                                               >  2:406788/1‑140 (MQ=255)                   gTCGGTCAGACGGTTTGCCGCCTGCTCTTCCCAGTTATCATCTTCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAgg                                                                                                             >  2:1001367/1‑140 (MQ=255)                     cGGTCAGACGGTTTGCCGCCTGCTCTTCCCAGTTATCATCTTCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTc                                                                                                           <  2:1015361/140‑1 (MQ=255)                     cGGTCAGACGGTTTGCCGCCTGCTCTTCCCAGTTATCATCTTCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTc                                                                                                           <  1:1001367/140‑1 (MQ=255)                        tCAGCCGGTTTGCCGCCTGCTCTTCCCAGTTATCATCTTCAATGCCGTCGGCAATGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGCCAGGTCAAAGGTCATg                                                                                                        <  2:206220/140‑1 (MQ=255)                                  tGCCGCCTGCTCTTCCCAGTTATCATCTTCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTgccg                                                                                              >  1:917690/1‑140 (MQ=255)                                  tGCCGCCTGCTCTTCCCAGTTATCATCTTCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTgccg                                                                                              <  1:181759/140‑1 (MQ=255)                                                 ccAGTTATCATCTTCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGAt                                                                                   <  2:409472/136‑1 (MQ=255)                                                 ccAGTTATCATCTTCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGAt                                                                                   >  1:409472/1‑136 (MQ=255)                                                  cAGTTATCATCTTCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCAt                                                                              <  1:802981/140‑1 (MQ=255)                                                        tcatcTTCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCacg                                                                        >  2:172000/1‑140 (MQ=255)                                                          atcTTCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCacgta                                                                      <  1:406788/140‑1 (MQ=255)                                                             ttCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCacg                                                                        >  2:897909/1‑135 (MQ=255)                                                             ttCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCacg                                                                        <  1:897909/135‑1 (MQ=255)                                                               cAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCacgtacgt                                                                   <  2:534124/138‑1 (MQ=255)                                                               cAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCacgtacgt                                                                   >  1:534124/1‑138 (MQ=255)                                                                  tGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTg                                                              >  1:139442/1‑140 (MQ=255)                                                                  tGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTg                                                              >  1:468329/1‑140 (MQ=255)                                                                      gTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCAGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTgg                                                          >  1:975931/1‑140 (MQ=255)                                                                             aaGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGcc                                                   <  2:468329/140‑1 (MQ=255)                                                                                    atgatTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCa                                            <  2:139442/140‑1 (MQ=255)                                                                                                 cAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCa                                            >  1:418275/1‑127 (MQ=255)                                                                                                 cAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCa                                            <  2:418275/127‑1 (MQ=255)                                                                                                    aTAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGa                            <  1:1092066/140‑1 (MQ=255)                                                                                                          cACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCAt                      >  1:853297/1‑140 (MQ=255)                                                                                                           aCCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGGTCACGGGCCACCATTTCGACTTCATc                     >  2:845430/1‑140 (MQ=255)                                                                                                                  ccATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCAttt                               >  1:75465/1‑123 (MQ=255)                                                                                                                  ccATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCAttt                               <  2:75465/123‑1 (MQ=255)                                                                                                                     tCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAAc           <  2:1087405/140‑1 (MQ=255)                                                                                                                      cGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAAcc          >  1:857497/1‑140 (MQ=255)                                                                                                                              tgTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTAGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGa  <  1:266436/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  GTCCAGACCCTGCATCGCGTCGGTCAGACGGTTTGCCGCCTGCTCTTCCCAGTTATCATCTTCAATGCCGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCGTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGA  >  NC_000913/3600085‑3600349

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑