Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 39,378 C→A 50.2% 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,3780CA50.2% ‑3.4 / 58.3 36D586Y (GAT→TAT) mecApenicillin binding protein 2 prime
Reads supporting (aligned to +/- strand):  ref base C (7/11);  new base A (8/10);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

ATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCA  >  CP000730/39242‑39511
                                                                                                                                        |                                                                                                                                     
attgCCATCATCATGTTTGGATTATCTTTATTATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATAtattct                                                                                                                                   <  1:375913/138‑1 (MQ=255)
 tagtcatcatcatGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTATGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATcttctt                                                                                                                                  <  2:90408/137‑1 (MQ=255)
    tcatcatcatGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTAt                                                                                                                               <  1:132705/140‑1 (MQ=255)
       tcatcatGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTAtgtg                                                                                                                            <  2:336893/141‑1 (MQ=255)
        catcatGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTAGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTACCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATgtgt                                                                                                                           <  1:132065/141‑1 (MQ=255)
                      tATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACtt                                                                                                             <  1:163731/141‑1 (MQ=255)
                       aTCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACttg                                                                                                            >  1:268388/1‑141 (MQ=255)
                       aTCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACttg                                                                                                            <  1:107659/141‑1 (MQ=255)
                       aTCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACttg                                                                                                            <  1:2606/141‑1 (MQ=255)
                          tttATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACttgttg                                                                                                         >  1:339613/1‑141 (MQ=255)
                              tCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATa                                                                                                     >  2:10948/1‑141 (MQ=255)
                                  atGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCAt                                                                                                 <  1:218643/141‑1 (MQ=255)
                                         aaaCCACCCACTTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTCACGGCTTTGGCCATTCAGTTTGCATAAGCTCTATAAATATCTTCTTTATGTGTTTTCTTTACGACTTGTTGCATACCATCAGTTaa                                                                                          <  1:189069/141‑1 (MQ=255)
                                            ccaccCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAg                                                                                       >  1:268435/1‑141 (MQ=255)
                                              accCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGAt                                                                                     >  1:331525/1‑141 (MQ=255)
                                                   aTTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGata                                                                                >  2:150404/1‑141 (MQ=255)
                                                   aTTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGata                                                                                >  2:330438/1‑141 (MQ=255)
                                                       gTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATAtttt                                                                            <  2:54161/141‑1 (MQ=255)
                                                          tGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCtt                                                                         >  2:285421/1‑141 (MQ=255)
                                                           gCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCttt                                                                        <  1:272641/141‑1 (MQ=255)
                                                                            tCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTttct                                                       >  2:367557/1‑141 (MQ=255)
                                                                                tttGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCa                                                   >  2:101961/1‑141 (MQ=255)
                                                                                   gAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGACATAATATTTTTCTTCCAAAc                                                <  2:244726/141‑1 (MQ=255)
                                                                                    aGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACt                                               <  2:118274/141‑1 (MQ=255)
                                                                                        cTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGt                                           <  1:10956/141‑1 (MQ=255)
                                                                                        cTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGt                                           <  1:10948/141‑1 (MQ=255)
                                                                                           cAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGtttt                                        <  1:189776/141‑1 (MQ=255)
                                                                                           cAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGtttt                                        <  1:201651/141‑1 (MQ=255)
                                                                                              tACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCg                                     <  1:348603/141‑1 (MQ=255)
                                                                                              tACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCg                                     >  1:37690/1‑141 (MQ=255)
                                                                                                ccGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCgtg                                   <  2:207368/141‑1 (MQ=255)
                                                                                                  ggATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTc                                 <  2:143810/141‑1 (MQ=255)
                                                                                                     tttGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGGAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCttt                              >  2:8783/1‑141 (MQ=255)
                                                                                                                    tGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTTTTTTATACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGc               >  2:91179/1‑141 (MQ=255)
                                                                                                                     gCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCg              >  1:68869/1‑141 (MQ=255)
                                                                                                                          aGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAat         <  1:150404/141‑1 (MQ=255)
                                                                                                                          aGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAat         <  1:330438/141‑1 (MQ=255)
                                                                                                                                 taAATATATTCTTTATGTGTTTTATTTACGACTTGATGCATAGCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCa  >  1:324285/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                     
ATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCA  >  CP000730/39242‑39511

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