Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 39,378 C→A 100% 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
*CP00073039,3780CA100.0% 71.3 / ‑6.5 32D586Y (GAT→TAT) mecApenicillin binding protein 2 prime
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (18/13);  minor base T (0/1);  total (18/14)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCATT  >  CP000730/39247‑39513
                                                                                                                                   |                                                                                                                                       
catcatcatGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATATATATATTCTTTAtt                                                                                                                                >  2:142929/1‑140 (MQ=255)
catcatcatGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGACCTATAAATATATTCTTTAtt                                                                                                                                >  2:38203/1‑140 (MQ=255)
   catcatGTTTGGATAATCTTTCTCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGCGTTCTGCAGTACCTGCTTTGCCAATTAAGTTTGCATAAGATCTATACATATATTCTTCATgtgt                                                                                                                             >  2:184571/1‑141 (MQ=255)
    atgatGTTTGGATTATCTTTATCATATGATATAAACAACCCAATTTGTCTGCAAGTTTGTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGtt                                                                                                                            >  1:180893/4‑141 (MQ=255)
    atcatGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGtt                                                                                                                            >  2:67064/1‑141 (MQ=255)
      catGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAAAATATTCTTTATGTGgttt                                                                                                                          >  2:92846/1‑141 (MQ=255)
       atGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGttttt                                                                                                                         >  2:77346/1‑140 (MQ=255)
             ggATTATTTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACg                                                                                                                   <  2:120601/141‑1 (MQ=255)
             ggATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACg                                                                                                                   <  2:101948/141‑1 (MQ=255)
                      ttATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGc                                                                                                          <  1:79109/141‑1 (MQ=255)
                       tATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCa                                                                                                         >  1:59199/1‑141 (MQ=255)
                           atatGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGTAGTACGGGATTTGCCAATTAAGTTTGCATGAGATCTATAAATGTATTCTTTCTGTGTTTTATTTAGGACTTGTTGCATAcc                                                                                                     >  1:92100/1‑141 (MQ=255)
                                atatAATCCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAg                                                                                                <  2:71006/141‑1 (MQ=255)
                                 tttaAACCACCCAATATGTCTTCCAGTTGCTCCTTGGTTCATTTTGAGTTCTGCAGTACCGGAGTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGGTTGATTTACGACTTGTGGCATACCATCAGt                                                                                               <  1:53220/139‑1 (MQ=255)
                                     aaCCACACAATTTGTCTGCCAGTTACACCTTGTTTCATGTTTAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAAt                                                                                           >  2:64438/1‑141 (MQ=255)
                                      accaccCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGAACTAAAGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATa                                                                                          >  1:50878/1‑141 (MQ=255)
                                      acaaccCAATTTGTGTGCCAGTTTCTCCTTGTTTCATCTTGAGGTCTGCAGTCCCTGATTTGCCAATTAACGTTGCCTAAGGTCTATAAATATATTCTTTATGTGATTTATTTACGACTTGTTACCTAGCATCAGttgaga                                                                                          >  2:118438/4‑137 (MQ=255)
                                                 tGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATAttt                                                                               <  1:47605/141‑1 (MQ=255)
                                                       ttaGTTTCGCCTTTTTTGATTGTGAGTTGTGCAGTACCGGATTTGCTAATTAAGTTTTCATAAGTTTTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTg                                                                         <  2:263523/139‑1 (MQ=255)
                                                        cAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTgg                                                                        >  2:86616/1‑141 (MQ=255)
                                                                         aTTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATTATATTTttcttc                                                       >  2:91944/1‑141 (MQ=255)
                                                                          ttttGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTcc                                                      >  1:133415/1‑141 (MQ=255)
                                                                              gAGTTCCGCAGTACCGTATTTGCCAATTAAGGTTGCATAAGATCTATAAAGATATTCTTTATGAGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGCAATAATCTTTTTCTTCCAAAc                                                  <  2:185363/141‑1 (MQ=255)
                                                                                gTTATGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAGACtt                                                >  2:31699/1‑141 (MQ=255)
                                                                                   cTGCAGTACCGGATTTTCCTATTAAGTTTGCATAAGATCTATAAATTTATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGt                                             <  2:216062/141‑1 (MQ=255)
                                                                                      cAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGtttt                                          >  1:158970/1‑141 (MQ=255)
                                                                                                 tttccAATTAATTTTGCATAATGTCTATTAATATATTCTTTATGTTTTTTATTTACGACTTTTTTCATATCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCtttt                               <  2:180893/138‑1 (MQ=255)
                                                                                                 ttGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGGTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCtttt                               <  1:67064/141‑1 (MQ=255)
                                                                                                              ttGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTg                  <  2:158970/141‑1 (MQ=255)
                                                                                                                           ataAATATATTCTGTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGcc     <  1:38203/141‑1 (MQ=255)
                                                                                                                            taAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCa    <  1:91466/141‑1 (MQ=255)
                                                                                                                              aaTATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCACACTTTGTTTTTCGTGTCTATTAATAAGTTGGGTGCGTTAATATTGCCatt  >  2:119071/1‑141 (MQ=255)
                                                                                                                                   |                                                                                                                                       
CATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCATT  >  CP000730/39247‑39513

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