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%	167.2 / NA	46	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 (20/26);  total (20/26)

ACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCA  >  NC_000913/3400638‑3400903                                                                                                                                        |                                                                                                                                   aCCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGgcgc                                                                                                                                 >  1:1416769/1‑139 (MQ=255)   cAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCAt                                                                                                                               <  2:806483/139‑1 (MQ=255)   cAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCAt                                                                                                                               <  2:70915/139‑1 (MQ=255)     aTTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATaa                                                                                                                             >  1:1182286/1‑139 (MQ=255)          ggCAGCCATCGGTTCTTCAATCAGGAAGACTTCACTGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTg                                                                                                                        <  2:974397/139‑1 (MQ=255)          ggCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTg                                                                                                                        <  2:1182286/139‑1 (MQ=255)             aGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTg                                                                                                                     >  2:1102147/1‑139 (MQ=255)                     gTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCAc                                                                                                             <  1:1082260/139‑1 (MQ=255)                         ttcAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGt                                                                                                         >  1:22645/1‑139 (MQ=255)                              tCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTtgat                                                                                                    >  2:716512/1‑139 (MQ=255)                               cAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTtgatg                                                                                                   <  1:825792/139‑1 (MQ=255)                                   aaGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGt                                                                                               <  1:716512/139‑1 (MQ=255)                                        ttCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCAc                                                                                                             >  1:223961/1‑120 (MQ=255)                                        ttCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCAc                                                                                                             <  2:223961/120‑1 (MQ=255)                                          cACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAg                                                                                        >  1:1343248/1‑139 (MQ=255)                                            cGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCa                                                                                                              <  1:1048187/115‑1 (MQ=255)                                            cGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCa                                                                                                              >  2:1048187/1‑115 (MQ=255)                                            cGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCa                                                                                      <  2:1416769/139‑1 (MQ=255)                                                                            tGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAAGCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                      <  1:966743/139‑1 (MQ=255)                                                                            tGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                      <  1:1505450/139‑1 (MQ=255)                                                                            tGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                      <  2:1343248/139‑1 (MQ=255)                                                                              cgcgGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCa                                                                                                              >  2:383873/1‑81 (MQ=255)                                                                              cgcgGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCa                                                                                                              <  1:383873/81‑1 (MQ=255)                                                                                             gggTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCAT‑AAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGa                                    >  2:1494925/1‑139 (MQ=255)                                                                                               gTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGaa                                   >  1:862461/1‑139 (MQ=255)                                                                                                tCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAAt                                  <  2:592094/139‑1 (MQ=255)                                                                                                   cgcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgc                               >  1:276294/1‑139 (MQ=255)                                                                                                    gcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCAtttttt                                                                                <  1:1456246/89‑1 (MQ=255)                                                                                                    gcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCAtttttt                                                                                >  2:1456246/1‑89 (MQ=255)                                                                                                    gcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgca                              >  1:1216061/1‑139 (MQ=255)                                                                                                    gcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATG‑AGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgcag                             <  2:22645/139‑1 (MQ=255)                                                                                                     ccAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgcag                             <  2:862461/139‑1 (MQ=255)                                                                                                     ccAACCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgcag                             <  2:1216061/139‑1 (MQ=255)                                                                                                        aCCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCtg                                                                                                                  >  1:513640/1‑51 (MQ=255)                                                                                                        aCCGGCACACAAACCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCtg                                                                                                                  <  2:513640/51‑1 (MQ=255)                                                                                                                    aCCAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCg              >  1:909074/1‑139 (MQ=255)                                                                                                                      cAGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTa            <  1:881173/139‑1 (MQ=255)                                                                                                                       aGAACGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTAc           <  1:1102147/139‑1 (MQ=255)                                                                                                                         aaCGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCAtt                                          <  1:601332/106‑1 (MQ=255)                                                                                                                         aaCGCGCGGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCAtt                                          >  2:601332/1‑106 (MQ=255)                                                                                                                             cgcgGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTtga                                                                                                     <  1:275749/43‑1 (MQ=255)                                                                                                                             cgcgGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTtga                                                                                                     >  2:275749/1‑43 (MQ=255)                                                                                                                              gcgGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGcc                                                    >  2:1533446/1‑91 (MQ=255)                                                                                                                              gcgGGCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGcc                                                    <  1:1533446/91‑1 (MQ=255)                                                                                                                                gggCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCAtttttt                                                                                >  2:1496276/1‑61 (MQ=255)                                                                                                                                gggCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCAtttttt                                                                                <  1:1496276/61‑1 (MQ=255)                                                                                                                                gggCTTGGGCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCa  >  2:1332665/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                   ACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCA  >  NC_000913/3400638‑3400903

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