Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 608,164 | C→A | 17.4% | 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 | 17.4% | 98.4 / 4.3 | 46 | intergenic (‑89/+82) | ydgF/dinB | putative amino acid permease/nuclease inhibitor |
Reads supporting (aligned to +/- strand): ref base C (21/15); new base A (2/7); total (24/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.06e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCTGCAACGTTTGCTGATTTATATTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTT‑CC‑C‑G‑AAAA‑A‑GCGGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGT > BS168_reference_genome/608027‑608293 | tctGCAACGTTTGCTGATTTATATTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATATATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTAc > 2:216181/1‑138 (MQ=255) tctGCAACGTTTGCTGATTTATATTGTCTTTCGTCATGTCGTCTCTAACTTTATGTGTTCCTCCATTGAAAATTCTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTAc > 7:218039/1‑138 (MQ=255) ttGCTGATTTATATTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTAc < 6:2068/128‑1 (MQ=255) ttGCTGATTTATATTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTAc > 5:2068/1‑128 (MQ=255) ttGCTGATTTATATTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGAt > 1:36507/1‑139 (MQ=255) gCTGATTTATATTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACCCTAAAATTTCTATCTTTTTTACTAAAAAAGGAtt > 1:211511/1‑138 (MQ=255) gCTGATTTATATTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGAtt > 3:280045/1‑138 (MQ=255) tttATATTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGaa < 1:2123/139‑1 (MQ=255) tttATATTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGaa < 5:243259/139‑1 (MQ=255) tataTTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTACATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATAATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGaaa < 7:346126/138‑1 (MQ=255) atatTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGaaaa < 1:38737/138‑1 (MQ=255) tGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAAtt > 6:341335/1‑139 (MQ=255) tCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTTTTGTGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTTCTAAAAATGGATTGATGAAAAGAAtttt > 3:264805/1‑139 (MQ=255) tCTTTCGTCATGTCGTCTGTCACTTTATGTGTTACTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAAtttt < 3:19740/139‑1 (MQ=255) tCTTTCGTCATGTCGTCTGTCACTATATGTGATACTACATTTAAAATTGTTTAATTATTCAGCAATATATTTTTAATTAATATTATAACACTAAAATTTATATCTTTTTTAATAAAAAAGTATTGATGAAAAGAAttt > 6:165441/1‑138 (MQ=255) tCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTa < 4:7569/136‑1 (MQ=255) tCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAAt < 4:280045/138‑1 (MQ=255) cGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTaaaa > 6:337547/1‑112 (MQ=255) cGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTaaaa < 5:337547/112‑1 (MQ=255) cGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTATATCTTTTTTACTAAAAAAGTATTGATGAAAAGAATTTTAAtgt > 7:87919/1‑139 (MQ=255) aTGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGAt < 6:167662/138‑1 (MQ=255) aTGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGAt < 2:233253/138‑1 (MQ=255) tGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTTTAACACTAAAATTTCTATCTTTTTTAATAAAAAAGGATTGATGAAAAGAATTTTAATGTGAttt > 4:190117/1‑139 (MQ=255) tATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATACCACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTT‑CC‑C‑G‑a > 3:318282/1‑138 (MQ=255) tATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTT‑CC‑C‑G‑a > 4:280649/1‑138 (MQ=255) tgtTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTT‑CC‑C‑G‑aaaa > 2:117932/1‑137 (MQ=255) ttCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTT‑CC‑C‑G‑aaaa‑aagc > 3:177776/1‑136 (MQ=255) ttCCTCCATTGAAAATTGTTTATTGAGACAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATATTTTTTACTAAAAAATGATTGATGAAATGAATTATAATGTGATTTGGACTATTT‑CC‑C‑T‑aaaa‑aagc > 1:305038/1‑136 (MQ=255) cctccATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTT‑CC‑C‑G‑aaaa‑aagcgg > 7:364713/1‑134 (MQ=255) cctccATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTT‑CA‑C‑G‑aaaa‑aagcgg > 3:89457/1‑134 (MQ=255) ctccTTTTAATATTTTTTATTGAGTCATCAATATATTTTTAATGATATTTAAAACAATAAAATTTCTATATTTTTTACTAAAAAAGGATTTATAAAAATAATTTTAATGTTATTTGGTATATTTCCC‑C‑A‑AAAA‑A‑GCGGc < 5:165441/139‑1 (MQ=255) ctccATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAAAACTAAAATTTCTATCTTTTTTAGTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTT‑CC‑C‑G‑aaaa‑aagcggg > 7:339956/1‑133 (MQ=255) tATTGAGTCACCAATCTATTTTTAATGATTATTATAACAATAAAATTTCTATATTTTTTACTAAAAAAGGATTGATTAAAATAATTTTAATGTTATTTGGACTATTT‑CC‑CTA‑AAAA‑A‑GCGGCTCCATTTAGGGAGCCg < 8:282002/138‑1 (MQ=255) gTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTT‑CC‑C‑GAAAAA‑A‑GCGGCTCCATTTAGGGAGCCGCtttttt > 6:64974/1‑139 (MQ=255) aaTCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTT‑CC‑C‑GAAAAA‑A‑GCGGCTCCATTTAGGGAGCCGCgttttttttgtt > 7:10175/1‑133 (MQ=255) aTTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTT‑CC‑C‑GAAAAA‑A‑GCGGCTCCATTTAGGGCGCCGCGTTTTTTTTCtttttta > 2:8894/1‑137 (MQ=255) taatttttataaTAAAAATAAAATTTCATTTTTTTTTAATAAAAATGGATTGATTAAAATAATTTTATTTTTATTTGTACTATTTCCC‑C‑A‑AAAA‑A‑GCGGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGttt < 2:36507/132‑1 (MQ=255) attattATAACACTAAAATTTCTATCTTTTTTACAAAAAAAGTATTCATAAAACCAATTTTAATGTGATTTGGACTATTG‑CC‑CTA‑AAAA‑A‑GCGGCGCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTc < 7:56066/138‑1 (MQ=255) ttatttTAAAACTAAAATTTCTATATTTTTAAATAAAAATGTATTGATGAAAATAATTTTAATGTGATTTTTACTATTT‑CCGT‑A‑AAAA‑A‑GCGGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCa < 8:218039/138‑1 (MQ=255) ttattaTAACACTAAATTTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTT‑CC‑C‑GAAAAA‑A‑GCGGCTACATTTAGGGAGACGCGTTTTTTTTGTTTTTTAGtttttttaa > 3:264067/1‑134 (MQ=255) ttaaAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAATAATTTTAATGTGATTTGGACTATTT‑CC‑C‑G‑AAAA‑ACGCGGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATaaa < 1:216181/139‑1 (MQ=255) taaaaCAATAATATTACTATTTTTTTTAAAAAAATATTATTTTTTAAAATATTTTTACTTGTTTTTTAATTATTT‑CC‑C‑C‑AAAA‑ACGCGGCTCCATTTAGGTATCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATaaa < 8:364713/135‑1 (MQ=255) cacTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTT‑CC‑C‑GAAAAA‑A‑GCGGCTCCATTTAGGGAGCCGCGTTTTTTTTGTTTTTTATTTTTTTTTATAAAAAAGTa > 2:234380/1‑139 (MQ=255) cacTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTT‑CC‑C‑GAAAAA‑A‑GCGGCTCCATTTAGGGAGCCGCGTTTTTTTTGTTTTTTATTTTTTTTGATAAAAAGGtt > 6:288535/1‑138 (MQ=255) aatAAAATTTATATTTTTTTTAAAAAAAAAGTATTTATAAAAATAATTTTATTGTGATTTGTACTATTC‑CC‑C‑C‑AAAA‑ACGCGGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTag < 5:288535/137‑1 (MQ=255) aaaaaatTTTATATTTTTTTAAAAAAAAATGTTTTTATTAAAATATTTTTAATGTTTTTTGCAATATTC‑CC‑C‑A‑AAAA‑ACGCGGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTag < 1:234380/136‑1 (MQ=255) aaaaTTTATATCTTTTTAAAAAAAAAAGGATTGATCAAAAGAATTTTAATGTGATTTGGACTATTG‑CC‑C‑C‑AAAACC‑GCGGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGc < 3:280649/139‑1 (MQ=255) tcaatatttttaACAAAAAACGGATAGAAGAAAAGAATTTAAATGTCATTTCGACTAATA‑CC‑C‑C‑AAAACA‑GCGGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTa < 4:258118/132‑1 (MQ=255) cTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTT‑CC‑A‑GAAAAA‑A‑GCGGCTCCATTTAGGGATCCGCTTTTTTATGCTTTTTTACTTTATGTAAAAAAAACGTAAATCCCGaaa > 8:297059/1‑138 (MQ=255) ttttttAATAAAAATGTTTTGATTAAAATAATTTTATTGTTATTTTGATTATTA‑CC‑C‑T‑AAAA‑ACGCGGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGt < 6:175513/138‑1 (MQ=255) | TCTGCAACGTTTGCTGATTTATATTGTCTTTCGTCATGTCGTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTT‑CC‑C‑G‑AAAA‑A‑GCGGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGT > BS168_reference_genome/608027‑608293 |
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: ‑ |