Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	935,625	A→G	Q801R (CAG→CGG)	ftsK →	cell division DNA translocase FtsK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	935,625	0	A	G	100.0%	156.2 / NA	47	Q801R (CAG→CGG)	ftsK	cell division DNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/26);  total (21/26)

CAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATT  >  NC_000913/935489‑935756                                                                                                                                         |                                                                                                                                    caacacCCATTCCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCgcggca                                                                                                                                  <  2:1351806/140‑1 (MQ=255)  aacaacCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCgcggcag                                                                                                                                 <  2:433045/140‑1 (MQ=255)          gTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATca                                                                                                                         <  1:1076828/140‑1 (MQ=255)            tCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATcagc                                                                                                                       >  2:591399/1‑140 (MQ=255)             cccccgcagccgcagTACCAGCACCCACAACAGCCGGTTGCGCCGCAGCAACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATcagca                                                                                                                      <  2:309090/137‑1 (MQ=21)              cgccgcAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGTGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATcagcag                                                                                                                     <  2:283012/140‑1 (MQ=255)              cgccgcAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGTGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATcagcag                                                                                                                     >  1:283012/1‑140 (MQ=255)              cgccgcAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATcagcag                                                                                                                     <  2:125068/140‑1 (MQ=255)              cgccgcAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATcagcag                                                                                                                     >  1:125068/1‑140 (MQ=255)                ccgcagccgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGcc                                                                                                                   >  2:1249962/1‑140 (MQ=255)                ccgcagccgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGcc                                                                                                                   >  1:1355018/1‑140 (MQ=255)                   cagccgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCcaca                                                                                                                <  1:1249962/140‑1 (MQ=255)                     gccgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCca                                                                                                                  >  1:154802/1‑136 (MQ=255)                     gccgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCca                                                                                                                  <  2:154802/136‑1 (MQ=255)                            tATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGt                                                                                                       <  1:1263255/140‑1 (MQ=255)                                 gcagcCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTgcgc                                                                                                  <  1:591399/140‑1 (MQ=255)                                 gcagcCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTgcgc                                                                                                  <  2:673381/140‑1 (MQ=255)                                    gcCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCcac                                                                                               <  2:1435655/140‑1 (MQ=255)                                     ccACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCcaca                                                                                              >  2:613137/1‑140 (MQ=255)                                      cacaACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGTCACAACAGCCAGTTgcgc                                                                                                  <  1:725616/135‑1 (MQ=255)                                      cacaACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGTCACAACAGCCAGTTgcgc                                                                                                  >  2:725616/1‑135 (MQ=255)                                      cacaACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACaa                                                                                             <  2:967936/140‑1 (MQ=255)                                      cacaACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACaa                                                                                             >  2:863760/1‑140 (MQ=255)                                      cacaACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACaa                                                                                             >  2:1067353/1‑140 (MQ=255)                                      cacaACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACaa                                                                                             >  1:967936/1‑140 (MQ=255)                                        caacaGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAAcc                                                                                           <  1:1067353/140‑1 (MQ=255)                                        caacaGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAAcc                                                                                           >  2:286184/1‑140 (MQ=255)                                        caacaGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAAcc                                                                                           >  2:1407486/1‑140 (MQ=255)                                           caGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCGACAACCGGTTGTGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCa                                                                                        <  1:1438678/140‑1 (MQ=255)                                              ccGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAgg                                                                                      >  1:1047637/1‑139 (MQ=255)                                              ccGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAgg                                                                                      <  2:1047637/139‑1 (MQ=255)                                                     cgccgcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCa                                                                                        >  1:560834/1‑130 (MQ=255)                                                     cgccgcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCa                                                                                        <  2:560834/130‑1 (MQ=255)                                                                    atCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTg                                                               <  2:25846/140‑1 (MQ=255)                                                                       agcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCt                                                                            >  2:1305052/1‑124 (MQ=255)                                                                       agcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCt                                                                            <  1:1305052/124‑1 (MQ=255)                                                                           gcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTa                                                        <  1:286184/140‑1 (MQ=255)                                                                           gcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTa                                                        >  1:1287278/1‑140 (MQ=255)                                                                                         ttGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgt                                          <  2:1287278/140‑1 (MQ=255)                                                                                           gcgcCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtcc                                        <  1:863760/140‑1 (MQ=255)                                                                                            cgcCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtccg                                       <  1:962396/140‑1 (MQ=255)                                                                                                  agcagcAATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCAt                                 >  2:672622/1‑140 (MQ=255)                                                                                                  accagcaATATCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCAt                                 <  2:264586/138‑1 (MQ=255)                                                                                                           atcatcAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACcgacg                        <  2:1355018/140‑1 (MQ=255)                                                                                                            tCAGCAGCCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACcgacga                       >  2:1431536/1‑140 (MQ=255)                                                                                                                  gcCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCa                                  <  2:1223348/123‑1 (MQ=255)                                                                                                                  gcCGCAACAACCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCa                                  >  1:1223348/1‑123 (MQ=255)                                                                                                                      caacaaCCGGTTGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGcc             <  2:390333/140‑1 (MQ=255)                                                                                                                                 tGCGCCGCGGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGAtt  >  1:262987/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                    CAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATT  >  NC_000913/935489‑935756

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