Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	39,378	C→T	D586N (GAT→AAT)	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	T	100.0%	113.9 / NA	35	D586N (GAT→AAT)	mecA	penicillin binding protein 2 prime
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/20);  total (15/20)

TTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCATTAT  >  CP000730/39239‑39515                                                                                                                                            |                                                                                                                                          ttAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTAt                                                                                                                                   <  1:81876/148‑1 (MQ=255)     tAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATgtgt                                                                                                                               <  2:288359/148‑1 (MQ=255)     tAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATgtgt                                                                                                                               <  1:135707/148‑1 (MQ=255)     tAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATgtgt                                                                                                                               <  1:350555/148‑1 (MQ=255)     tAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATgtgt                                                                                                                               <  2:411608/148‑1 (MQ=255)            catcatGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTtttattt                                                                                                                        <  2:25282/148‑1 (MQ=255)            catcatGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTtttattt                                                                                                                        <  2:344842/148‑1 (MQ=255)            catcatGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTtttattt                                                                                                                        <  2:214211/148‑1 (MQ=255)                                  tCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAg                                                                                                  >  2:226766/1‑148 (MQ=255)                                  tCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTCGAGTACTGCAGTACCGGATTTGCCAATTAAGTTTGCATACGATCTATAAATACTTTCTTTATGTGTTTTATTTACGCCTTGTTGCCTACCATCAg                                                                                                  >  2:23309/1‑148 (MQ=255)                                         atatAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACTATCAGTTAATat                                                                                           >  2:96013/1‑147 (MQ=255)                                         atatAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAg                                                                                           >  2:93014/1‑148 (MQ=255)                                              aaCCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGa                                                                                      >  2:273955/1‑148 (MQ=255)                                                ccaccCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGata                                                                                    >  2:273502/1‑148 (MQ=255)                                                ccaccCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGata                                                                                    >  1:141518/1‑148 (MQ=255)                                                       aTTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCtt                                                                             <  2:96067/148‑1 (MQ=255)                                                        tttGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCttt                                                                            >  2:388402/1‑148 (MQ=255)                                                              tGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAaat                                                                      >  1:109877/1‑148 (MQ=255)                                                              tGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAaat                                                                      >  2:272319/1‑148 (MQ=255)                                                              tGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAaat                                                                      >  1:293620/1‑148 (MQ=255)                                                                  aGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAataata                                                                  >  2:180760/1‑148 (MQ=255)                                                                                  aTTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATACGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATAAGTTAATAGATTGATATTTTCTTTGGATATAATGTTTTTCTTCCATACtt                                                  >  2:166075/1‑148 (MQ=255)                                                                                   ttttGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACttt                                                 <  2:260239/148‑1 (MQ=255)                                                                                   ttttGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACttt                                                 <  2:73032/148‑1 (MQ=255)                                                                                   ttttGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACttt                                                 <  2:72963/148‑1 (MQ=255)                                                                                                aGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCt                                    <  1:52131/148‑1 (MQ=255)                                                                                                aGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCt                                    <  1:225730/148‑1 (MQ=255)                                                                                                 gTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCtt                                   <  2:38942/148‑1 (MQ=255)                                                                                                      ggATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTtaat                              <  2:281194/148‑1 (MQ=255)                                                                                                           tGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTg                         <  2:49675/148‑1 (MQ=255)                                                                                                           tGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTg                         <  2:223206/148‑1 (MQ=255)                                                                                                                       ttGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAat             >  1:21854/1‑148 (MQ=255)                                                                                                                          cATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATAtt          <  2:64103/148‑1 (MQ=255)                                                                                                                              aGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCa      <  1:273502/148‑1 (MQ=255)                                                                                                                                  cTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCattat  >  1:326496/1‑148 (MQ=255)                                                                                                                                            |                                                                                                                                          TTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCATTAT  >  CP000730/39239‑39515

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