Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 408,486 | A→C | 15.6% | R134R (CGT→CGG) | yczE ← | integral inner membrane protein regulating antibiotic production |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 408,486 | 0 | A | C | 15.6% | 124.9 / 2.4 | 63 | R134R (CGT→CGG) | yczE | integral inner membrane protein regulating antibiotic production |
Reads supporting (aligned to +/- strand): ref base A (33/19); new base C (2/9); total (35/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.13e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTGCCAAAACCGATCGGTCCGCCCATGCCCCATGCCGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACGATGATCGAGCCTGTGTAGGTCGAGAGGGCAGGCAGAATGAAATTGAAAAAATCTATAAA > BS168_reference_genome/408364‑408623 | gggCCAAAACCGCCCGGGCCGCCCAAGCCCCATCCCGCAAAAAAATTGGTTAATTCCATGCCGTGCCCCACCAAACCCACACCCCAGCGGGTTTTTTCTGTAATCAGCAACAT‑CAGTCACACACGCGGCCCCGCGcc < 4:53749/135‑1 (MQ=255) gCCAAAACCGATCGGTCCGCCCATGCCCCATGCCGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAgg > 4:269311/1‑139 (MQ=255) gCCAAAACCGATCGGTCCGCCCATGCCCCATGCCGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAgg > 5:63058/1‑139 (MQ=255) ccAAAACCGATCGGTCCGCCCATGCCCCATGCCGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGc < 3:415849/139‑1 (MQ=255) aCCGATCGGTCCGCCCATGCCCCATGCCGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCctgct < 3:77412/139‑1 (MQ=255) cGATCGGCCCGCCCATGCCCCATGCCCCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCCCGCGGCCCCGCGCCAAGGCCTGCTGa < 5:131469/139‑1 (MQ=255) cttgccccaggccCCAAAAAAAAGGTTAAATTCCAGCCCGTTCCGCACCATTTGCACATTCCAGCGGTTTTTTTTTTTAATCTGCATCAT‑CTGTCAAACCCGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCg < 6:21860/137‑1 (MQ=255) gccccatgccGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAGTCAGCATCTT‑CAGTGAATCACGCGGCCCCGCGCCACGGCCTGCTGATACGTGAGCAccgccg > 8:154532/1‑139 (MQ=255) ccccatgccGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCCGCCCCCCGCCAAGGCCTGCTGATACAAAAACAccacccc > 7:154834/1‑134 (MQ=255) atgccGCAAACAAAATTGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGt > 3:394499/1‑138 (MQ=255) atgccGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGCCGCCGt > 5:329981/1‑138 (MQ=255) atgccGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTATGTAATAAGCATCAT‑AAGATAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACAccgccgccga > 3:6397/1‑137 (MQ=255) tgccGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGccc > 7:191242/1‑99 (MQ=255) tgccGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGccc < 8:191242/99‑1 (MQ=255) tgccGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGcc > 4:399023/1‑111 (MQ=255) tgccGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGcc < 3:399023/111‑1 (MQ=255) tgccGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTaa > 4:173432/1‑139 (MQ=255) tgccGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTaa > 3:195305/1‑139 (MQ=255) tgccGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTaa > 7:307557/1‑139 (MQ=255) gccGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCGGGTTTTTTCTGTAATCAGAATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTaa < 6:382779/138‑1 (MQ=255) ccGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAAc < 7:224403/138‑1 (MQ=255) caaacaaaATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTaa > 7:31459/1‑134 (MQ=255) caaacaaaATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTaa < 8:31459/134‑1 (MQ=255) caaaATGGTTAATTCCATGCCGTCCCGCACCCATTGCACATTCCCGCCGGTTTTTTCTTTAATCAGAAGCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGc < 4:421425/139‑1 (MQ=255) aattgtttAATTCTATGCCGTTCCGCCCCATTTCCAAATTCCTGCGTTTTTTTTGTTTAATCAGCTACAT‑CTGTAAAACACGCGGCCCCGCGCCAAGGCCTGCTGATAAATAAACACCGACGCCGTAACCAATCAGcac < 2:103560/133‑1 (MQ=255) aactGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCGGGTGTTTTCTGAAATCAGCATCAG‑CAGTGAACCCCGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGcac < 5:18863/136‑1 (MQ=255) aaaTGGTTATTTCCAGGCCGTGCCTCACCCATTGCACATTCCCGCCTGGTTTTGCTGTAACCAGCATCAT‑CAGTCCAGCCCGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGcac < 6:12770/139‑1 (MQ=255) aaaTGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGcac < 3:269311/139‑1 (MQ=255) aaaTGGTTAATCCCAGGCCGTCCCCCCCCAATGCCACAATCCCGCGGTTTTTTTCGTCAATCACCAACAA‑CCGGCAATCCCGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGcac < 4:236519/139‑1 (MQ=255) aaTTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGcc > 6:27712/1‑137 (MQ=255) aaTTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGcc > 3:226783/1‑137 (MQ=255) ttCCATGCCGTGCCGCCCCCTTTCCAATTTCCCGCGGTTTTTTTGTTTAATC‑CCAACATACTGTTAAGCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTaga < 4:226783/139‑1 (MQ=255) cATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTACCCAATCAGCACCACGCCTaga > 5:405875/1‑136 (MQ=255) tGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTagagag > 5:70452/1‑137 (MQ=255) tGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGGAACCAATCAGCACCACGCCTagagag > 1:345685/1‑137 (MQ=255) ttCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTTATCAGCATCAT‑CAGTGAATCCCGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGACAAGCCga > 7:165434/1‑139 (MQ=255) ttCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACga > 3:323737/1‑139 (MQ=255) ttCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACga > 7:128336/1‑139 (MQ=255) ttCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACga > 3:379031/1‑139 (MQ=255) ttCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATAAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTCATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGAc > 1:217445/1‑137 (MQ=255) gCACCATTTGAACATGCATGCCGGTTTTTTCTTTATTACGAATCAT‑CTGTTATTCCGGCGCCCCCGCGCCTAGGCCGGCTGTTACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACgttgat < 6:326962/139‑1 (MQ=255) ccATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGAAGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACGATGATCGAg > 2:52098/1‑139 (MQ=255) caaTTGCCAATCCCAGCCGGTTTTTTCTGAAATCCGAATCAT‑CCGTGAACCCCGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACGATGATCGAg < 4:76111/137‑1 (MQ=255) aTTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCCCGCGGCCCCGCGCCAAGGCCTGCTGCTACATAAACACCGCCGCCGTACCCAATCAGCCCCACGCCTAGAGAGCAGACGATGATCGCGcc > 5:62750/1‑139 (MQ=255) tGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACgatgat > 7:393347/1‑131 (MQ=255) tGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACgatgat < 8:393347/131‑1 (MQ=255) ctttCCAGCGGGTTTTTTCTGTAATCAGCATAAT‑CTGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACGATGATCGAGCCTGTGTa < 8:286168/137‑1 (MQ=255) ttCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGcac > 7:312951/1‑101 (MQ=255) ttCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGcac < 8:312951/101‑1 (MQ=255) ttCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACGATGATCGAGCCTGTGTTgg > 4:384032/1‑139 (MQ=255) ttCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCGGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACGATGATCGAGCCTGTGTAgg > 3:29229/1‑139 (MQ=255) ttCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCAGCTGATCCATATACACCGACGCCTTAACCAATCAGCACCACGCCTAGAGAGTAGACGATGATCGAGCCTGTGTAgg > 8:81072/1‑139 (MQ=255) tCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACGATGATCGAGCCTGTGTAGGt > 3:76111/1‑139 (MQ=255) aGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGCCGATGATCGAGCCTGTGTAGGTCga > 4:9450/1‑139 (MQ=255) gCCGGTTTTTTCGTAAAACAGCCCCAT‑CCGGGACCCCCGCGGCCCCGCGCCAAGGCCTGCTGATACATAACCACCGACGCCGTAACCAATCAGCACCCCGCCTAGAGAGAAGACGATGATCGAGCCTGTGTAGGTCgag < 2:345685/139‑1 (MQ=255) cat‑caGTGACTCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGCCGATGATCGAGCCTGTGTAGGTCGAGAGGGCAGGCAGAATGAAATTGaaa < 6:80830/139‑1 (MQ=255) cat‑caGTGAATCCCGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACGATGATCGAGCCTGTGTAGGTCGAGAGGGCAGGCAGAATGAAATTGaaa < 8:128336/139‑1 (MQ=255) cat‑caGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACGATGATCGAGCCTGTGTAGGTCGAGAGGGCAGGCAGAATGAAATTGaa > 1:251248/1‑138 (MQ=255) t‑caGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACGATGATCGAGCCTGTGTAGGTCGAGAGGGCAGGCAGAATGAAATTGaaaa > 7:303722/1‑138 (MQ=255) aaTCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACGATGATCGAGCCTGTGTAGGTCGAGAGGGCAGGCAGAATGAAATTGAAAAAATCta > 2:309094/1‑138 (MQ=255) cACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACGATGATCGAGCCTGTGTAGGTCGAGAGGGCAGGCAGAATGAAATTGAAAAAATCTATaaa < 4:195305/139‑1 (MQ=255) cACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACGATGATCGAGCCTGTGTAGGTCGAGAGGGCAGGCAGAATGAAATTGAAAAAATCTATaaa < 3:173432/139‑1 (MQ=255) cACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACGATGATCGAGCCTGTGTAGGTCGAGAGGGCAGGCAGAATGAAATTGAAAAAATCTATaaa < 1:52098/139‑1 (MQ=255) | GTGCCAAAACCGATCGGTCCGCCCATGCCCCATGCCGCAAACAAAATGGTTAATTCCATGCCGTTCCGCACCCATTGCACATTCCAGCCGGTTTTTTCTGTAATCAGCATCAT‑CAGTGAATCACGCGGCCCCGCGCCAAGGCCTGCTGATACATAAACACCGACGCCGTAACCAATCAGCACCACGCCTAGAGAGAAGACGATGATCGAGCCTGTGTAGGTCGAGAGGGCAGGCAGAATGAAATTGAAAAAATCTATAAA > BS168_reference_genome/408364‑408623 |
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: ‑ |