Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	39,378	C→A	D586Y (GAT→TAT)	mecA ←	penicillin binding protein 2 prime

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	39,378	0	C	A	100.0%	64.6 / NA	19	D586Y (GAT→TAT)	mecA	penicillin binding protein 2 prime
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (10/9);  total (10/9)

AATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCAT  >  CP000730/39241‑39512                                                                                                                                          |                                                                                                                                       aaTAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATAtattc                                                                                                                                     >  1:9379/1‑141 (MQ=255)     gCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTa                                                                                                                                 <  2:221991/141‑1 (MQ=255)                           cTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACttgtt                                                                                                           <  1:56362/141‑1 (MQ=255)                                  atatGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATaat                                                                                                    >  2:89291/1‑139 (MQ=255)                                            aaCCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAAt                                                                                          >  2:93864/1‑141 (MQ=255)                                                accCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGAt                                                                                      >  2:47721/1‑141 (MQ=255)                                                 cccAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAGTTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCTTACCTTCAGTTAATATAtt                                                                                     >  2:319173/1‑141 (MQ=255)                                                     aTTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGata                                                                                 >  1:290520/1‑141 (MQ=255)                                                        tGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATAttt                                                                              >  1:130510/1‑141 (MQ=255)                                                        tGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATAttt                                                                              <  1:203140/141‑1 (MQ=255)                                                                                    tGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCaaa                                                  <  2:54552/141‑1 (MQ=255)                                                                                      aGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACt                                                <  1:108019/141‑1 (MQ=255)                                                                                          cTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGt                                            >  1:198852/1‑141 (MQ=255)                                                                                              aGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGttttt                                        <  1:66093/141‑1 (MQ=255)                                                                                                           ccAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTtaat                           <  2:417524/141‑1 (MQ=255)                                                                                                                 aaGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGa                     >  2:96282/1‑141 (MQ=255)                                                                                                                            aGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAat          >  1:310428/1‑141 (MQ=255)                                                                                                                               tCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATAtt       <  2:132783/141‑1 (MQ=255)                                                                                                                                  ataAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGcc    <  1:369991/141‑1 (MQ=255)                                                                                                                                    aaaTATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCat  >  1:375819/1‑141 (MQ=255)                                                                                                                                          |                                                                                                                                       AATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCAT  >  CP000730/39241‑39512

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑