Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,525,382	A→G	T838A (ACG→GCG)	mrcA →	penicillin‑binding protein 1a, murein transglycosylase and transpeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,525,382	0	A	G	100.0%	126.4 / NA	36	T838A (ACG→GCG)	mrcA	penicillin‑binding protein 1a, murein transglycosylase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/21);  total (15/21)

GCCGGAGCAGCCGCTGACGCCGCCACCGGGTATTGTGACGGTGAATATCGATCGCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTC  >  NC_000913/3525234‑3525527                                                                                                                                                     |                                                                                                                                                  gCCGGAGCAGCCGCTGACGCCGCCACCGGGTATTGTGACGGTGAATATCGATCGCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCg                                                                                                                                                 >  2:54925/1‑151 (MQ=255)    ggAGCAGCCGCTGACGCCGCCACCGGGTATTGTGACGGTGAATATCGATCGCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGAcc                                                                                                                                              >  2:62223/1‑151 (MQ=255)      agcagcCGCTGACGCCGCCACCGGGTATTGTGACGGTGAATATCGATCGCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCAt                                                                                                                                            <  2:647060/151‑1 (MQ=255)               tGACGCCGCCACCGGGTATTGTGACGGTGAATATCGATCGCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATa                                                                                                                                    >  2:423557/1‑150 (MQ=255)                 aCGCCGCCACCGGGTATTGTGACGGTGAATATCGATCGCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATg                                                                                                                                 <  1:261294/151‑1 (MQ=255)                 aCGCCGCCACCGGGTATTGTGACGGTGAATATCGATCGCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATg                                                                                                                                 >  1:274464/1‑151 (MQ=255)                              gTATTGTGACGGTGAATATCGATCGCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGa                                                                                                                    <  2:274464/151‑1 (MQ=255)                                aTTGTGACGGTGAATATCGATCGCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGaa                                                                                                                   <  2:69596/150‑1 (MQ=255)                                   gtgACGGTGAATATCGATCGCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTg                                                                                                                <  2:296247/150‑1 (MQ=255)                                      aCGGTGAATATCGATCGCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTc                                                                                                             <  2:521282/150‑1 (MQ=255)                                      aCGGTGAATATCGATCGCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCt                                                                                                            <  2:327354/151‑1 (MQ=255)                                        ggTGAATATCGATCGCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTg                                                                                                           >  1:279731/1‑150 (MQ=255)                                             atatCGATCGCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTaaa                                                                                                     <  1:323703/151‑1 (MQ=255)                                                 cgatcgCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTACACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGGTCTGATTAAAAAgg                                                                                                 >  1:161770/1‑151 (MQ=255)                                                   atcgCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGCACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGgcg                                                                                               >  2:491860/1‑151 (MQ=255)                                                       cagcaCCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTc                                                                                           >  2:45936/1‑151 (MQ=255)                                                                   gTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCa                                                                               <  2:362570/151‑1 (MQ=255)                                                                    ttAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTc                                                                                >  1:364596/1‑149 (MQ=255)                                                                     tAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAg                                                                              <  1:387227/150‑1 (MQ=255)                                                                     tAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAg                                                                              <  2:445353/150‑1 (MQ=255)                                                                             ggtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGAc                                                                     <  2:79468/151‑1 (MQ=255)                                                                             ggtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGAc                                                                     <  1:342302/151‑1 (MQ=255)                                                                                     cAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGc                                                              >  2:353490/1‑150 (MQ=255)                                                                                                      tCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGAt                                             >  1:521166/1‑150 (MQ=255)                                                                                                          cGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGc                                        <  1:301349/151‑1 (MQ=255)                                                                                                                   gCAGCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAAcgcg                               >  1:277555/1‑151 (MQ=255)                                                                                                                      gCCGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTa                            <  2:482849/151‑1 (MQ=255)                                                                                                                       ccGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTa                            <  1:353487/150‑1 (MQ=255)                                                                                                                        cGACACAACAGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATc                          <  2:213976/151‑1 (MQ=255)                                                                                                                              aacaGGCAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCt                     >  2:661857/1‑150 (MQ=255)                                                                                                                                    cAGTGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACaaaa              <  1:339847/151‑1 (MQ=255)                                                                                                                                       tGCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCg           >  2:252308/1‑151 (MQ=255)                                                                                                                                        gCACGAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGt          <  2:265094/151‑1 (MQ=255)                                                                                                                                            gAGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGaa       >  1:325864/1‑150 (MQ=255)                                                                                                                                             aGGTGGGAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGaaaa     <  2:293208/151‑1 (MQ=255)                                                                                                                                                 gggAGCGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTc  <  2:279730/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GCCGGAGCAGCCGCTGACGCCGCCACCGGGTATTGTGACGGTGAATATCGATCGCAGCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTC  >  NC_000913/3525234‑3525527

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