Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,786,954	A→G	Q39Q (CAA→CAG)	envC →	murein hydrolase activator EnvC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,786,954	0	A	G	100.0%	174.9 / NA	50	Q39Q (CAA→CAG)	envC	murein hydrolase activator EnvC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/20);  total (30/20)

GAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAA  >  NC_000913/3786821‑3787090                                                                                                                                      |                                                                                                                                         gAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCaaa                                                                                                                                    <  1:162851/140‑1 (MQ=255)       cGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTAtt                                                                                                                              <  1:678935/140‑1 (MQ=255)        gATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTc                                                                                                                             <  1:499745/140‑1 (MQ=255)         aTTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCa                                                                                                                            >  1:573271/1‑140 (MQ=255)         aTTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCa                                                                                                                            >  1:795971/1‑140 (MQ=255)                 cATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACa                                                                                                                    >  1:309664/1‑140 (MQ=255)                 cATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACa                                                                                                                    >  1:828761/1‑140 (MQ=255)                  aTGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACAt                                                                                                                   >  1:195216/1‑140 (MQ=255)                  aTGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACAt                                                                                                                   <  1:814/140‑1 (MQ=255)                    gACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATcg                                                                                                                 <  2:795971/140‑1 (MQ=255)                       acGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGACCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATcgccg                                                                                                              >  2:768094/1‑140 (MQ=255)                          ggCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATcg                                                                                                                 >  2:51671/1‑134 (MQ=255)                          ggCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATcg                                                                                                                 >  2:744200/1‑134 (MQ=255)                          ggCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATcg                                                                                                                 <  1:51671/134‑1 (MQ=255)                          ggCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATcg                                                                                                                 <  1:744200/134‑1 (MQ=255)                           gCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCaaa                                                                                                          >  1:863253/1‑140 (MQ=255)                               tGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAagag                                                                                                      >  1:311149/1‑140 (MQ=255)                                 aaaCCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAgcg                                                                                                    <  1:465750/140‑1 (MQ=255)                                      gcgcAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCgg                                                                                               <  2:214478/140‑1 (MQ=255)                                       cgcAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGt                                                                                              >  2:845725/1‑140 (MQ=255)                                         cAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTAc                                                                                            >  2:432905/1‑140 (MQ=255)                                                  aaTCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGc                                                                                   <  1:263420/140‑1 (MQ=255)                                                         cccATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGcaacaaca                                                                            >  2:620106/1‑140 (MQ=255)                                                                  tACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGcaacaaca                                                                            >  1:577750/1‑131 (MQ=255)                                                                  tACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGcaacaaca                                                                            <  2:577750/131‑1 (MQ=255)                                                                   aCGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGc                                                                  <  1:845725/140‑1 (MQ=255)                                                                    cGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGcc                                                                 >  1:135956/1‑140 (MQ=255)                                                                     gCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCt                                                                <  2:309664/140‑1 (MQ=255)                                                                     gCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCt                                                                <  2:828761/140‑1 (MQ=255)                                                                          cGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCg                                                           >  1:495534/1‑140 (MQ=255)                                                                                     cTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGaaa                                                <  2:178574/140‑1 (MQ=255)                                                                                     cTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGaaa                                                >  1:178574/1‑140 (MQ=255)                                                                                      tGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGaaaa                                               >  2:68704/1‑140 (MQ=255)                                                                                       ggCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGaaaaa                                              >  1:532719/1‑140 (MQ=255)                                                                                        gCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAg                                             >  1:242623/1‑140 (MQ=255)                                                                                             ttgttgTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGga                                        >  1:39435/1‑140 (MQ=255)                                                                                                  gtgCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCa                                           >  1:408476/1‑132 (MQ=255)                                                                                                  gtgCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCa                                           <  2:408476/132‑1 (MQ=255)                                                                                                  gtgCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGaagaag                                   >  1:823408/1‑140 (MQ=255)                                                                                                  gtgCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGaagaag                                   >  2:76192/1‑140 (MQ=255)                                                                                                  gtgCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGaagaag                                   >  1:8513/1‑140 (MQ=255)                                                                                                        cTTTTCCGCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGAtct                             >  1:26240/1‑140 (MQ=255)                                                                                                               gCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGc                      <  1:768094/140‑1 (MQ=255)                                                                                                               gCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGc                      <  1:432905/140‑1 (MQ=255)                                                                                                               gCCCACGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGc                      <  2:760765/140‑1 (MQ=255)                                                                                                                   aCGCGGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCAcc                  >  1:809793/1‑140 (MQ=255)                                                                                                                    cgcgGATGAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCAccc                 >  2:793446/1‑140 (MQ=255)                                                                                                                           gAGCGTGACCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCt          >  2:507879/1‑140 (MQ=255)                                                                                                                                  aCCAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGaa   <  2:823408/140‑1 (MQ=255)                                                                                                                                   ccAGCTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGaaa  >  2:211806/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                         GAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGGAAGAAGCGATCTCTGAAGCCACCCGTAAGCTGCGCGAAA  >  NC_000913/3786821‑3787090

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