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% 74.6 / NA 25D586Y (GAT→TAT) mecApenicillin binding protein 2 prime
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/9);  total (16/9)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCAT  >  CP000730/39249‑39512
                                                                                                                                 |                                                                                                                                      
acatcatGTTTGGATTATCTTTATCATATGAAAGTATCCTCCCAAATTAACTTCCATTATCCCCTTATTTCAGGGTGACTTCTTCATTACCTGATTTTCCAATTAAGTTTTCATAAGATCTATAAATATATTCTTTAtgtg                                                                                                                             <  2:95798/140‑1 (MQ=255)
   tcatGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTATGGTATATTCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGttt                                                                                                                          >  1:79652/1‑141 (MQ=255)
               tATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACtt                                                                                                              >  1:47807/1‑141 (MQ=255)
                               tataAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGt                                                                                              <  2:22926/141‑1 (MQ=255)
                                       accCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGCTTTGCCAATTAAGTTTGCATACGATCTATCACCACATTCTCTATGTCTTTAAGTTGCGACTTGTAGCATACCATCACTTAACAGAt                                                                                      >  2:270120/1‑141 (MQ=255)
                                                     ccAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTg                                                                        >  1:257108/1‑141 (MQ=255)
                                                            ctcCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAaataa                                                                 >  2:200546/1‑141 (MQ=255)
                                                              ccTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGCTCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCCGTTAAGAGAATGATATCTTCTTTGGAaataatt                                                               >  2:153949/1‑140 (MQ=255)
                                                                   tttCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATAttttt                                                          >  2:206644/1‑141 (MQ=255)
                                                                             ggTTCTGCAGGTCCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACt                                                <  1:28539/140‑1 (MQ=255)
                                                                                     aGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAACATATTCTTTATGTGTTCTATTTACGACTTGTTGCATATCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGttttt                                        <  2:141847/141‑1 (MQ=255)
                                                                                       tACGGGATTTGGGAATTAAGTTTGCATAAGATCTATAAATATATTATATAAGAGATTTATATACGACTAGTTGCATACTATCAGTTAATAGTTTGATATTTTCTTTGGAAATAATATTTTTATTACAAACTTTGttttttg                                      >  2:56677/1‑139 (MQ=255)
                                                                                        cccGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCgt                                     <  2:269496/140‑1 (MQ=255)
                                                                                             aTTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCtt                                >  2:112575/1‑141 (MQ=255)
                                                                                                    aaTTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCATTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTtaataa                         <  2:126450/141‑1 (MQ=255)
                                                                                                        aaGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATGAGTGa                     >  1:266084/1‑141 (MQ=255)
                                                                                                         aGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAg                    <  2:91493/141‑1 (MQ=255)
                                                                                                          gTTTGCATAAGATCTATAAATATATTCTTTGTGTGTTTTATTTACGACTTGTTGCATACCATCAGTTTATAGTTTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAgt                   >  1:260977/1‑140 (MQ=255)
                                                                                                              gCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCg               >  2:1582/1‑141 (MQ=255)
                                                                                                              gCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTAGTGGCTTTTAATAAGTGAGGAGCg               >  2:177382/1‑141 (MQ=255)
                                                                                                               cATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGt              <  2:257108/141‑1 (MQ=255)
                                                                                                                aTAAGATCTATAAATAAATTCCTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTCGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGtt             >  1:147246/1‑141 (MQ=255)
                                                                                                                  aaGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTTAGGTGCGTTat           >  2:224501/1‑140 (MQ=255)
                                                                                                                           aaaTATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTTCCat  <  1:194594/141‑1 (MQ=255)
                                                                                                                           aaaTATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCat  >  1:95538/1‑141 (MQ=255)
                                                                                                                                 |                                                                                                                                      
TCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCAT  >  CP000730/39249‑39512

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