Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	934,042	T→G	D273E (GAT→GAG)	ftsK →	cell division DNA translocase FtsK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	934,042	0	T	G	100.0%	91.8 / NA	26	D273E (GAT→GAG)	ftsK	cell division DNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (6/20);  total (6/20)

GCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGC  >  NC_000913/933916‑934164                                                                                                                               |                                                                                                                           gCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGAc                                                                                                                      <  2:397949/133‑1 (MQ=255)  cATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACg                                                                                                                     <  1:497069/133‑1 (MQ=255)    tGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGaa                                                                                                                   <  2:371295/133‑1 (MQ=255)    tGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGaa                                                                                                                   <  2:131987/133‑1 (MQ=255)       aTCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAagag                                                                                                                <  1:140235/133‑1 (MQ=255)       aTCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAagag                                                                                                                <  1:399720/133‑1 (MQ=255)       aTCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAagag                                                                                                                >  1:411340/1‑133 (MQ=255)       aTCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAagag                                                                                                                <  2:411340/133‑1 (MQ=255)       aTCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAagag                                                                                                                <  2:392434/133‑1 (MQ=255)       aTCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAagag                                                                                                                <  2:345866/133‑1 (MQ=255)                tGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAAGAGATTACCTac                                                                                                       <  1:426347/133‑1 (MQ=255)                tGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGGTCTCCGGTAAGCGGATAGAGGATGACGAAGAGATTACCTac                                                                                                       >  2:189356/1‑133 (MQ=255)                             cgcgGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAAGAGATTACCTACACTGCa                                                                                                 >  2:751/1‑126 (MQ=255)                             cgcgGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAAGAGATTACCTACACTGCa                                                                                                 <  1:751/126‑1 (MQ=255)                             cgcgGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAAGAGATTACCTACACTGCACGCGgtg                                                                                          >  2:510459/1‑133 (MQ=255)                                               cgtAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGgacg                                                                        <  2:350502/133‑1 (MQ=255)                                                     cGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTcc                                                                  <  1:510459/133‑1 (MQ=255)                                                                  aaTTCATTCATCCGACGGCGCGCCAAACAGCCGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCaa                                                     <  1:193107/133‑1 (MQ=255)                                                                            tCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACg                                           <  2:148763/133‑1 (MQ=255)                                                                                aTGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGc                                       <  2:128010/133‑1 (MQ=255)                                                                                aTGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGc                                       <  2:360210/133‑1 (MQ=255)                                                                                           aaCAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCGAGAATATGAc                            <  2:57666/133‑1 (MQ=255)                                                                                                  gCTGCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACg                     <  1:500860/133‑1 (MQ=255)                                                                                                     gCGTTGTTCTCCGGTAAGCGGATGGAGGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATc                  <  1:587045/133‑1 (MQ=255)                                                                                                                cGGTAAGCGGATGGAGGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAAc       >  2:43373/1‑133 (MQ=255)                                                                                                                     aGCGGATGGAGGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTgc  >  2:488396/1‑133 (MQ=255)                                                                                                                               |                                                                                                                           GCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGC  >  NC_000913/933916‑934164

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