Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 2,818,348 | G→T | 24.4% | P5Q (CCA→CAA) | yrzK ← | conserved protein of unknown function |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 2,818,348 | 0 | G | T | 24.4% | 109.7 / 7.2 | 53 | P5Q (CCA→CAA) | yrzK | conserved protein of unknown function |
Reads supporting (aligned to +/- strand): ref base G (11/28); new base T (12/2); total (23/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.20e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCCTCAGCGCTGCTTCTTTGGACATAAGACCAAATTCAGATGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGCTTAGTATAAACAAATCATTTTATTTTTATTTAAAAATCAGAAAGGACTTTACGATATTAAACAAGAATGCTTTTTTAATGATCAT > BS168_reference_genome/2818214‑2818477 | tCCTCAGCGCTGCTTCTTTGGACATAAGACCAAATTCAGATGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGgat < 1:43905/138‑1 (MQ=255) tCCTCAGCGCTGCTTCTTTGGACATAAGACCAAATTCAGATGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGgat < 6:323489/138‑1 (MQ=255) cgcTGCTTCTTTGGACATAAGACCAAATTCAGATGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGGAATTCCC‑GTGCCCGTTTTTGATGGTAtttg > 3:292221/1‑138 (MQ=255) ttGGACATAAGACCAAATTCAGATGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATGTTTTGAGGTTTTTTAATGATTTgga > 4:320254/1‑137 (MQ=255) ggACATAAGACCAAATTCAGATGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGcac < 5:165820/138‑1 (MQ=255) cATAAGACCAAATTCAGATGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGACC‑CTCCCAATTTTTGATTGTTTTTAATGTTTTTGACACCt > 7:123541/1‑139 (MQ=255) aaGACCAAATTCAGATGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCGTGTTTGGATGAAATTTCATGATTTGCACACCtttt > 6:59007/1‑139 (MQ=255) ccAAATTCAGATGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGAt < 1:291653/139‑1 (MQ=255) ccAAATTCAGATGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGAt < 5:60101/139‑1 (MQ=255) ccAAATTCAGATGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGAt < 3:36635/139‑1 (MQ=255) ccAAATTCAGATGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGAt < 3:320254/139‑1 (MQ=255) aataaTTTAGATGCTGTTTCAGTAAGATAACCAGTCGTATCGGGATCTGCGTGTTCGCTGTTATGCAGATTCGCTTCATCGTCGTTAAATGCAC‑GTCCCATTTTTGTATGATATTTCATGATTTGCAGAACTTTTTGa < 4:237301/135‑1 (MQ=255) aTTCAGATGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATTCAC‑GTCCCATTTTTTGATTATATTTAATTATTTTCCAACCATTTTTTTTGTa > 8:71356/1‑139 (MQ=255) aTTCAGATGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGAAC‑GTCCCAATTTTGGTTTATATTTCATGAATTGCACACCATTTTGaattta > 4:279219/1‑134 (MQ=255) gATGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGGAATTCAC‑GGGCCAATTTTGGAGTATATTTAATGTTTTTCGCACCGTTTTTATTGTtttttc > 8:313502/1‑137 (MQ=255) tGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGc < 1:184990/139‑1 (MQ=255) tGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGc < 4:57608/139‑1 (MQ=255) tGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGc < 5:305926/139‑1 (MQ=255) tGCTGTTTCCGTAAGATAAACAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATTCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGc < 3:308830/139‑1 (MQ=255) cTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGCt < 1:131580/138‑1 (MQ=255) cTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGCt < 3:342794/138‑1 (MQ=255) cTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGCt < 7:71356/138‑1 (MQ=255) cTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGCt < 6:186082/138‑1 (MQ=255) tAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCAGTTTTTGTTGATATTTAATGATTTGCACACCATTTTTATTTTATTTACGCTTAGTATaaa > 8:222023/1‑138 (MQ=255) tCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAAAGCAC‑GTACCATTTTTGGATGATATTTAATGATTTTCACACCTTTTTTTTTGTATTTACGCTTAGTATAAACAAATAattttat > 4:86807/1‑139 (MQ=255) cGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGGAATTCAC‑CTCCCCTTTTTTGTTGTTTTTTTATTTTTTTGGCACCTTTTTTTTTGTATTTTAGCTTAGTATAAACAAATAaatttttt > 3:202494/1‑139 (MQ=255) aTCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTTGATGATATTTCATGATTTGcac > 5:253448/1‑84 (MQ=255) aTCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGcac < 6:253448/84‑1 (MQ=255) aTCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTAATTATTTGCCCACCTTTTTGATTGTATTTAAGCTTTGTATAAACAAATAATTTTATTtttattta > 4:267314/1‑138 (MQ=255) tGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATTAAC‑ATACCCTTTTTTTAATATATTTAATTATTTTCAAACATTTTTTATTGTTTTTAAAATTAATATAAAAAAAAAAATTTATTTTTAATTtaaa > 4:192039/1‑138 (MQ=255) cGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGGAATGGGCGGGCCCGTTTTTGG‑TGTTTTTTGCTGGTTTGCGCCCCGGTTTGTTTTGTTTTTCGCTTTGTGTAAACAAATCATTTTTTTTTTTTTTTAAAAt > 5:46702/1‑138 (MQ=255) aTCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTTGATGATATTTCATGATTTGCACACCTTTTTGAt > 2:56578/1‑86 (MQ=255) aTCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGAt < 1:56578/86‑1 (MQ=255) aTCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGCTTAGTATAAACAAATCATTTTATTTTTATTTAAAAATCAGaa < 8:44315/139‑1 (MQ=255) aTCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGCTTAGTATAAACAAATCATTTTATTTTTATTTAAAAATCAGa < 3:279219/138‑1 (MQ=255) tCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGc > 8:50512/1‑96 (MQ=255) tCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGc < 7:50512/96‑1 (MQ=255) cTGTAATGCAGATTCGCTTCACCGTCGTGAAAAGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGCTTAGTATAAACAAATCATTTTATTTTTATTTAAAAATCAGAAAGGa < 3:86807/139‑1 (MQ=255) gTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGAt < 4:144085/80‑1 (MQ=255) gTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGAt > 3:144085/1‑80 (MQ=255) aTGCAGATTCGCTTCTCCGTAGTGATATGGAC‑GTTCCATTTTTTGTTGATATTTCATGATTTGCACAGCTTTTTGATTGTATTTACGATTATTATAAACAAATCATTTTATTTTTATTTAAAAATCAGAAAGTACTTTa < 4:115685/139‑1 (MQ=255) aTGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTTGATGATATATAATGATTTGCACACCTTTTTGATTGTATTTACGCTTAGTATAAACAAATCATTTTATTTTTATTTAAAAATCAGAAAGGACTTTa < 3:192039/139‑1 (MQ=255) aTGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGCTTAGTATAAACAAATCATTTTATTTTTATTTAAAAATCAGAAAGGACTTTa < 7:46194/139‑1 (MQ=255) tGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGCTTAGTATAAACAAATCATTTTATTTTTATTTAAAAATCAGAAAGGACTTTAc < 7:313502/139‑1 (MQ=255) aTTCGCTTCACCGTCGTGTAATGCAC‑CTACCAATTTTTGATGTTATTTAATTAATTTCACACCTTATTTTTTGTATTTAAGCTATTTTAAAACAAATAATTTTATTTTTATTTAAAAAAAAAAAATGACTTTActaaaa > 4:207633/1‑134 (MQ=255) aTTCGCTTCACCGTCGTGGAATGCAC‑GGGCCCTTTTTTGTTGATATTTGATGATTTTGACACCTTTTTGATTGTATTTTCGCTTAGTATAAACAAATCATTTTATTTTTTTTTAAAAATAATCAAGGACATTACGattt > 2:260625/1‑137 (MQ=255) aTTCGCTTCACCGTCGTGAAATTGGC‑GTACCATTTTTTGATGATATTTCAATAATTGCAAACCATTTTTATTTTATTTAAGAATAATATAAACAAAAAAATTTAATTTTTTTATAAAATAAGAAAAGACATTACGatat > 4:349208/1‑139 (MQ=255) cAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGc > 4:173589/1‑59 (MQ=255) cAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGc < 3:173589/59‑1 (MQ=255) cAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGCTTAGTATAAACAAATCATTTTATTTTTATTTAAAAATCAGAAAGGACTTTa > 7:339487/1‑110 (MQ=255) cAC‑GTACCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGCTTAGTATAAACAAATCATTTTATTTTTATTTAAAAATCAGAAAGGACTTTa < 8:339487/110‑1 (MQ=255) c‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGCTTAGTATAAACAAATCATTTTATTTTTATTTAAAAATCAGAAAGGACTTTACGATATTAAACAAGAATGCTTTTTTAATGat < 7:228477/139‑1 (MQ=255) aTTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGCTTAGTATAAACAAATCATTTTATTTTTATTTAAAAATCAGAAAGGACTTTACGATATTAAACAAGAATGCTTTTTTAATGatcat > 7:118020/1‑136 (MQ=255) | TCCTCAGCGCTGCTTCTTTGGACATAAGACCAAATTCAGATGCTGTTTCCGTAAGATAACCAGTCGTATCGGGATCTGCGTGATCGCTGTAATGCAGATTCGCTTCACCGTCGTGAAATGCAC‑GTCCCATTTTTGGATGATATTTCATGATTTGCACACCTTTTTGATTGTATTTACGCTTAGTATAAACAAATCATTTTATTTTTATTTAAAAATCAGAAAGGACTTTACGATATTAAACAAGAATGCTTTTTTAATGATCAT > BS168_reference_genome/2818214‑2818477 |
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: ‑ |