| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | USA300TCH1516_ALE | 1,797,772 | G→C | 100% | R7G (CGC→GGC) | phoR ← | Alkaline phosphatase synthesis sensor protein PhoR |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | USA300TCH1516_ALE | 1,797,772 | 0 | G | C | 95.5% | 59.3 / ‑3.7 | 22 | R7G (CGC→GGC) | phoR | Alkaline phosphatase synthesis sensor protein PhoR |
| Reads supporting (aligned to +/- strand): ref base G (0/0); major base C (11/10); minor base T (1/0); total (12/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
ATAAATATATGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACA > USA300TCH1516_ALE/1797635‑1797898 | ataaataTATGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCcgt > 1:5299/1‑140 (MQ=255) cATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAAGAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCAtt < 2:266831/140‑1 (MQ=255) cATCTGTTGCTAGAGCTTGTTCGTGTTGTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCAtt < 1:48764/140‑1 (MQ=255) gCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCtt < 2:332518/140‑1 (MQ=255) gTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTATAGGTGTTTccg > 2:277564/1‑139 (MQ=255) tCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATAATTGTTCTTTAGGTCTTTCCAAt > 1:280881/1‑140 (MQ=255) cGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAAtt > 2:184951/1‑140 (MQ=255) cGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGGGGTGAAACTTCATCATTGTTCTTTCGGTCTTTCCAAtt > 2:18069/1‑140 (MQ=255) gtTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCTGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTtat > 1:310788/1‑140 (MQ=255) atatGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGcac < 2:186801/140‑1 (MQ=255) tgtgATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAAGAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGt < 1:351139/140‑1 (MQ=255) aTGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTta < 2:291024/140‑1 (MQ=255) ggAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGACTTTCCAATTTATAGCCTAAGCACCGCACAGATTtaat > 1:111419/1‑140 (MQ=255) gAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTtaata < 1:59526/140‑1 (MQ=255) cGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCCCAGTTTTATTAAttt > 2:60090/1‑140 (MQ=255) tAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGAtt > 1:53178/1‑140 (MQ=255) tAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGAtt > 1:18352/1‑140 (MQ=255) aaaaATAATAATTATTAAACTGCTTATAAATAACATTAACCCGTGGTAAAACGTAATCATTATTCTATAGGTCTTCCCAATTTATAGCCTAAGCCACGCCAAGTTTTAATAAGTTATCGCTTTTTAGGATTAtctactca > 2:37934/1‑138 (MQ=255) aCTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACt < 1:179559/140‑1 (MQ=255) tGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGa < 2:5299/140‑1 (MQ=255) cTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGTATTATCTTCTAATTTATCTCTTAAATGACTGata > 1:78734/1‑140 (MQ=255) aTCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGt < 1:213139/140‑1 (MQ=255) aaTAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGCTTATCTTCTAATTTATCTCTTAAATGACTGATATGTACa > 2:193738/1‑140 (MQ=255) | ATAAATATATGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACA > USA300TCH1516_ALE/1797635‑1797898 |
| 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: ‑ |