| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | BS168_reference_genome | 2,910,929 | A→G | 6.4% | intergenic (‑183/+187) | lysCA ← / ← uvrC | aspartokinase II alpha subunit (aa 1‑>408)/excinuclease ABC (subunit C) |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | BS168_reference_genome | 2,910,929 | 0 | A | G | 6.4% | 189.7 / 3.4 | 63 | intergenic (‑183/+187) | lysCA/uvrC | aspartokinase II alpha subunit (aa 1‑>408)/excinuclease ABC (subunit C) |
| Reads supporting (aligned to +/- strand): ref base A (30/28); new base G (2/2); total (33/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01 | |||||||||||
AAAAAGCAATGAGAGGAATACTCTCATTGCTTATCAATTAATCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCGCACCTCTATCTTCACCATATAACACGATC > BS168_reference_genome/2910811‑2911064 | aaaaaGCAATGAGAGGAATACTCTCATTGCTTATCAATTAATCTTCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGa > 3:271221/1‑139 (MQ=255)aaaaaGCAATGAGAGGAATACTCTCATTGCTTATCAATTAATCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGa < 3:10570/139‑1 (MQ=255) gCAATGAGAGGAATACTCTCATTGCTTATCAATTAATCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATgggg > 3:81250/1‑139 (MQ=255) aTACTCTCATTGCTTATCAATTAATCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCt > 2:144887/1‑138 (MQ=255) tACTCTCATTGCTTATCAATTAATCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAg < 4:189044/108‑1 (MQ=255) tACTCTCATTGCTTATCAATTAATCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTAATTCAATGCACGGCCAg > 3:189044/1‑108 (MQ=255) aTTGCTTATCAATTAATCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACg > 5:160361/1‑139 (MQ=255) aTTGCTTATCAATTAATCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACg > 3:65415/1‑139 (MQ=255) aTTGCTTATCAATTAATCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACg > 1:39963/1‑139 (MQ=255) aTTGCTTATCAATTAATCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACg > 3:216589/1‑139 (MQ=255) gCTTATCAATTAATCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCt < 1:169975/138‑1 (MQ=255) tATCAATTAATCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCt < 8:8830/135‑1 (MQ=255) tATCAATTAATCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCt > 7:8830/1‑135 (MQ=255) tATCAATTAATCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTcccc > 8:297196/1‑139 (MQ=255) aTTAATCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCa < 8:197247/139‑1 (MQ=255) aaTCTTCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGc < 4:271221/138‑1 (MQ=255) aaTCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATAACGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGcc > 6:291374/1‑139 (MQ=255) catAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCGCGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTc > 3:232151/1‑139 (MQ=255) aTTATGATCAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCTATGCCCGGCAAGCAAATAATACCGATGGGGTTGTATTTGCTTCGGCGCCGCTCCCCTTTCAGCGTTGTTCACAGCATc > 6:307827/1‑139 (MQ=255) aTTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCa > 2:256771/1‑72 (MQ=255) aTTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCa < 1:256771/72‑1 (MQ=255) ttATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCGCGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATc < 4:232151/138‑1 (MQ=255) tATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCa > 5:55221/1‑139 (MQ=255) tATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCa > 1:48395/1‑139 (MQ=255) tATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCa > 5:232481/1‑139 (MQ=255) caacGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTtc < 4:216589/139‑1 (MQ=255) caacGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTtc > 6:210100/1‑139 (MQ=255) acGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTtctc > 6:61756/1‑139 (MQ=255) cGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTtctc > 6:202461/1‑138 (MQ=255) gCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCa < 4:65415/139‑1 (MQ=255) gCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCa < 2:191008/139‑1 (MQ=255) gCCCACCAAGAAAATGATTTCTTGACAGCCTTACATTTATTAAAAGCACGGCCAGCAAATAATACCGATGTGGTATTATTTGCTTCGGCGACGCTCGCATTTCAGCGTTTTTGACCGATTCCATCTTTCTCCAAAGGCa < 5:307827/139‑1 (MQ=255) ccAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCaaa < 2:62430/54‑1 (MQ=255) ccAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCaaa > 1:62430/1‑54 (MQ=255) aaGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCtt < 6:232481/139‑1 (MQ=255) aTTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCg < 2:39963/138‑1 (MQ=255) aTTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCg < 3:102244/138‑1 (MQ=255) aTTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCg < 5:202461/138‑1 (MQ=255) cTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCGCAcc < 6:55221/138‑1 (MQ=255) gACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCGCACct < 1:144887/136‑1 (MQ=255) gCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCGCACCTCTATCtt < 5:61756/139‑1 (MQ=255) cTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCGCACCTCTATCTTCa < 5:291374/139‑1 (MQ=255) cTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCGCACCTCTATCTTCa < 4:81250/139‑1 (MQ=255) tACATTTATTCAATGCGCGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCa < 1:261492/108‑1 (MQ=255) tACATTTATTCAATGCGCGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCa > 2:261492/1‑108 (MQ=255) aTTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGCCGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCAAAAGGCATACTCTTGAAGTTCGCACCGCTATCTTCACCATATaa > 8:138886/1‑138 (MQ=255) aTTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATc > 5:239785/1‑83 (MQ=255) aTTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATc < 6:239785/83‑1 (MQ=255) aTTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCGCACct > 3:130476/1‑121 (MQ=255) aTTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCGCACct < 4:130476/121‑1 (MQ=255) aTTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCGCACCTCTATCTTCa > 1:181036/1‑130 (MQ=255) aTTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCGCACCTCTATCTTCa < 8:15657/130‑1 (MQ=255) aTTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCGCACCTCTATCTTCa > 7:15657/1‑130 (MQ=255) aTTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCGCACCTCTATCTTCa < 2:181036/130‑1 (MQ=255) aTGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGTCATACTCTTGAAGTTCGCACCTATATCTTCACCATATAACAAGa > 6:191283/1‑138 (MQ=255) aTGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCGCACCTCTATCTTCACCATATAACACGAt > 1:76145/1‑139 (MQ=255) tGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCgg > 1:208526/1‑44 (MQ=255) tGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCgg < 2:208526/44‑1 (MQ=255) tGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGa > 3:186908/1‑74 (MQ=255) tGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGa < 4:186908/74‑1 (MQ=255) tGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCa < 4:289571/95‑1 (MQ=255) tGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCa > 3:289571/1‑95 (MQ=255) tGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCGCACCTCTATCTTCACCATATAACACGATc > 4:289619/1‑139 (MQ=255) | AAAAAGCAATGAGAGGAATACTCTCATTGCTTATCAATTAATCATCATAAATGATTATGAACAACGAGATAGCCCTCCAAGAAAATGATTTCTTGACAGCCTTACATTTATTCAATGCACGGCCAGCAAATAATACCGATGGGGTTTTATTTGCTTCGGCGACGCTCCCCTTTCAGCCTTTTTCACAGAATCCATCTTTCTCCAAAGGCATACTCTTGAAGTTCGCACCTCTATCTTCACCATATAACACGATC > BS168_reference_genome/2910811‑2911064 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |