Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 3,658,834 | G→T | 37.6% | pseudogene (59/486 nt) | tuaAn ← | putative undecaprenyl‑phosphate N‑acetylgalactosaminyl‑1‑phosphate transferase; N‑terminal part of TuaA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 3,658,834 | 0 | G | T | 37.6% | 21.5 / 3.6 | 33 | pseudogene (59/486 nt) | tuaAn | putative undecaprenyl‑phosphate N‑acetylgalactosaminyl‑1‑phosphate transferase; N‑terminal part of TuaA |
Reads supporting (aligned to +/- strand): ref base G (14/4); new base T (3/7); total (20/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.04e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCTCCAAACAAATCAAAATTGAAAATACCGCAATCATTGGCAGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGCAGACGGATGTTAATTGGA > BS168_reference_genome/3658703‑3658969 | tctcCAAACAAATCAAAATTGAAAATACCGCAATCATTGGCATCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGTAt < 7:270602/139‑1 (MQ=255) tcCAAACAAATCAAAATTGAAAATACCGCAATCATTGGCAGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTGTTTTGAGATAATGTAAAAAAAAGTAttt > 3:159376/1‑139 (MQ=255) tcCAAACAAATCAAAATTGAAAATACCGCAATCATTGGCAGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTGTTCATAAAAAGTAAAAAAatttattt > 7:267580/1‑139 (MQ=255) aaTCAAAATTGAAAATACCGCAATCATTGGCAGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGTATTTGTTGAACt < 2:231545/139‑1 (MQ=255) aaaaTACCGCAATCATTGGCAGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAATAATGAATTTGTTGAACTTAGaaataac > 8:27559/1‑133 (MQ=255) aGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGTATTTGTTGAACTGAGAAATGACGACTCACATTCATGCTTTtct > 4:214780/1‑139 (MQ=255) aGCCCAATCCCCAATCCAATCAGCCCAAAAAAAATATCAATTCCGCGTTTTATTCCAAAAAACCGCCCCTTCTTTTCAGAAAATCAAAAACAATTTATTTGTTCAACTGACAATCCACGACTCACATTCATGCTTTtct < 6:124084/139‑1 (MQ=255) cgcAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCGTTTGAGATTATTTTAAAGAATGTATTTTTTGAACTGAGAATTGACGACTCACATTCATGGTTTTctct > 1:391962/1‑139 (MQ=255) cgcAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCGTTTGAGATAATGTAAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCttttgtgt > 3:275137/1‑135 (MQ=255) cgcAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCGTTTCAGATAATGTAAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTctct > 5:405378/1‑139 (MQ=255) gcATTCGCCAATCCAATGACCGCAAACCAAAAATCCATTACGCGGTTTATTGAAAAATAACGACCCGTCTTTTCATATAATTTAAAACAATGTATTTGTTGAACTGATAATTCACGACTCACATTCATGCTTTTCTCTg < 6:178315/139‑1 (MQ=255) aaTCGCAAATCCAATGAGCGCAAACAAAAAAGACATACCGCGCTTTATTCCCAAATACCGCGCCGGCTTTACAGAAAATCTAAAATAATGAATTTGTTCAACTTAGAAATCACGACTCACATTCATGCTTTTCTCTGca < 3:187887/139‑1 (MQ=255) aatccaatGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTGTTTTAAGATAATGTAAAATAATGTATTTTTTGTACTGAGAATTCACGACTCACATTCATGCTTTTGTCTGcactcacg > 8:396288/1‑137 (MQ=255) atGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCGTTTGAGATAATGTAAAAGAATGTATTTTTTGGACTGAGAATTCACGCCTCCCATTCATGCTTTTCTCTGCACACACGGacccac > 8:278644/1‑139 (MQ=255) cgcAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGTATTTTTTGAACTGAGAAATAACGACTCACAATCATGCTTTTATCTGCAATAAATGACACACACCt > 3:420080/1‑138 (MQ=255) aaaaCAAAATACAAATCACGCGATTTATTGAAAAATCACGAGCTGTATTTTCATATAATTTAAAATATTTTTTTTTTTTAACTGATAAAGCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCAc < 1:204529/138‑1 (MQ=255) acccAAATATCCTTTCCGGGTTTAATTCAAAAATAACGATCGTTCTTTTCCTATAATTAAAAATAAGTTATTTGTAGACCTGAGAATCCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGa < 2:318695/137‑1 (MQ=255) atatCCATTACGCGCTTTATTCCAAAATAACGAGCCTTCTTATAATATAATTTAAAATAATTTATTTGTTGAACTGACAATCCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGAtttgtt < 3:239288/139‑1 (MQ=255) cATTACGCGCTTTATTGCCAAATAACGCGCCGTCTTTTAAGATAATGTAAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTACCTGAACCACACCTACCGAGTTTTTTGGGa > 1:238334/1‑139 (MQ=255) ttACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATg > 3:213495/1‑139 (MQ=255) ttACGCGCTTTATTGCCAAATAACGAGCCGTCGTTTAAGATAATTTAAAAAAATGTATTTTTTTAACTGAGAATTAACGACTAACATTAATGCTTTTATATGCACTCACTGACACACCCCACACCATATGTTTTGGAtt > 6:32127/1‑138 (MQ=255) tACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATTTAAAAGAATGTATTTGTTGAACTTAGAATTAACGACTCACATTAATGCTTTTCTCTGCACACACTGACACACAACGCACAAATTGTTTTGGATGt > 8:139915/1‑139 (MQ=255) cgcTTTATTGCAAAAAAACGACCCGTATTTTAATAAAATTTAAAATATTGTTTTTTTTTAACTGATAATACACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTatg < 4:232917/139‑1 (MQ=255) cgatttATTCCAAAATACCGCGCCGTCTTTAAATATAAGCTAAAAGAATGTTTTTGTTGAACTCAAAACTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTatg < 4:293654/136‑1 (MQ=255) aaTAACGAGCCGTCTTTTCAGATAATGTAAAATAATGTATTTGTTGAACTGAGAATTCACGACTCACATTAAAGCTTTTGTCTGCACTCAATGACACACACCACACGATTTTTTTGGGATGTTATGATGACCTATGc > 3:155515/1‑137 (MQ=255) aaTAACGAGCCGTCTTTTAAGATAATGTAAAAGAATGTATTTGTTTAAATGAGAATTAACGACTCACATTCATGCTTTTATCTGCACTAACAGACACAAACCACACGATTTGTTTGGGATGTTTTTATGACCTAttc > 4:372990/1‑135 (MQ=255) aaaCGCGCCGTCTTTTCTAATATTTTAAAATAATTTTTTTGTTTATCTGAAAATCCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCcagt < 7:104089/138‑1 (MQ=255) cgCCTTCTTTTCATATAATTTAAAATAATTTTTTTGTTGAACTGAGAAAGCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGtt < 7:396288/138‑1 (MQ=255) ccttCTTTTCTTATAATTTAAAAGAAGGTTTTTGTTGAAGTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTa < 2:238334/136‑1 (MQ=255) tgttttGAGATAATTTAAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCACCTTCATGGTTTTCTCTGCACTCCCTGACACCCACCTCACGATTTTTTTGGGATGTTTTGATTACCTATGGCGGTCAGTTTAg > 4:271014/3‑137 (MQ=255) aTGTAAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGCAGACGGATGTTaa > 3:345377/1‑139 (MQ=255) aaaGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACTCACCTCACTATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGCAGACTGATGTTAATTgta > 7:113848/1‑137 (MQ=255) aaaGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGCAGACGGATGTTAATTGGa > 6:253492/1‑139 (MQ=255) | TCTCCAAACAAATCAAAATTGAAAATACCGCAATCATTGGCAGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGCAGACGGATGTTAATTGGA > BS168_reference_genome/3658703‑3658969 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |