Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 2,369,488 | T→G | 6.4% | G293G (GGA→GGC) | tyrA ← | prephenate dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 2,369,488 | 0 | T | G | 6.4% | 174.3 / 3.3 | 58 | G293G (GGA→GGC) | tyrA | prephenate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (26/28); new base G (2/2); total (28/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 = 8.66e-01 |
GTTAATATCCTCTCGTGTTTCGATAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCT > BS168_reference_genome/2369350‑2369623 | gTTAATATCCTCTCGTGTTTCGATAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCt < 8:254790/139‑1 (MQ=255) ttAATATCCTCTCGTGTTTCGATAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTc > 1:159568/1‑139 (MQ=255) ttAATATCCTCTCGTGTTTCGATAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTc > 1:6330/1‑139 (MQ=255) ctcGTGTTTCGATAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGcc > 8:342345/1‑138 (MQ=255) tttCGATAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCgg > 6:29569/1‑139 (MQ=255) tttCGATAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGATCCCTTCTGCCGAAGCgg > 7:26859/1‑139 (MQ=255) aTAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGcc > 2:40467/1‑139 (MQ=255) aTGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGt < 2:159568/139‑1 (MQ=255) gATATTCGTGATACTGATGCGCTCCGCTGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGa < 6:262072/139‑1 (MQ=255) gATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGa < 1:301912/139‑1 (MQ=255) gATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGa < 1:40467/139‑1 (MQ=255) tCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATc > 8:137810/1‑139 (MQ=255) tGATACTGATGCGCTCCGCCGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGAGCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGCTaaaactt > 8:11994/1‑134 (MQ=255) tACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGc < 2:6330/139‑1 (MQ=255) gATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGtt > 6:296701/1‑139 (MQ=255) aTGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTgg < 1:39911/133‑1 (MQ=255) aTGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTgg > 2:39911/1‑133 (MQ=255) aTGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGttt > 3:164553/1‑139 (MQ=255) cTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCGGTTTTtaaa > 4:48299/1‑139 (MQ=255) atagccgTTATCTCGGAAATCACACCGGGATGATCGGGAACATCCACATATAAATCGTAAAATGCAGGTATCGCGCCCTTCTGCCGAAGCGGCAGTCCGTCGCGATAATCCTTGGCCGTTTTAAAATAACGAAATAGAt > 8:338355/7‑139 (MQ=255) aTCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTaaaa > 2:1863/1‑126 (MQ=255) aTCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTaaaa < 1:1863/126‑1 (MQ=255) tCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGAtt < 6:129740/139‑1 (MQ=255) cGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGAttt > 5:296230/1‑139 (MQ=255) ttATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCg < 5:29569/138‑1 (MQ=255) cGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGGATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCAtcttc < 8:151748/139‑1 (MQ=255) ggATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCAtcttc < 2:28048/138‑1 (MQ=255) ggATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCAtcttc < 7:360178/138‑1 (MQ=255) ggATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCAtc < 8:67941/135‑1 (MQ=255) gATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCATGGCGGTTTTAAAAAAAAGAAATAGATTATACGCAtctaatt > 8:360178/1‑135 (MQ=255) gATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCAGTGGCTGTTTTAAAAAAACGAAATATATTTTTCGCATcttctt > 1:170815/1‑139 (MQ=255) atatTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTg < 4:226755/139‑1 (MQ=255) ttACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGtt < 7:373822/138‑1 (MQ=255) ttACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTc < 8:30378/139‑1 (MQ=255) cccGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACat < 5:296701/139‑1 (MQ=255) cccGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACat < 1:364590/139‑1 (MQ=255) cccGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACat < 7:81909/139‑1 (MQ=255) ggATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCGTGGCTGTTTTTAAATAACGAAATAGATTTTCCGCATATTATTGTTCTACatatt > 6:329295/1‑138 (MQ=255) tgatCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGt > 6:127024/1‑70 (MQ=255) tgatCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGt < 5:127024/70‑1 (MQ=255) tgatCGGGAACATCCACATATAAATCGTAAAATGCAGGTATCGCGCCCTTCTGCCGAAGCGGCAGTCCGTCGCGATa > 1:271903/1‑77 (MQ=255) tgatCGGGAACATCCACATATAAATCGTAAAATGCAGGTATCGCGCCCTTCTGCCGAAGCGGCAGTCCGTCGCGATa < 2:271903/77‑1 (MQ=255) gTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTgg < 7:220841/78‑1 (MQ=255) gTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTgg > 8:220841/1‑78 (MQ=255) caTGTAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAAttt > 3:397396/1‑138 (MQ=255) atatAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCa < 5:249861/139‑1 (MQ=255) atatAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAAACCTTGGGTGTTTTAAAATAACGAAATAGATTTTACGCATCGTCTTGTACTACAAAATTACGGAGATTGTGAAATTCa > 4:307587/1‑139 (MQ=255) aTCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGaaa > 7:116529/1‑81 (MQ=255) aTCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGaaa < 8:116529/81‑1 (MQ=255) aTCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATc > 1:24845/1‑139 (MQ=255) tAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGCGTATCCGTTGGGGGTTTAAAAAAACGAAATAGATTTTACGCATCGTGTTGTTATACATATGTTCGGTTCTTGTCAATTTCACGACTCCCCt > 8:276584/1‑139 (MQ=255) tAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCAtcttc < 4:372027/95‑1 (MQ=255) tAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCAtcttc > 3:372027/1‑95 (MQ=255) tAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCa > 2:272227/1‑137 (MQ=255) tAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCa < 1:272227/137‑1 (MQ=255) cAGGTAGCGCGCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaa < 6:49289/139‑1 (MQ=255) aGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaaa < 7:342345/139‑1 (MQ=255) ggTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAAc < 1:406003/139‑1 (MQ=255) cGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCt < 1:379904/139‑1 (MQ=255) | GTTAATATCCTCTCGTGTTTCGATAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCT > BS168_reference_genome/2369350‑2369623 |
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: ‑ |