Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 166,037:1 | +T | 100% | intergenic (+4/‑27) | trnI‑Gly → / → trnI‑Arg | tRNA‑Gly/tRNA‑Arg |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 166,037 | 1 | . | T | 97.0% | 79.4 / ‑3.6 | 33 | intergenic (+4/‑27) | trnI‑Gly/trnI‑Arg | tRNA‑Gly/tRNA‑Arg |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base T (13/19); total (14/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.24e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
ATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGC‑TTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTA‑CCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCATATCTTTTAATAGAATAGATAGGAAATCG > BS168_reference_genome/165906‑166169 | gacagttattaCCAAGGTTTCTCATAAGGAGAAAGC‑TTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTaa < 1:82602/133‑1 (MQ=255) tattatTACCAAGGTTTCTCATAAGGAGAAAGC‑TTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACggg < 3:188274/139‑1 (MQ=255) cAAGGTTTCTCATAAGGAGAAAGCTTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACGGGATATTGt > 5:198546/1‑138 (MQ=255) aaGGTTTCTCATAAGGAGAAAGC‑TTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACGGGATAtt < 6:102825/134‑1 (MQ=255) aaGGTTTCTCATAAGGAGAAAGC‑TTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACGGGATAtt > 5:102825/1‑134 (MQ=255) aaGGTTTCTCATAAGGAGAAAGC‑TTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACGGGATATTGTaaa > 3:109331/1‑139 (MQ=255) aGGTTTCTCATAAGGAGAAAGC‑TTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACGGGATATTGTAaat < 7:129643/139‑1 (MQ=255) aaaaGGAGAAAGC‑TTTATTTATTGCGATATGCGTAAGTATTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGTTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACGGGATATTGTAaattgaattg < 4:216629/135‑1 (MQ=21) aGGAGAAAGC‑TTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACGGGATATTGTAAATTGAATTGCGccc > 6:190629/1‑139 (MQ=255) gaAAGC‑TTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAg < 3:225399/139‑1 (MQ=255) gaAAGC‑TTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGGTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAg < 5:164336/139‑1 (MQ=255) aaaGC‑TTTTTTTATTGCGATATGCGGATGTAGTTCAGTGGTAGAACACGACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAg < 4:225264/138‑1 (MQ=255) aaGC‑TTTTTTTATTGCTATATGCGGTAGTAGTTCAGTGGTATAACAGCTACTTGCCTAGTTGGGGGTCGCGGGTTCGAATCCCGTCTTCATCTCCA‑TTTATCCTAAGGGGATTTTGTAAATTGAATTGCGCCCGTAGCt < 4:30424/139‑1 (MQ=21) agcttttttttATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCt < 6:198546/135‑1 (MQ=255) c‑TTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCaa < 4:109331/139‑1 (MQ=255) ttGCGTTATGCGGAAGTAGTTCAGTGGTAGAACTCCAGCTTGCCAAGGTGGGGGTCGCGGGATCGAATCCCGTCTTGCGCTCCA‑TTTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATaga < 2:149200/139‑1 (MQ=21) gtTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGGCGTCCGCTCCATTTTT‑CCTAACGGGATATTTTAAATTGAATTGCGCCCGTAGAGCAAATGGATAGAGAGTTTGACTAAGGagc > 2:162538/1‑136 (MQ=21) gtTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCCTCGTCCGCTCCA‑TTTATCCAAACGGGATATTGTAAATTGAATTGCGCCCGTTCCTCAAATGGTTAGAGCGTGTGACTTCGGATc > 7:62471/1‑138 (MQ=21) gTAGAACACCACCTTGCCAAGGGGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGtt < 5:154360/138‑1 (MQ=255) caccTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGa < 5:190629/139‑1 (MQ=255) caccTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCGTCCGCGCCA‑TTTATGCGAACGGGATATTGTAAATTTAATTTCGCCCGTTGCGCAATTGGTTAGAGCGGTGGGCTACGGATCAAAAGGGTATGGGGTCgc > 2:37652/1‑138 (MQ=21) ccTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGAct > 7:179160/1‑139 (MQ=255) ccTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTACGCTCAA‑ATTATACTAACGGGATTTTGTTAATTTAATTTCGCCCCTAGCTGAATTGGAATGAGCGTTTGACAACGGATCAAAAGGTTAGGGGTTCgcag > 6:46890/1‑136 (MQ=21) ccAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATACCAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCctct > 6:80858/1‑139 (MQ=255) ccAAGGTGGGGGTCGCGGGGTCGAATCCCGGCGTCAGCTCCA‑ATTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCCAATGGATAGAGCGTTTGGCTACGGCTCAAAATGTTAGGGGTTGGActcctgg > 2:34070/1‑137 (MQ=37) aGGTGGGGGTCGCGGGTTCGAATCCCGGCTTCCGCGCCA‑TTTATCCGTACGGGATATTGTAAATTTAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGGCGACGGATCAAAAGGGTTGGGGTTCGCCTCCTCTCgg > 8:8411/1‑139 (MQ=21) cgcgGGTTCGAATCCCGTCTTCCGCTCCA‑TTTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCata < 1:220848/139‑1 (MQ=255) gggTTCGAATGCGGTCTTCCGCTCCA‑TTTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGATCGACTCCTCTCGGGCGCGCCATATCt < 5:53912/139‑1 (MQ=21) ccGCTCCA‑TTTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAgggg < 2:70070/93‑1 (MQ=255) ccGCTCCA‑TTTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAgggg > 1:70070/1‑93 (MQ=255) ccGCTCCA‑TTTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCGCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCATATCTTTTAATAGAatagat < 6:105335/136‑1 (MQ=255) cA‑TTTATCCTAACGGGTTATTGTTAATTGAATTGCGCCCGTAGCTCAATTGGATTGAGCGTTTGACTACGGATCAAAAGGTTAGGGGATCGACTCCTCTCGGGCGCGCCATATCTTTTAATAGATTAGATAGGAAATCg < 4:165795/139‑1 (MQ=14) cA‑TTTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCATATCTTTTAATAGAATAGATAGGAAATCg > 5:186816/1‑139 (MQ=255) | ATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGC‑TTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCA‑TTTA‑CCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCATATCTTTTAATAGAATAGATAGGAAATCG > BS168_reference_genome/165906‑166169 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |