Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,400,773	A→G	R105R (CGT→CGC)	mreB ←	cell wall structural complex MreBCD, actin‑like component MreB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,400,773	0	A	G	100.0%	171.9 / NA	47	R105R (CGT→CGC)	mreB	cell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/26);  total (21/26)

CAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAG  >  NC_000913/3400640‑3400901                                                                                                                                      |                                                                                                                                 cAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCAt                                                                                                                             <  2:643228/139‑1 (MQ=255) cAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCAt                                                                                                                             <  1:618061/139‑1 (MQ=255)   aTTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATaa                                                                                                                           >  2:96827/1‑139 (MQ=255)   aTTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATaa                                                                                                                           >  1:1138573/1‑139 (MQ=255)   aTTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATaa                                                                                                                           >  2:36010/1‑139 (MQ=255)     tGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAg                                                                                                                         <  1:36010/139‑1 (MQ=255)        ggCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTg                                                                                                                      >  2:620963/1‑139 (MQ=255)         gCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGt                                                                                                                     >  1:89499/1‑139 (MQ=255)           aGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTg                                                                                                                   <  1:171845/139‑1 (MQ=255)            gCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGc                                                                                                                  >  2:596295/1‑139 (MQ=255)             ccATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCa                                                                                                                              >  2:1281253/1‑126 (MQ=255)             ccATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCa                                                                                                                              <  1:1281253/126‑1 (MQ=255)                tCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCtgtg                                                                                                              <  1:620963/139‑1 (MQ=255)                  ggTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCa                                                                                                            <  1:487126/139‑1 (MQ=255)                  ggTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCa                                                                                                            <  2:1138573/139‑1 (MQ=255)                   gTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCAc                                                                                                           <  2:371617/139‑1 (MQ=255)                       ttcAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACtt                                                                                                         <  1:768713/137‑1 (MQ=255)                       ttcAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACtt                                                                                                         >  2:768713/1‑137 (MQ=255)                       ttcAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGt                                                                                                       >  2:173526/1‑139 (MQ=255)                             cAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTtgatg                                                                                                 <  2:51874/139‑1 (MQ=255)                             cAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTtgatg                                                                                                 <  2:146849/139‑1 (MQ=255)                                      ttCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTgg                                                                                        <  1:96827/139‑1 (MQ=255)                                                ccAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGTTGGAGCAtttttt                                                                              <  2:881609/139‑1 (MQ=255)                                                   gcgcCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAg                                                                           >  1:619404/1‑139 (MQ=255)                                                                     cGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGa                                                         >  1:1091074/1‑139 (MQ=255)                                                                          tGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                    <  1:596295/139‑1 (MQ=255)                                                                          tGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                    <  2:736962/139‑1 (MQ=255)                                                                            cgcgGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGcc                                                  >  2:140294/1‑139 (MQ=255)                                                                                   ttCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCt                                                                                                                 <  1:507128/69‑1 (MQ=255)                                                                                   ttCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCt                                                                                                                 >  2:507128/1‑69 (MQ=255)                                                                                     cAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCAt                                         <  2:89499/139‑1 (MQ=255)                                                                                           gggTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCg                                   >  2:350796/1‑139 (MQ=255)                                                                                            ggTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGa                                  <  2:1091074/139‑1 (MQ=255)                                                                                             gTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGtt                                                                                                      >  1:877762/1‑70 (MQ=255)                                                                                             gTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGa                                                         <  2:962989/115‑1 (MQ=255)                                                                                             gTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGa                                                         >  1:962989/1‑115 (MQ=255)                                                                                             gTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGCTGCTGTGCACTTGtt                                                                                                      <  2:877762/70‑1 (MQ=255)                                                                                                 cgcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCa                                                                                                            >  2:484712/1‑60 (MQ=255)                                                                                                 cgcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCa                                                                                                            <  1:484712/60‑1 (MQ=255)                                                                                                   ccAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgcag                           <  1:37040/139‑1 (MQ=255)                                                                                                    cAACCGGCACACAAGCCAGACCGCGCGGGCTTGGGCGTTTAAAGCTGTTGCTGTGCACTTGTTGGATGAAGTGCTGGATCTTTTTTTCAGTCACGAAGAAGTCCGCGATAACGCCGTCTTTCATTGGGCGAATGgcagc                          <  2:258194/139‑1 (MQ=255)                                                                                                        cGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAata                      <  2:90805/139‑1 (MQ=255)                                                                                                           cacacaAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTg                   >  1:435951/1‑139 (MQ=255)                                                                                                                  aCCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                    <  1:362160/99‑1 (MQ=255)                                                                                                                  aCCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                    >  2:362160/1‑99 (MQ=255)                                                                                                                  aCCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCg            >  2:748915/1‑139 (MQ=255)                                                                                                                            gcgGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAg  >  1:1121299/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                 CAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAG  >  NC_000913/3400640‑3400901

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