| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | CP000730 | 39,378 | C→A | 100% | 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 | ||
| * | CP000730 | 39,378 | 0 | C | A | 84.6% | 25.0 / 0.0 | 13 | D586Y (GAT→TAT) | mecA | penicillin binding protein 2 prime | 
| Reads supporting (aligned to +/- strand): ref base C (0/0); major base A (5/6); minor base T (1/1); total (6/7) | |||||||||||
| 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.39e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
GTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCC  >  CP000730/39256‑39510                                                                                                                          |                                                                                                                                    gTTTGGATTGTCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCACTTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTtttatt                                                                                                                    >  2:10421/1‑141 (MQ=255)                    atGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCAt                                                                                                >  1:156010/1‑141 (MQ=255)                                                   ttCTACTTGTTTCATTTTGAGTTCTGCAAAGCCGGGTTTGCCACTTACGTGTGCATAAGAGCTATAAAAAGATGCTTTATGTGTTTTATTTGCGACGTGTTGCATACATTCAGTTACTAGAGTGATGTTGTCTTTGGAgat                                                                 >  2:152024/1‑141 (MQ=255)                                                          tGTTTCTTTTCGAGTGCTGCCGTGCCGGATGTGCCACTTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATAttt                                                          <  2:57052/141‑1 (MQ=255)                                                              tCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTttct                                                      <  2:44552/141‑1 (MQ=255)                                                                         tCTGCAGTACCGTATTTGCCAATTAAGTTTGCATAAGAACTATAAATTTATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTg                                           <  1:244297/141‑1 (MQ=255)                                                                             cAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCAAAAGTTCGtctt                                       >  2:108915/1‑141 (MQ=255)                                                                               gTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTc                                     <  1:115421/141‑1 (MQ=255)                                                                                   cGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCgtgt                                 >  1:49547/1‑141 (MQ=255)                                                                                                               tCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATAtt     <  2:160328/141‑1 (MQ=255)                                                                                                                 tataAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGc   <  2:5211/141‑1 (MQ=255)                                                                                                                  ataAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGcc  >  2:106392/1‑141 (MQ=255)                                                                                                                  aaaaaTATATTCCTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATAATTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAAGTGCGTTAATATTGcc  <  1:95523/139‑1 (MQ=255)                                                                                                                          |                                                                                                                                    GTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCC  >  CP000730/39256‑39510 | 
| 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: ‑ |