Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 2,408,513 | T→G | 25.4% | K317N (AAA→AAC) | gmmB ← | subunit of the platform for genome maintenance machinery |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 2,408,513 | 0 | T | G | 25.4% | 42.2 / 2.3 | 30 | K317N (AAA→AAC) | gmmB | subunit of the platform for genome maintenance machinery |
Reads supporting (aligned to +/- strand): ref base T (3/19); new base G (6/2); total (9/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.14e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AAAAACTGATATAACGTTTGCTGTAATTTAGTCATATTGCTGCACCTGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTGCGATTTGAAGCTGGTCTTCTGCTGATACATAGTCATTAATCATAAACGCAGGCTCATGGAGAGAAATCTCGACAATGTGGTCTT > BS168_reference_genome/2408378‑2408610 | aaaaaCTGATATAACGTTTGCTGTAATTTAGTCATATTGCTGCACCTGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTg > 6:267507/1‑139 (MQ=255) aaaCTGATATAACGTTTGCTGTAATTTAGTCATATTGCTGCACCTGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGc < 5:209316/139‑1 (MQ=255) aaaCTGATATAACGTTTGCTGTAATTTAGTCATATTGCTGCACCTGCTTTGGCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTTATTTTATTCGTCCGCATCCGGTTGGc < 1:183857/139‑1 (MQ=255) cTGATATAACGTTTGCTGTAATTTAGTCATATTGCTGCACCTGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCaaa < 8:53038/139‑1 (MQ=255) atatAACGTTTGCTGTAATTTAGTCATATTGCTGCACCTGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTTTTAGTCCGCCTCCGTTTTGGCAActc > 4:107014/1‑139 (MQ=255) ataACGTTTGCTGTAATTTAGTCATATTGCTGCACCTGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTg < 1:100100/139‑1 (MQ=255) gTTTGCTGTAATTTAGTCATATTGCTGCACCTGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTGCGAtt < 4:124898/139‑1 (MQ=255) ttttCTGTAATTTAGTCATATTGCTGCACCAGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTGCGAttt < 6:248837/139‑1 (MQ=255) tAATTTAGTCATATTGCTGCACCTGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTGCGATTTGAAGc < 3:217136/137‑1 (MQ=255) tAATTTAGTCATATTGCTGCACCTGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTGCGATTTGAAGc < 4:27970/137‑1 (MQ=255) gTCATATTGCTGCACCTGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAAATTGTTTGTTTTTATTCGTCCGCCTGCGGTTGGCAAACACTGGGATTTGACGCGGGtcgtct > 1:199526/1‑139 (MQ=255) tCATATTGCTGCACCTGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTGCGATTTGAAGCTGGTCTTctg < 5:267507/139‑1 (MQ=255) tCATATTGCTGCACCTGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTGCGATTTGAAGCTGGTCTTctg < 8:151669/139‑1 (MQ=255) gatatTGCTGCACCTGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTGCGATTTGAAGCTGGTCTTctgc < 5:162069/138‑1 (MQ=255) tatTTCTGCACCTGCTTTGTCTATGCCTGTCTGATTTGAAAGTAGCTGAACCGGTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGGTTGGCAAAATCTGCGATTTGAAGCTGGTCTTctgctg < 6:187134/139‑1 (MQ=255) ctgcACCTGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTGCGATTTGAAGCTGGTCTTCTGCTGATACa < 3:25659/139‑1 (MQ=255) ccTGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGGTTGGCAAACTCTGCGATTTGAAGCTGGTCTTGGGCTGATACATAGTc > 6:162069/1‑139 (MQ=255) gacaggTCACTGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCTAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTGCGATTTGAAGCTGGTCTTCTGCTGATACATAGTCAtt < 7:122578/134‑1 (MQ=255) gCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGGATTTGTTTGTTTTTTTTCGTCCGCATACGTTTGGCAAACTCTGCGATTTTAAGCTGGTGTTCTGGTGTTTCATATTCAtt > 2:97030/1‑139 (MQ=255) gCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTTTTAGTCCGCCTCCGGTTGGCAAACACTGCGATTTGAAGCTGGTCTTGTGGTGATACATattaat > 5:215639/1‑133 (MQ=255) cTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCTAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTGCGATTTGAAGCTGGTCTTCTGCTGATACATAGTCATTa < 6:22650/139‑1 (MQ=255) gTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTTTTAGTCCGCCTCCGGTTGGCCAACTCTTCGCGTTGAAGCTGGTCTTCTGGTGGTACATAGTCAtttaaca > 7:248270/1‑134 (MQ=255) ccAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTGCGATTTGAAGCTGGTCTTCTGCTGATACATAGTCATTAATCATAAACg < 2:66970/138‑1 (MQ=255) cAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTGCGATTTGAAGCTGGTCTTCTGCTGATACATAGTCATTAATCATAAACGCa < 1:171275/139‑1 (MQ=255) tCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGGATTTGTTTGTTTTTTTTAGTGCGCGTCCGGTTGGCAAACTCTGGGATTTGAAGGTGGTCGTGGGGTGTTAAATAGTCAGTAATAATAAAACCAGGc > 6:86233/1‑139 (MQ=255) gCTGAACCGCTGTTCCAGCCCGTCGCTAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTGCGATTTGAAGCTGGTCTTCTGCTGATACATAGTCATTAATCATAAACGCAGGCTCATGGAGAGAAAt < 7:275872/138‑1 (MQ=255) cTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGGCCGCATCCGTTTGGCAAACTCTGCGATTTGAAGCTGGTCTTCTGCTGATACATAGTCATTAATCATAAACGCAGGCTCATGGAGAGAAATCTCGACAAtg < 2:232309/139‑1 (MQ=255) gTTCCCGCCCGGCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTGCGATTTGCAGCTGGTCTTCTGCTGATACATAGTCATTAATCATAAACGCAGGCTCATGGAGAGAAATCTCGACAAtgtg < 4:272792/139‑1 (MQ=255) ttCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTTGCAAACTCTGCGATTTGAAGCTGCTCTTCTGGTGATACATAGTCATTAATCATAAACGCAGGCTCATGGACAGAAATCTCGACAAtgtg < 1:136213/138‑1 (MQ=255) gCCCGTCGCGGATTTGTTTTATTTTTTTAGTCCGCATCCGGGTGGCAAACTCTGCGATTTTAAGGTGGTGTTGTGGGGATACATATTCATTAATCATAAACGCAGGCTCATGGCTAGAAATCTCGACAAAGGGGTGtt > 6:228258/1‑138 (MQ=255) | AAAAACTGATATAACGTTTGCTGTAATTTAGTCATATTGCTGCACCTGCTTTGTCAATGCCAGTCTGATTTGAAAGTAGCTGAACCGCTGTTCCAGCCCGTCGCGAATTTGTTTGATTTTATTCGTCCGCATCCGTTTGGCAAACTCTGCGATTTGAAGCTGGTCTTCTGCTGATACATAGTCATTAATCATAAACGCAGGCTCATGGAGAGAAATCTCGACAATGTGGTCTT > BS168_reference_genome/2408378‑2408610 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |