Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 39,195 T→C 62.2% K647E (AAA→GAA)  mecA ← penicillin binding protein 2 prime

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073039,1950TC62.2% 19.3 / 28.0 32K647E (AAA→GAA) mecApenicillin binding protein 2 prime
Reads supporting (aligned to +/- strand):  ref base T (9/3);  new base C (8/12);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 7.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.89e-01

AGCACTTATTAATAATTCATAATAAAACAGTGAAGCAACCATCGTTACGGATTGCTTCACTGTTTTGTTATTCATCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTG  >  CP000730/39060‑39332
                                                                                                                                       |                                                                                                                                         
aGCAATTATTAATTATTCATAATAAAACAGTGAAGCAACCATCGTTACGGATTGCTTCACTGTTTTGTTATTCATCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTTGGCAt                                                                                                                                      <  2:284965/141‑1 (MQ=255)
     aaaataagaaTTCATAATAAATCAGTGAAGAAACCATCGTTACGGATTGATTCACTGTTTTGTTATTCATCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTTGGCATTGtag                                                                                                                                 <  1:93683/137‑1 (MQ=255)
                  aaaataaAACAGTGAAGCATCCATCGTTACGGATTACATCAATGTTTTGTTATTCATCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGATTTGTCTTCATTGTAGCTAGCCATTCCtt                                                                                                                    <  1:94399/139‑1 (MQ=255)
                       aaaaCAGGGAAGCAACCATCGTTACGGATTACTTCACTGTTTTGTTATTCATCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTCGGCATTGTAGCTAGCCATTCCTTTATCt                                                                                                               <  1:155145/141‑1 (MQ=255)
                            aGTGAAGCAACCATCGTTACGGATTGCTTCACTGTTTTGTTATTCATCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTCGGCATTGTAGCTAGCCATTCCTTTATCTTGTAc                                                                                                          >  1:8388/1‑141 (MQ=255)
                                   cAACCATCGTTACGGATTGCTTCACTGTTTTGTTATTCATCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTa                                                                                                   <  2:74352/141‑1 (MQ=255)
                                           gtaaCGGATTGCTTCACTGTTTTGTTATTCACCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTCGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAAt                                                                                           <  1:133355/138‑1 (MQ=255)
                                           gTTACGGATTGCTTCACTGTTTTGTTATTCATCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGCGATTTCGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAAt                                                                                           >  1:275704/1‑141 (MQ=255)
                                           gTTACGGATTGCTTCACTGTTTTGCTATTCATCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTGTACCTGAGATTTCGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAAt                                                                                           >  1:288265/1‑141 (MQ=255)
                                                    tGCTCCACTTTTTTGTTATTTATCTAAATCGTATTTTTTATTACCCTGCTCATATAGCTCATCATACACTTTACCTGAGATTTCGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCcatcat                                                                                  <  2:95746/141‑1 (MQ=255)
                                                           ctattttGTTATTCATCTATATCGTATTTTTTATTACCGTTCTAATATAGCTCATCATACACTTTACCTGAGATTTCGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGttt                                                                           <  1:84529/138‑1 (MQ=255)
                                                               tttGTTATTCATCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTCGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATATTTAACATTAATAGCCATCATCATGTTTGGAt                                                                       >  1:51890/1‑141 (MQ=255)
                                                                ttGTTATTCATCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGAtt                                                                      >  1:142979/1‑141 (MQ=255)
                                                                       tCATCTATATCGTATTTTTTATTACGGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTCGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCTTTATCTTGTTTCGATTATCtttt                                                               >  1:371477/1‑140 (MQ=255)
                                                                                  ggaTTTTTTAGTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTCGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGAATTATCTTTATCATATGatat                                                    <  2:143104/139‑1 (MQ=255)
                                                                                     ttttttATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATaaa                                                 >  1:142139/1‑141 (MQ=255)
                                                                                         ttattaCCGTTCTCATATACCGCATCATACACTTTACCTGAGATTTCGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAaccac                                             <  2:13775/141‑1 (MQ=255)
                                                                                                       atatAGCTCATCATACACTTTACCTGAGATTTTGGCATTGTAGCTAGACATTCCTTTATCTTGTACATCTTTAACATTAATAGCCGTCATCATGTTTGGATTATCTTTATCGTATCATATAAACCACCCACTTTGTCTGcc                               >  2:386876/1‑141 (MQ=255)
                                                                                                         ataGCTCATGCTACACTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAg                             >  1:290070/1‑141 (MQ=255)
                                                                                                           aGCTCATCATACACTTTACCTGAGATTTCGGCATTGTAGCTAGCCATTCCTTCATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTACCAGtt                           <  2:248697/141‑1 (MQ=255)
                                                                                                                atcatACACTTTACCTGAGATTTCGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTcc                      >  1:204308/1‑141 (MQ=255)
                                                                                                                   ttACACTTTACCTGGGATTTCGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTg                   <  2:1055/140‑1 (MQ=255)
                                                                                                                    tgcacTTTACCTGAGGTTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAGCATTAATAGCCATCATCATGTTTGGCTTATCATTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGt                  >  1:186975/3‑141 (MQ=255)
                                                                                                                     acacTTTACCTGAGATTTTGGCCTTGTAACTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTGGGATTATCTTTATCATATGATATAAACCACGCAATTTGTCTGCCAGTTTCTGCTTGtt                 >  1:29482/1‑141 (MQ=255)
                                                                                                                     acacTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTGCTTTATCTTGTATATATTTCCAATTCATAGCCCTCGTCATGTTTGGATTATCTTTAGCATATTATATAAACCACCCAATTTGTTCGCCAGTTTCTACTTGtt                 >  2:229832/1‑141 (MQ=255)
                                                                                                                     acacTTTACCTGAGATTTCGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATGGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTTTCTGCCAGTTTCTCCTTGtt                 >  1:170886/1‑141 (MQ=255)
                                                                                                                      accgtcgccTGAGATTTCGGCATTGTAGCTAGCCATTCCTTTATTATGTACATCTTTAACATTAATAGCCATCATCATGTTTGGGTTATCTTTATCATATGATAGAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGgtt                <  2:4330/134‑1 (MQ=255)
                                                                                                                       acTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTTTCATATGATATAAACCACCCAATTTGTCTGCCGGTTTCTCCTTGtttt               >  2:270932/1‑140 (MQ=255)
                                                                                                                       acTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTCCCACCTTCTCCTTCTTTc               >  2:333878/1‑141 (MQ=255)
                                                                                                                           tACCGGCGATTTCGGCCTTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTACCATTAATAGCCATCATCATGTTTGGATTAACTCTATAATACGATCTCAACCACCCAATTAGTCTACCATTTTCTCCTTGTTTCtttt           <  2:331242/141‑1 (MQ=255)
                                                                                                                                gagaTTTCGGCATTGTAGCTAGCCTTTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCTTTTTgtgt      >  1:192116/1‑141 (MQ=255)
                                                                                                                                  gaTTTCGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTc    <  1:29841/141‑1 (MQ=255)
                                                                                                                                    tttcGGTATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTTTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTg  <  2:345967/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                         
AGCACTTATTAATAATTCATAATAAAACAGTGAAGCAACCATCGTTACGGATTGCTTCACTGTTTTGTTATTCATCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTTATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTG  >  CP000730/39060‑39332

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: