Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 1,359,996 | T→C | 7.4% | F181F (TTT→TTC) | proCC → | redundant pyrroline‑5‑carboxylate reductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 1,359,996 | 0 | T | C | 7.4% | 154.6 / 4.0 | 54 | F181F (TTT→TTC) | proCC | redundant pyrroline‑5‑carboxylate reductase |
Reads supporting (aligned to +/- strand): ref base T (19/31); new base C (2/2); total (21/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.38e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01 |
AGCACTGGCTGCAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCGAGGAGACGCTGGAACGT > BS168_reference_genome/1359858‑1360134 | aGCACTGGCTGCAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCttt > 6:12571/1‑139 (MQ=255) cACTGGCTGCAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTga < 5:338876/139‑1 (MQ=255) tGGCTGCAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTgaag > 3:202498/1‑138 (MQ=255) ggCTGCAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTgaagaa < 1:158534/139‑1 (MQ=255) gCTGCAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCGGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAAc < 7:273445/139‑1 (MQ=255) gCTGCAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTgaagaa > 7:131883/1‑138 (MQ=255) cAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTTCAATCTTTGAAGAACTgg > 3:231971/1‑138 (MQ=255) cAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTgg > 6:152429/1‑138 (MQ=255) agagAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTga < 4:28257/129‑1 (MQ=255) agagAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTga > 3:28257/1‑129 (MQ=255) agagAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCa > 6:122235/1‑139 (MQ=255) aaaaaaGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTgg < 1:368422/139‑2 (MQ=255) aaaaaGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTc > 8:305598/1‑139 (MQ=255) acaGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTg < 1:123169/137‑1 (MQ=255) gTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGt > 8:401135/1‑136 (MQ=255) gTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAg < 2:52425/139‑1 (MQ=255) gTCTGAATGAATTGGTATAGGTGTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAg < 2:151887/139‑1 (MQ=255) atTGTTTTCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCGAGCAATTTAACGAGCTCGGCACCGGGATTCATAGGTGCAATCTTTGAAGAACTGGCAATGTCTGCCGTCAGAAACagcagc < 7:53805/139‑1 (MQ=255) atTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACagcagc < 3:184324/139‑1 (MQ=255) atTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACagcagc < 7:305598/139‑1 (MQ=255) atTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACagcagc < 3:382363/139‑1 (MQ=255) atTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACagcagc < 4:186581/139‑1 (MQ=255) atTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACagcagc < 5:152429/139‑1 (MQ=255) ttGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCGCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACagcagc < 3:291938/138‑1 (MQ=255) tATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTCGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTc > 5:72648/1‑139 (MQ=255) tATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTCGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTc > 5:66677/1‑139 (MQ=255) tCAGCCGTGTACGGGAAATCAAAGTATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCGATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAaagaag < 7:11343/138‑1 (MQ=255) tCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCa > 2:49459/1‑132 (MQ=255) tCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCa < 1:49459/132‑1 (MQ=255) gtgtACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAAATGGCACTGTCTGCCGTCAGAAACAGCAGCCTCTCCAAACAAGAAGCCtt > 5:354493/1‑139 (MQ=255) tgtACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCttt > 3:103215/1‑139 (MQ=255) ggAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGAtttt < 5:12571/139‑1 (MQ=255) aTCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGa > 1:10136/1‑139 (MQ=255) aTCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCg < 1:293662/139‑1 (MQ=255) aTCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATGCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCg < 3:399676/139‑1 (MQ=255) tCGAATTTAGATGTGGCCATCAAATTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGGCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAATCCTTTGATTTTCTGATTCATTCGc < 2:54808/139‑1 (MQ=255) aTTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCtata < 4:103215/139‑1 (MQ=255) gATGAGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGaa < 2:10136/139‑1 (MQ=255) gCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTga < 6:85924/50‑1 (MQ=255) gCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTga > 5:85924/1‑50 (MQ=255) caCGAGCTCGGCACCGGTATTCATAGAGGCCATCTTTTAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGAt < 4:231971/137‑1 (MQ=255) gAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGaa < 5:122235/139‑1 (MQ=255) gCTCGGCACCGGGATTCATAGCTGCAATCTATGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGaaa > 7:84375/1‑138 (MQ=255) tCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACagcagc > 1:286317/1‑66 (MQ=255) tCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACagcagc < 2:286317/66‑1 (MQ=255) cACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAatat < 6:398901/139‑1 (MQ=255) tCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTc > 1:296969/1‑137 (MQ=255) aTAGCTGCAATCTTCGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGTTTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCgagg < 6:66677/139‑1 (MQ=255) aTAGCTGCAATCTTCGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGTTTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCgagg < 6:72648/139‑1 (MQ=255) tAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCgagga < 7:401135/139‑1 (MQ=255) gCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCgaggag < 1:28569/138‑1 (MQ=255) cAATCTTTGAAGCACTGTCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTTATTTTCTGATTCATTCGCTATACGTAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCGAGGAGACGCt > 7:38860/1‑139 (MQ=255) cAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGattcatt > 5:387890/1‑81 (MQ=255) cAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGattcatt < 6:387890/81‑1 (MQ=255) cAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCGAGGAGACGCt > 7:293278/1‑139 (MQ=255) tGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCGAGGAGACGCTGGAACgt > 1:283010/1‑139 (MQ=255) | AGCACTGGCTGCAGAGAAAAAAGAACGTCTGAATGAATTGTTATCGGTTTTCAGCCGTGTACGGGAAATCAAAGAATCGAATTTAGATGTGGCCAGCAATTTAACGAGCTCGGCACCGGGATTCATAGCTGCAATCTTTGAAGAACTGGCACTGTCTGCCGTCAGAAACAGCAGCCTGTCCAAAGAAGAAGCCTTTGATTTTCTGATTCATTCGCTATACGGAACTGGCAAGATGCTGATCGAAAAAAATATGTCTTTCGAGGAGACGCTGGAACGT > BS168_reference_genome/1359858‑1360134 |
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: ‑ |