Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 880,492 | T→C | 28.7% | intergenic (‑12/‑235) | gstB ← / → dacC | glutathione S‑transferase GstB/D‑alanyl‑D‑alanine carboxypeptidase DacC |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 880,492 | 0 | T | C | 28.7% | 8.7 / 6.6 | 14 | intergenic (‑12/‑235) | gstB/dacC | glutathione S‑transferase GstB/D‑alanyl‑D‑alanine carboxypeptidase DacC |
Reads supporting (aligned to +/- strand): ref base T (6/4); new base C (2/2); total (8/6) | |||||||||||
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.08e-02 |
ACTGTGCCGCCAGATAGCGGACAATGGCGTTTGATTCCCAAAGAATGAGATCACTTTCGTCGTCACGCAACAACGGCACCAGCCCGTTAGGGTTCATCGCCAGAAAATCAGCATCGTGATTTATCCCAAACTCACGGCCCGCGAGAATTTGCTCATAAGGTAGTTCCAGTTCTTCGAGCGTCAGCAATACTTTTTTTACGTTAGTTGAATTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCA‑GCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTAACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAACTTACGGTTTCTTTAAGTTTGTGTGTGCGTTATTAATCACCAAACTTATCATACGGCGATATAACGTATTTTTTTTGAATGGATACTCGGGTGGCATTTATGACGCAATACTCCTCTCTCCTTC > NC_000913/880246‑880751 | acggtgCCGCCAGAAACGGGACAATGGCGTTTGATTCCAAAAAAATGAGATCACTTTCGTCGTCACGCAACAACCGCACCAGCCCGTTAGGGTTCATCGCCAGAAATTCAGCATCGTGATTTATCCCAAACTCACGGCCCGCGAGAATTTGCTCATAAGGTAGTTCCAGTTCTTCGAGCGTCAGCAATACTTTTTTTACGTTAGTTGAATTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCA‑GCCACAGGCTGAACAg < 1:167647/271‑1 (MQ=255) ccAGATAGCGGACAATGGCGTTTGATTCCCAAAGAATGAGATCACTTTCGTCGTCACGCAACAACGGCACCAGCCCGTTAGGGTTCATCGCCAGAAAATCAGCATCGTGATTTATCCCAAACTCACGGCACGCGAGAATTTGCTCATAAGGTAGTTCCAGTTCTTCGAGNGTCAGCAATACTTTTTTTACGTTTGTTGAATTATTCCGACCCCACAGCTTAATCATAATCACCCCTTTTCCGACAGGGACGCCCCAGGCGGAAAAgagcacaggg > 2:167646/1‑266 (MQ=255) tcgtcACGCAACAACGGCACCAGCCCGTTAGGGTTCATCGCCAGAAAATCAGCATCGTGATTTATCCCAAACTCACGGCCCGCGAGAATTTTCTCATAGGGTAGTTCCAGTTCTTCGAGCGTCAGCAATACTTTTTTTACGTTAGTTGAATTATTCCGACCCCACCGCGTATTCCTACTCCCCCCTTTTCCGACAGGCA‑Gcccaagga > 1:220799/1‑202 (MQ=255) tttttgACGTTAGTTGAATTATTCCGACCCCACAGCGTAATCATACTCACCCCGTTTCCGACAGGCA‑GCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTAACTTTTAAGCAACAGCTCGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAACTTACGGt < 1:4513/199‑1 (MQ=255) aCCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCA‑GCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTAACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGGATAAACTTTAAAAAGTTTACCAACTTa > 1:388925/1‑169 (MQ=255) aCCCCACAGCGTAATCATACTCACCCCTTCTCCGACAGGCA‑GCCACAGGCTGAACAGGACTcc > 1:231312/1‑63 (MQ=255) aCCCCACAGCGTAATCATACTCACCCCTTCTCCGACAGGCA‑GCCACAGGCTGAACAGGACTcc > 1:231319/1‑63 (MQ=255) cacaGCGTAATCATACTCACCCCTTCTCCGACAGGCA‑GCCACAGGCTGAACAGGACTCCGATg < 2:231311/63‑1 (MQ=255) cacaGCGTAATCATACTCACCCCTTCTCCGACAGGCA‑GCCACAGGCTGAACAGGACTCCGATg < 2:231318/63‑1 (MQ=255) acaGCGTAATCATACTCACCCCTTTTCCGACAGGCA‑GCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTAACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAAATTTAAAAAATTTACAAACTTACGGTTTCTTTAAGTTTGTGTGTGCGCTATTAATCACCAAACTTATCATCCGGCGATATAACGTATTTTTTTTGAATGGATACCCGGGTGGTATTTATGACTCAATActcc > 1:88998/1‑274 (MQ=255) tCATACTCACCCCTTTTCCGACAGGCA‑GCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTAACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAACTTACGGTTTCTTTAAGTTTGTGTGTGCGTTATTAATCACCAAACTTATCATACGGCGATATAACGTATTTTTTTTGAATGGATACTCGGGTGGCATTTATGACGCAATACTCCTCTCTCCTTc > 1:30601/1‑275 (MQ=255) cATACTCACCCCTTTTCCGAAAGGGA‑GCCACATGCTGAACAAGACTACGATGGTGAAGTAAAGTTAACTTGTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAACTTACGGTTTCTTTAAGTTTGTGTGTGCGTCAGTAATCACCAAACTTATCATACGGCGATATAACGTATTTTTTTTGAATGGATACTCGGGTGGCATTTATGACg < 1:174992/255‑1 (MQ=255) cTCACCCCTTTTCCGACAGGCA‑GCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTAACTTTTAAGCAACAGCTGGCAAAAAAATGATACCAGAAGCAGCGAAGCGGCCGGCTATTGCATAAAc > 2:12443/1‑129 (MQ=255) ccccTTTCCCGACAGGCA‑GCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTAACTTTTAAGCAACAGCTGGCAAAAAAAAGATACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTa < 1:12443/129‑1 (MQ=255) | ACTGTGCCGCCAGATAGCGGACAATGGCGTTTGATTCCCAAAGAATGAGATCACTTTCGTCGTCACGCAACAACGGCACCAGCCCGTTAGGGTTCATCGCCAGAAAATCAGCATCGTGATTTATCCCAAACTCACGGCCCGCGAGAATTTGCTCATAAGGTAGTTCCAGTTCTTCGAGCGTCAGCAATACTTTTTTTACGTTAGTTGAATTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCA‑GCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTAACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAACTTACGGTTTCTTTAAGTTTGTGTGTGCGTTATTAATCACCAAACTTATCATACGGCGATATAACGTATTTTTTTTGAATGGATACTCGGGTGGCATTTATGACGCAATACTCCTCTCTCCTTC > NC_000913/880246‑880751 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |