Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 1,023,150 | C→A | 20.8% | intergenic (+50/‑200) | yhdF → / → bcaP | putative NAD(P)‑dependent dehydrogenase/branched‑chain amino acid transporter or methylthioribose‑degradation metabolite exporter |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 1,023,150 | 0 | C | A | 20.8% | 89.8 / 2.2 | 39 | intergenic (+50/‑200) | yhdF/bcaP | putative NAD(P)‑dependent dehydrogenase/branched‑chain amino acid transporter or methylthioribose‑degradation metabolite exporter |
Reads supporting (aligned to +/- strand): ref base C (9/20); new base A (4/3); total (16/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.25e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATG > BS168_reference_genome/1023026‑1023277 | ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTTAACGTTATCAAAGAAAAACAACCACAACAGCGGGTTTTTTTTGTGGT‑GATAAAATTAAAAAATAAAAAAtt > 5:60752/1‑139 (MQ=255) ctgGCGGCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAAttt < 6:175793/139‑1 (MQ=255) gCGTCTGACGAATATTACTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTATTTTAAAGAAATAAAACAAAAAACACCACAAGAGCGGGTTTTTTTTGTGGTGGA‑AAAATTTAAAAAAAAAAAAatttg > 1:134950/1‑137 (MQ=255) aTCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTaa < 1:37031/139‑1 (MQ=255) ggCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCAttt > 1:116673/1‑139 (MQ=255) cAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAACGTAATCAACGAAAAAAAAGCTCAAGAGCTGGTTTTTTTTGTGGT‑GA‑AAAAATTCAAAATAAAAAATTTTTCGATTTTTATAAAAATTTTTTTTTAAttttt > 5:154448/1‑139 (MQ=255) aCCATTCATGTGAATGGCGGCCGTGTTTTTTGAAACGAAACAAAGAAAAACCAACCCAAACGCAGGTTTTTTTTGTGGT‑GA‑TAAAATTCAAAAAAAAAAAATTTTAGATTTTTGTAAAAAATTTTTTTTGAttttttg > 6:122631/1‑136 (MQ=255) tCATGTGAATGGCGGCCGTTTTATTTAAACGTAATAAAAGAAAAACCAACTCAAGAGCTGGGTTTTTGTGTGGG‑GG‑AAAATTTTAAAATAAAAAAATTTTCGATTTTTTTAAAAAATTTTTTTTAATTTTTTGAAtttt > 4:127299/1‑139 (MQ=255) tCATGTGAATGGCGGCCGTGTTTTTTAAAAGTAATAAAACAAAAAAAAACGCAAAAGCGGGGTTTTTTTGTGGT‑GT‑AAAATTTTAAAATAAAAAAATTTTGGATTTTTTGAAAAAATTTTTTTTGATTTTTTGAAtttt > 2:206040/1‑139 (MQ=255) cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt < 1:206040/138‑1 (MQ=255) cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt < 4:104809/138‑1 (MQ=255) cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt < 2:134950/138‑1 (MQ=255) cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt < 2:78977/138‑1 (MQ=255) aTGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat < 8:219874/139‑1 (MQ=255) aTGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat < 3:209577/139‑1 (MQ=255) tgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata < 3:127299/139‑1 (MQ=255) tgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata < 1:174737/139‑1 (MQ=255) tgAATGGCGGCCGTTTTATTTAAACGTAATAAAACAAAAAACAGCGAAAAAGCAGGGTTTTTTTGTGGT‑GT‑AAAAATTAAAAATAAAAAAATTTTCGATTTTTGTATAAAATTTTTTTTTATTTTTTAAAtttttttta > 2:174737/1‑134 (MQ=255) gAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTTGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata < 7:64448/138‑1 (MQ=255) aaTGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatatata < 1:13494/139‑1 (MQ=255) tGGCGGCCGTTTTATTTAAACGTAATAAACGAAAAAACAGCTCAATAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTAAAAATAAAAAATTTTTCGATTTTTCTAACAATTTTTTTTTGATTTTTTGAATTTTTTTTataa > 7:219874/1‑138 (MQ=255) tGGCGGCCGTGTTATTTAAACGTAATAAAACAAAAACCAGCTCAAGAGCTGGTTTTTTTTTTGGT‑GA‑AAAATTTGAAAATTAAAAAATTTTGGATTTTTTTAACAAATTTTTTTTCATTTTTTAAATTTTTTATttaa > 4:209577/1‑138 (MQ=255) ccGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg < 2:116673/139‑1 (MQ=255) ccGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg < 7:130706/139‑1 (MQ=255) tATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGga < 5:122631/138‑1 (MQ=255) aTCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGAGTTTTGTAACAATTTTTTTTTAATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGAtt > 5:32716/1‑139 (MQ=255) gAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt < 1:103636/100‑1 (MQ=255) gAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt > 2:103636/1‑100 (MQ=255) aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc < 3:21973/139‑1 (MQ=255) ggTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTGTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGTATGATTCGCGACCAGTACAAATATAGTAGaa > 4:182858/1‑138 (MQ=255) gTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata > 5:223189/1‑89 (MQ=255) gTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata < 6:223189/89‑1 (MQ=255) ttttttGTGGT‑TA‑AAAATTTCAAAATAAAAAATTTGTCTATTTTTCTAACATTTTTGTTTTAATTTTTGAAATTTTATATATAATATTGAGTAGGAGAAAGGAATTATTCGCGACCAGTACAAATATAGTAGaaaaaa < 8:90253/138‑1 (MQ=255) tttgttttgt‑TA‑AAAATTTAAAAATAAAAAATTTTGCTATTTTTATAAAAATTTTTTTTTATTTTTTGCAATTTTATATATAATATTGAGTAGGACAAAGTAATCATTCGCGACCAGTACAAATATAGTAGaaaaaaa < 1:98551/138‑1 (MQ=255) tgggtt‑gA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTTCAATTTTAAATATAATATTGAGTCGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaat < 1:137843/134‑1 (MQ=255) tggt‑gA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGTAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGtataa > 1:31642/1‑139 (MQ=255) tggt‑gA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaataa > 7:40915/1‑139 (MQ=255) ttttta‑aaaaTTTAAAAAAAAAAAATTTGTCTTTTTTTATAAAATTTTTGTTTTAATTTTTTAAATTTTATATATAATTTTGAGTAGTATAAAGTAATTATCCGCGACCAGTACAAATATAGTAGAAAAAAAAGaata < 6:32716/133‑1 (MQ=255) gtta‑aaaaTTTAAAAATAAAAAATTTTTCGATTTTTCTAAAATTTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGAATCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATg < 1:34170/136‑1 (MQ=255) | CTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGT‑GA‑AAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATG > BS168_reference_genome/1023026‑1023277 |
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: ‑ |