Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP000730 | 1,646,189 | G→C | 100% | A155G (GCA→GGA) | pepP ← | M24B family Xaa‑Pro dipeptidase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 1,646,189 | 0 | G | C | 91.5% | 32.8 / ‑5.1 | 11 | A155G (GCA→GGA) | pepP | M24B family Xaa‑Pro dipeptidase |
Reads supporting (aligned to +/- strand): ref base G (0/0); major base C (4/6); minor base A (1/0); total (5/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.55e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
TTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGA > CP000730/1646054‑1646319 | ttATCACTCGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCtttt < 1:23114/141‑1 (MQ=255) cACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTacaa < 1:12161/141‑1 (MQ=255) cacCATGTGGTAATGCGCCTCTATGACCATATGCTACAATCGTATCGAATTATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTaaac < 1:131755/141‑2 (MQ=255) aTGACCAGATGCTACAATCGTATCGAATGATGGTACATTTGCTCATAATTATAGCATTTTGCTTTTTAATATTGCTTTTAATTCTTTTTAAGTCATGCCTCCTTTTAAAACAGTTAAAATATATTAATATGTTTAATcaaa > 1:96930/1‑140 (MQ=255) gATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTa > 1:22490/1‑141 (MQ=255) gCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTAAGTCATGGCTATTTTTATCACCGTTGTAATATATTCATATGTTTCATCTACAATATTAgct > 1:82701/1‑141 (MQ=255) cAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCtt < 1:85318/141‑1 (MQ=255) aTCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAAt > 1:110663/1‑141 (MQ=255) ttGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAAtt < 2:76648/141‑1 (MQ=255) cTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTa < 2:134251/141‑1 (MQ=255) aGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGa > 2:80913/1‑141 (MQ=255) | TTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGA > CP000730/1646054‑1646319 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |