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%	187.7 / NA	52	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 (25/27);  total (25/27)

CCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCAT  >  NC_000913/3400639‑3400904                                                                                                                                       |                                                                                                                                    ccAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCa                                                                                                                                 <  2:551681/139‑1 (MQ=255)  cAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCAt                                                                                                                                <  1:1772200/139‑1 (MQ=255)    aTTGCGGCAGCCATCGGTTCTTCAATCAGGGAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATaa                                                                                                                              >  1:793978/1‑139 (MQ=255)    aTTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATaa                                                                                                                              <  2:982292/139‑1 (MQ=255)      tGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAg                                                                                                                            <  2:793978/139‑1 (MQ=255)      tGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAg                                                                                                                            <  2:1924568/139‑1 (MQ=255)         ggCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTg                                                                                                                         <  1:2031125/139‑1 (MQ=255)         ggCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTg                                                                                                                         <  2:341250/139‑1 (MQ=255)            aGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTg                                                                                                                      >  2:2010130/1‑139 (MQ=255)              ccATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCt                                                                                                                    <  1:2010130/139‑1 (MQ=255)                    gTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCAc                                                                                                              <  2:1449412/139‑1 (MQ=255)                    gTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCAc                                                                                                              <  1:303017/139‑1 (MQ=255)                       cttcAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTg                                                                                                           >  2:107611/1‑139 (MQ=255)                              cAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTtgatg                                                                                                    <  1:107611/139‑1 (MQ=255)                                 gAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAg                                                                                                 >  1:1231551/1‑139 (MQ=255)                                      cTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCAc                                                                                                              <  2:413021/121‑1 (MQ=255)                                      cTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCAc                                                                                                              >  1:413021/1‑121 (MQ=255)                                         cACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAg                                                                                         >  1:1873985/1‑139 (MQ=255)                                           cGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAg                                                                                         >  1:47786/1‑137 (MQ=255)                                           cGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAg                                                                                         <  2:47786/137‑1 (MQ=255)                                                    gcgcCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCAtttttt                                                                                 <  1:887843/136‑1 (MQ=255)                                                    gcgcCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCAtttttt                                                                                 >  2:887843/1‑136 (MQ=255)                                                            cgcgGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACgaag                                                                      >  1:1345298/1‑139 (MQ=255)                                                                           tGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                       >  2:1993131/1‑139 (MQ=255)                                                                           tGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                       >  2:1060310/1‑139 (MQ=255)                                                                           tGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                       <  2:1231551/139‑1 (MQ=255)                                                                                     tCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGGGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATCTTTTCAGTCACGAAGCAGTCGGCGATAACGCCCTCTTTCa                                             <  2:1255886/139‑1 (MQ=255)                                                                                            gggTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTg                                                                                                                         <  2:2000434/56‑1 (MQ=255)                                                                                            gggTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTg                                                                                                                         >  1:2000434/1‑56 (MQ=255)                                                                                             ggTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGa                                     <  1:1060310/139‑1 (MQ=255)                                                                                              gTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGaa                                    >  1:477140/1‑139 (MQ=255)                                                                                                 gcgcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                       >  2:230109/1‑117 (MQ=255)                                                                                                 gcgcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                       <  1:230109/117‑1 (MQ=255)                                                                                                  cgcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                       <  1:592753/116‑1 (MQ=255)                                                                                                  cgcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                       >  2:592753/1‑116 (MQ=255)                                                                                                   gcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCAt                                                                                                                                >  2:838047/1‑42 (MQ=255)                                                                                                   gcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCAt                                                                                                                                <  1:838047/42‑1 (MQ=255)                                                                                                   gcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgca                               <  1:13324/139‑1 (MQ=255)                                                                                                              cacaAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGcc                    <  1:678525/139‑1 (MQ=255)                                                                                                                   aCCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCtgtg                                                                                                                 <  2:382417/41‑1 (MQ=255)                                                                                                                   aCCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCtgtg                                                                                                                 >  1:382417/1‑41 (MQ=255)                                                                                                                   aCCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCg               >  2:661278/1‑139 (MQ=255)                                                                                                                      aGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTAc            >  1:1541573/1‑139 (MQ=255)                                                                                                                        aaCGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAg                                                                              >  1:1001184/1‑71 (MQ=255)                                                                                                                        aaCGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAg                                                                              <  2:1001184/71‑1 (MQ=255)                                                                                                                        aaCGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACgg          >  1:1711267/1‑139 (MQ=255)                                                                                                                         aCGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGc         >  2:180759/1‑139 (MQ=255)                                                                                                                           gcgcgGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGccc       <  2:477140/139‑1 (MQ=255)                                                                                                                             gcgGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGTCCGGCGTACGGCCCAg     >  2:1839523/1‑139 (MQ=255)                                                                                                                             gcgGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAg     >  1:509632/1‑139 (MQ=255)                                                                                                                             gcgGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAg     >  2:906463/1‑139 (MQ=255)                                                                                                                                ggCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCAt  <  1:906463/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                    CCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCAT  >  NC_000913/3400639‑3400904

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: ‑