Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 608,164 | C→A | 29.5% | intergenic (‑89/+82) | ydgF ← / ← dinB | putative amino acid permease/nuclease inhibitor |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 608,164 | 0 | C | A | 29.5% | 46.9 / 7.3 | 29 | intergenic (‑89/+82) | ydgF/dinB | putative amino acid permease/nuclease inhibitor |
Reads supporting (aligned to +/- strand): ref base C (17/3); new base A (2/7); total (19/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.14e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCTGCAACGTTTGCTGATTTATATTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTC‑CC‑G‑AAAA‑A‑GCGGCTCCATTTAGGGAGCCGC‑T‑TTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCG > BS168_reference_genome/608027‑608295 | tctGCAACGTTTGCTGATTTATATTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTAATTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTAc > 7:162578/1‑138 (MQ=255) cGTTTGCTGATTTATATTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCCATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTTAAAAAg > 4:79695/1‑139 (MQ=255) tttGCTGATTTATATTGTCTTTCGTCATGTCGTCTGTCACTTTTTGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTTATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGa > 3:183372/1‑139 (MQ=255) tttGCTGATTTATATTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGa > 2:245935/1‑139 (MQ=255) ttGCTGATTTATATTGTCTTTCGTCATGTCGTCTGTGACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATTACCCTAAAATTTCTATCTTTTTTACTAAAAAAGGAt > 8:19035/1‑139 (MQ=255) gCTGATTTATATTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATATATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGAtt > 1:274494/1‑138 (MQ=255) atataTTGTCTTTGGTAATGTCGTCTGTCACTTTATGTGTTGCTCCATTGTATATTGTTTTTTGTTTGAGCATTCTATTTTTAATGTTTTTTATAACACTAATATTTCTTTCTTTTTTACTAAAAAAGGATTGATGaa < 7:156166/137‑1 (MQ=255) tCTTTCGTCATGTCGTCAGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTTAGTAAGCAATCTATTTTTAATGATTATTATAACAATAAAATTTAAAAATTTTTTAATAAAAAAAGATTTATGAAAAGAAttt > 2:9563/1‑138 (MQ=255) tttCGTCATGTCGTCTGTCAATTTATGTGTTCCTCCATTGAAAATTGTTTATTGATTCAGCAATCTATTTTTAATGATTATTATAACAATAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAATAATTTTaa > 5:279136/1‑139 (MQ=255) tCGTCATGTCGTCTGTCACTTTATGTGTTACTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAAt < 1:9563/138‑1 (MQ=255) cGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAATGATTGATGAAAAGAATTTTAAtgt > 8:39717/1‑139 (MQ=255) cGTCATGTCCTCGGAAAATTTATGTTTTGCTCCATAGAAAAATAATTATTAAGTCATAAATCTATATTTAATGATTATTATAAAAAAAAAAAAACAATATTTTTTAATAAAAAATTATTGATTAAAAAAATTTaaaat > 8:156166/1‑136 (MQ=255) gtcgtcTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTg > 6:303523/1‑139 (MQ=255) cctccATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTC‑CC‑G‑aaaa‑aagcgg > 5:286899/1‑134 (MQ=255) cAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAATAATTTTAATGTGATTTGGACTATTTC‑CC‑GAAAAA‑A‑GCGGCTCCATTTAGGGAGCAGC‑G‑TTTTTAt > 1:30107/1‑139 (MQ=255) cAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTC‑CC‑GAAAAA‑A‑GCGGCTCCATTTAGGGAGCCGC‑g‑ttttttt > 3:54538/1‑137 (MQ=255) ttttAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTC‑CC‑GAAAAA‑A‑GCGGCTCCATTTAGGGAGCCGC‑G‑TTTTTTTTATTTTTAAGt > 8:84402/1‑139 (MQ=255) ttAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTC‑CC‑GAAAAA‑A‑GCGGCTCCATTTAGGGAGCCGC‑T‑TTTTTTTTGTTTTTTAtttt > 3:200792/1‑139 (MQ=255) ttattaTAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTC‑CC‑GAAAAA‑A‑GCGGCTCCATTTAGGGAGCCGC‑TGTTTTTTTGTTTTTTATGttttttga > 7:128646/1‑135 (MQ=255) ttattaTAACAATAAAATTTCTATCTTTTTTAATAAAAAAGGATTTATCAAAATAATTTTAATGTGATTTGCACTATTTC‑CCTA‑AAAA‑A‑GCGGCTCCATTTAGGGAGCCGC‑T‑TTTTTATGCTTTTTCAGTTTTTGTCa < 5:303523/138‑1 (MQ=255) tattatAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTC‑CC‑GAAAAA‑A‑GCGGCTCCATTTAGGTAGACGCTT‑TTTTTTTTATTTTTAAGttttttgat > 5:123362/1‑134 (MQ=255) aattaaAACAAAAAAATTTATATATTTTTAAAAAAAACATGATTTATCAAAATATTTTATATGTTATTTGGAATATCCC‑CC‑C‑AAAACC‑GCGGCTCCATTTAGGGAGCCGC‑T‑TTTTTATGCTTTTTCAGTTTTTGTCATa < 2:181064/138‑1 (MQ=255) tttatAAAAAAAAAATTTTTATATTTTAAAATAAAAATGTTTTTATAAAAATATTTTTATTGTTATTTGCAAAATACCACC‑A‑AAAA‑C‑GCGGCTCCATTTAGGGAGCCGC‑T‑TTTTTATGCTTTTTCAGTTTTTGTCATa < 5:213617/137‑1 (MQ=255) ttaaAACAATAAAATTTATTTATTTTTTAATAAAAAAGTTTTTATGAAAAGAATTTTAATGGGTTTTGGACTATTTC‑CC‑T‑AAAA‑ACGCGGCTCCATTTAGGGAGCCGC‑T‑TTTTTATGCTTTTTCAGTTTTTGTCATaa < 6:279136/138‑1 (MQ=255) aaaaaaaaaTTTATATATTTTTTAATAAAAATTTATTGTTTAAAATATTTTAAATTGCATTGGGCCTATTTC‑CC‑C‑AAAA‑ACGCGGCTCCATTTAGGGAGCCGC‑T‑TTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGt < 2:30107/134‑1 (MQ=255) aaaaaatTTTATATTTTTTTAAAAAAAAATTTATTTATGAAATTATTTTAAATGTTTTTTGTAATTTTAC‑CC‑C‑AAAACA‑GCGGCTCCATTTAGGGAGCCGC‑T‑TTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTag < 5:187657/136‑1 (MQ=255) aaaTTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTC‑CC‑GAAAAA‑A‑GCGGCTCCATTTAGGGAGCCGC‑T‑TTTTTATTATTTTTTAGTTTTTTTAATAAAAAGGTATAGcc > 3:306623/1‑139 (MQ=255) aaaTTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTC‑CC‑GAAAAA‑A‑GCGGCTCCATTTAGGGAGCCGC‑G‑TTTTTTTTGTTTTTTATTTTTTTTGATAAAAATGTAGAGcc > 3:39119/1‑139 (MQ=255) tttttaAATAAAAATTTATTTATTAAATTAATTTTATTGTTTTTTGAACTATTCC‑CC‑T‑AAAA‑ACGCGGCTCCATTTAGGGAGCCGC‑T‑TTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGt < 4:39119/138‑1 (MQ=255) tttaACAAAAAAAGTATTTATCAAAATAATTTTAATGTGATTTGCACTATTTC‑CC‑T‑AAAA‑ACGCGGCTCCATTTAGGGAGCCGC‑T‑TTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCg < 4:200792/138‑1 (MQ=255) | TCTGCAACGTTTGCTGATTTATATTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTC‑CC‑G‑AAAA‑A‑GCGGCTCCATTTAGGGAGCCGC‑T‑TTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCG > BS168_reference_genome/608027‑608295 |
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: ‑ |