Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 39,378 C→T 100% 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
*CP00073039,3780CT74.1% 47.0 / 16.5 31D586N (GAT→AAT) mecApenicillin binding protein 2 prime
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (9/13);  minor base A (7/1);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 3.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCATT  >  CP000730/39241‑39513
                                                                                                                                         |                                                                                                                                       
aaTAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATAttttc                                                                                                                                      >  2:9092/1‑141 (MQ=255)
 aTAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAATCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATAttttct                                                                                                                                     >  2:258842/1‑141 (MQ=255)
            catGTTTGGATTATCTTTATCATATGATATAAACCACAAAATTTGTCTGCCAGGTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGtttt                                                                                                                          >  2:210599/1‑141 (MQ=255)
              tGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCGTAAGATCTATAAATATTTTCTTTATGTGTtttat                                                                                                                        >  2:167520/1‑141 (MQ=255)
              tGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTtttat                                                                                                                        >  2:55175/1‑141 (MQ=255)
               gTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTtttatt                                                                                                                       <  2:236066/141‑1 (MQ=255)
                tttGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTtttattt                                                                                                                      <  1:160421/141‑1 (MQ=255)
                              aTCATATGATATAAATCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCAt                                                                                                        <  1:258842/141‑1 (MQ=255)
                                    tGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATc                                                                                                  >  1:209402/1‑141 (MQ=255)
                                              caccCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGa                                                                                        >  2:139817/1‑141 (MQ=255)
                                                                  ttCTCCTTGTTTCGTTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAaat                                                                    <  2:182785/141‑1 (MQ=255)
                                                                  ttCTCCTTGTTTCATTTTGAGTTCTGCACTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAaat                                                                    >  1:98684/1‑141 (MQ=255)
                                                                    ctcCATGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAaataa                                                                  >  1:191271/1‑141 (MQ=255)
                                                                              cATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTttctt                                                        <  2:236687/141‑1 (MQ=255)
                                                                                     aGTTCTTCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATGTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACt                                                 <  1:201240/141‑1 (MQ=255)
                                                                                     aGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACt                                                 <  1:9092/141‑1 (MQ=255)
                                                                                            cAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGtttt                                          >  1:117736/1‑141 (MQ=255)
                                                                                             aGTATCGGATTTTCTAATTAAGTTTGTATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGttttt                                         <  2:243286/141‑1 (MQ=255)
                                                                                                  cGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGGTTTATTAACGACTTGTTGCATACCATCAGTTAATAGCTTGATATTTTGTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCgcgt                                    >  2:121976/1‑141 (MQ=255)
                                                                                                    gATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATCATATTTTTCTACCAAACTATGtgtaccgattct                                  >  2:122682/1‑130 (MQ=255)
                                                                                                    gATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAAAAATATTTTTCTTCCAAACTTTGTTTTTTGTGTCt                                  >  1:160396/1‑141 (MQ=255)
                                                                                                     aTTTGCCAATTAAGTTTGCATAAGAGCTATAAATATATTCTTTAGGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCtt                                 >  1:153601/1‑141 (MQ=255)
                                                                                                      tttGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGATAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTAGTTTCttt                                >  2:13083/1‑141 (MQ=255)
                                                                                                       ttGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCtttt                               <  1:96872/141‑1 (MQ=255)
                                                                                                          ccAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTAAGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTAGTTTTTCGGGTCTTTtaat                            >  1:78180/1‑141 (MQ=255)
                                                                                                                   tttGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGt                   >  1:275644/1‑141 (MQ=255)
                                                                                                                      gCATAAGATCTATAAATATTTTCTTTATTTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCg                <  2:102171/141‑1 (MQ=255)
                                                                                                                        aTAAGCTCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGtt              <  2:176135/141‑1 (MQ=255)
                                                                                                                        aTAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGtt              <  2:176121/141‑1 (MQ=255)
                                                                                                                        aTAAAATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGtt              <  2:265136/141‑1 (MQ=255)
                                                                                                                          aagtttaataAATTTTTTCCTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTaa            <  2:78934/134‑1 (MQ=255)
                                                                                                                             aTCTATAAATATTTTCTTTATGTGTGTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTACAAATCGACTTTTTTCGTTGCTTTTAATAAGTGAGGTTCGGTAatat         >  2:268257/1‑141 (MQ=255)
                                                                                                                                    aaTATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCatt  <  1:67747/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
AATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCATT  >  CP000730/39241‑39513

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