Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 39,246 C→T 50.0% A630T (GCT→ACT)  mecA ← penicillin binding protein 2 prime

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073039,2460CT50.0% ‑5.8 / 20.6 18A630T (GCT→ACT) mecApenicillin binding protein 2 prime
Reads supporting (aligned to +/- strand):  ref base C (2/7);  new base T (5/4);  total (7/11)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GATTGCTTCACTGTTTTGTTATTCATCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTA  >  CP000730/39109‑39385
                                                                                                                                          |                                                                                                                                           
gATTGCTTCACTGTTTTGTTATTCATCTATATCGTATTATTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCcat                                                                                                                                          >  1:221019/1‑141 (MQ=255)
                                cGTAGTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGata                                                                                                          <  2:264900/141‑1 (MQ=255)
                                              cggTTCTCAGATATCTTATCATACACTTTACCTGAGATCTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTACTAGTCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAAtt                                                                                            <  2:16905/139‑1 (MQ=255)
                                                    tcATATAGCTCATCCTACCCTTTACCTGAGATTTTGGCTTCGCAGGCTTCTTGTTCTTTATCGTGTACATCATTAACATTAATAGTCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTg                                                                                      <  1:161777/141‑1 (MQ=255)
                                                                tcatACACTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGTCATCATCATGTTTGGACTATCTTTATCGTATGATATAAACCACCCAATTTGTCGACCAGTTtctcat                                                                          >  2:236429/1‑139 (MQ=255)
                                                                  ctACACTTTTCCTGAGATTTTGGCATTGTAGCTAGCCATAACTTTAATCCACAACCTTTAATGACTAATAGTCATCATTATGTTTCTAGTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTg                                                                        <  2:332277/140‑1 (MQ=255)
                                                                    acacTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGtt                                                                      <  1:189813/141‑1 (MQ=255)
                                                                    acacTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGtt                                                                      <  1:255979/141‑1 (MQ=255)
                                                                               gagaTTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGTCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGt                                                           >  2:117440/1‑141 (MQ=255)
                                                                                                    gcCATTCCTTTATCTTGTACATCTTTAACATTAATAGTCATCATCATGTTTGGAGTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCa                                      >  1:229218/1‑141 (MQ=255)
                                                                                                      cATTCCTTTATCTTGTACATCTTTAACATTAATAGTCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAAt                                    <  1:117440/141‑1 (MQ=255)
                                                                                                      cATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAAt                                    <  2:26007/141‑1 (MQ=255)
                                                                                                        ttCTTTTACCTTGTACATCTTCAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTa                                  <  1:91423/141‑1 (MQ=255)
                                                                                                                 tCTTGTACATCTTTAACATTAATAGTCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCa                          >  2:293478/1‑141 (MQ=255)
                                                                                                                 tCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCa                          <  2:339303/141‑1 (MQ=255)
                                                                                                                                     aaTAGTCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATAtcttc      >  2:283892/1‑141 (MQ=255)
                                                                                                                                      aTAGCGATCATCATGTTTGGATTATCCTTATCATATGATATAAACCTCCCCATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCACTAACGGATTTTCCAATTAAGTTTGCATAAGATCTATAAATAtattct     <  2:226490/141‑1 (MQ=255)
                                                                                                                                         gCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTa  >  1:104871/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                           
GATTGCTTCACTGTTTTGTTATTCATCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTA  >  CP000730/39109‑39385

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: