Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 1,744,607 | T→C | 100% | D22G (GAT→GGT) | apt ← | Adenine phosphoribosyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 1,744,607 | 0 | T | C | 95.8% | 68.5 / ‑3.6 | 24 | D22G (GAT→GGT) | apt | Adenine phosphoribosyltransferase |
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (16/7); total (16/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.33e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-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. |
GCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACTTTATCCAATTTTTTATTTCTGAAAAATCTTG > USA300TCH1516_ALE/1744472‑1744740 | gCTACAGGACAGCCAATGATACAGCCACGCGATTAAGGTCCTACAACGAGAACAACCTCTCTGTCTTTTTCGTATTCTAAAAAATTATCTGTTGCATAGCCATATGCTTCACCATTAACCATAATTGTAGTAATACCCtta > 2:131639/1‑140 (MQ=255) gCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTCCATAGCCATATGCTTCACCATTCTCCATAATTGTAGTAATACCCTTg > 2:37124/1‑141 (MQ=255) gACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTaa > 1:265216/1‑141 (MQ=255) aTAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACTCTTGAAACTAACCCCCGGtttg > 2:107967/1‑140 (MQ=255) ttCAGGTCCTACAACGATATCTACATCTCTGTTTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAAc < 2:281166/141‑1 (MQ=255) cGATATCAACATCTCTGTCTTTTGCATATTATACAATTTTAGCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTAtt > 1:163407/1‑141 (MQ=255) aaCATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAATTTCTGATACGTATTGCTTTaa > 1:67538/1‑141 (MQ=255) tCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAat < 1:181868/141‑1 (MQ=255) tATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCTATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTCTTGCTTTAAATCCATTAATATTttctcc > 2:113220/1‑141 (MQ=255) tATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTtcctcc > 1:178348/1‑141 (MQ=255) ttCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGa < 2:143703/140‑1 (MQ=255) cccgTTGCATCGCCATACGCTTCACCATTATCGATCATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGAc < 1:131639/138‑1 (MQ=255) cTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATAGGTATTGCTTTAAATCCATCAATATTTCCTCCTAAATTCATCACGACa > 2:275385/1‑141 (MQ=255) tGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTTATACGTATTGCTTTAAATCCATTAATATTTCCTCGTAAATTGCTGActtcaa > 2:557/1‑136 (MQ=255) tGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACaa > 2:150103/1‑141 (MQ=255) tGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATtg < 2:306381/141‑1 (MQ=255) ccATATGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACttt > 1:175438/1‑141 (MQ=255) tatGCTTCACCATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACTTTATc > 1:134696/1‑141 (MQ=255) cATTATCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGCCAACCGTGACGTTATCCAACTTTTTa > 1:274720/1‑141 (MQ=255) aTCCATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACTTTATCCAATTTTTTAtttt > 1:58617/1‑140 (MQ=255) ccATAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACTTTATCCAATTTTTTATTTCTg < 2:113669/141‑1 (MQ=255) tAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACTTTATCCAATTTTTTATTTCTGaaa < 2:183187/141‑1 (MQ=255) tAATTGTAGTAATACCCTTGAAACTAACACCTGGTTTCGGACAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACTTTATCCAATTTTTTATTTCTGaaa < 2:265216/141‑1 (MQ=255) agtaATACCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACTTTATCCACTTTTTTATTTCTGAAAAATCTTg > 2:197141/1‑141 (MQ=255) | GCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCTCCTAAATTGCTCACGACAATTGTGACTTTATCCAATTTTTTATTTCTGAAAAATCTTG > USA300TCH1516_ALE/1744472‑1744740 |
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: ‑ |