| 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 | ||
| * | CP000730 | 39,246 | 0 | C | T | 50.0% | ‑5.8 / 20.6 | 18 | A630T (GCT→ACT) | mecA | penicillin 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 | |||||||||||
GATTGCTTCACTGTTTTGTTATTCATCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTT‑ATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTA  >  CP000730/39109‑39385                                                                                                                                          |                                                                                                                                           gATTGCTTCACTGTTTTGTTATTCATCTATATCGTATTATTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTT‑ATCTTGTACATCTTTAACATTAATAGCcat                                                                                                                                          >  1:221019/1‑141 (MQ=255)                                cGTAGTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTT‑ATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGata                                                                                                          <  2:264900/141‑1 (MQ=255)                                              cggTTCTCAGATATCTTATCATACACTTTACCTGAGATCTTGGCATTGTAGCTAGCCATTCCTTT‑ATCTTGTACATCTTTAACATTACTAGTCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAAtt                                                                                            <  2:16905/139‑1 (MQ=255)                                                    tcATATAGCTCATCCTACCCTTTACCTGAGATTTTGGCTTCGCAGGCTTCTTGTTCTTT‑ATCGTGTACATCATTAACATTAATAGTCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTg                                                                                      <  1:161777/141‑1 (MQ=255)                                                                tcatACACTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTT‑ATCTTGTACATCTTTAACATTAATAGTCATCATCATGTTTGGACTATCTTTATCGTATGATATAAACCACCCAATTTGTCGACCAGTTtctcat                                                                          >  2:236429/1‑139 (MQ=255)                                                                  ctACACTTTTCCTGAGATTTTGGCATTGTAGCTAGCCATAACTTTAATCCACAACCT‑TTAATGACTAATAGTCATCATTATGTTTCTAGTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTg                                                                        <  2:332277/140‑1 (MQ=255)                                                                    acacTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTT‑ATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGtt                                                                      <  1:189813/141‑1 (MQ=255)                                                                    acacTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTT‑ATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGtt                                                                      <  1:255979/141‑1 (MQ=255)                                                                               gagaTTTTGGCATTGTAGCTAGCCATTCCTTT‑ATCTTGTACATCTTTAACATTAATAGTCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGt                                                           >  2:117440/1‑141 (MQ=255)                                                                                                    gcCATTCCTTT‑ATCTTGTACATCTTTAACATTAATAGTCATCATCATGTTTGGAGTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCa                                      >  1:229218/1‑141 (MQ=255)                                                                                                      cATTCCTTT‑ATCTTGTACATCTTTAACATTAATAGTCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAAt                                    <  1:117440/141‑1 (MQ=255)                                                                                                      cATTCCTTT‑ATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAAt                                    <  2:26007/141‑1 (MQ=255)                                                                                                        ttCTTTT‑ACCTTGTACATCTTCAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTa                                  <  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)                                                                                                                                          |                                                                                                                                           GATTGCTTCACTGTTTTGTTATTCATCTATATCGTATTTTTTATTACCGTTCTCATATAGCTCATCATACACTTTACCTGAGATTTTGGCATTGTAGCTAGCCATTCCTTT‑ATCTTGTACATCTTTAACATTAATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTA  >  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: ‑ |