Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP006881 | 3,725,205:1 | +T | 100% | intergenic (‑37/‑492) | U712_19205 ← / → U712_19210 | Uncharacterized protein ywcI/Minor extracellular protease vpr |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP006881 | 3,725,205 | 1 | . | T | 84.8% | 60.0 / 11.9 | 33 | intergenic (‑37/‑492) | U712_19205/U712_19210 | Uncharacterized protein ywcI/Minor extracellular protease vpr |
Reads supporting (aligned to +/- strand): ref base . (0/5); new base T (18/7); total (21/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.56e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. |
TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑C‑T‑TTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCACGGGA‑AAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATAT > CP006881/3725068‑3725337 | tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACCCCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑C‑t > 7:188762/1‑139 (MQ=255) tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑C‑t > 3:293952/1‑139 (MQ=255) gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGC‑C‑t‑t > 1:258828/1‑139 (MQ=255) gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑C‑t‑t > 4:258373/1‑139 (MQ=255) gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGGAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGGGGGTTTTTGGCAGGCGGG‑C‑t‑t > 4:294838/1‑139 (MQ=255) aaTAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTTGCAGGCAGC‑C‑t‑ttttagc > 7:323965/1‑135 (MQ=255) aaTAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑C‑t‑ttttagc > 1:27181/1‑135 (MQ=255) accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCCCCCCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACCTCCCTCCGCTGCTTTTTGGCAGGCCGC‑CTT‑TTTA‑GCAGccc > 7:207266/4‑139 (MQ=255) accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATCCAAGCAAACGTTTCATAGAAAACCATCCCTCCGCTTGTTTTTGGCAGGCAGC‑CTT‑TTTa‑gcagc > 6:127958/4‑137 (MQ=255) accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑TTTA‑GCAGccc > 3:270498/4‑139 (MQ=255) cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑TTTA‑GCAGCCCGt > 3:286278/2‑139 (MQ=255) ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATC‑AAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACTTCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑TTTA‑GCAGCCCGTTTtc > 8:96541/1‑139 (MQ=255) cAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑C‑TATTTA‑GCAGCCCGTTTtct < 2:4181/139‑1 (MQ=255) aCCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑TTTA‑GCCGCCCGTTTTCTCAGCCGCAGCCCGCAAGA‑AAAGACgg > 6:74602/1‑139 (MQ=255) cccAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑TTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGA‑AAAGACGGc > 7:157968/1‑139 (MQ=255) ttCATCAAAAAGACCATCCACCCCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCCGC‑CGT‑TTTA‑GCCGCCCGGTTTCTCAGCCGCCGCCCGCAAGC‑AAAGACGGCCGAtt > 6:72144/1‑139 (MQ=255) ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGC‑CTT‑TTTA‑TCAGCCCGTTTTCTCAGCCCCACCCCGCAAGA‑AAAGACGGCCGAtt > 1:245536/1‑139 (MQ=255) ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGC‑CTT‑TTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGA‑AAAGACGGCCGAtt > 8:141976/1‑139 (MQ=255) caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑TTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGA‑AAAGACGGCCGAtttcttt < 3:258373/139‑1 (MQ=255) caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑TTTA‑GCAGACCGTTTTCTCAGCCGCCGCCCCCAAGA‑AAAGACGGCCGAtttcgtt > 5:244875/1‑139 (MQ=255) caAAAAGAACATCCACCCCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTCGCAGGCAGC‑CGT‑TTTA‑GCAGCCCGTTTTCTCAGCCGCAGACCGCAAGA‑AAATACGACCGAtatattt > 7:42767/1‑139 (MQ=255) acCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑TTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGA‑AAAGACGGCCGATTTCTTTTCTCGCAAACTAGCgg > 2:135/1‑139 (MQ=255) ccGAAAAGATACACGCAAACTTTTAATACAAAAACAGCCCTCCGTTTTTTTTTGGCCGGCCTC‑CTT‑TTTA‑CCCGCCCGTTTTCCCCGCCGCAGCCCCCAACA‑AACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGct < 6:244875/139‑1 (MQ=255) tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CGT‑TTTA‑GCAGCCCGTTTTCTCATCCGCAGCCCGCAAGA‑AAAGACGGCCGATTTGTTTTTTAGCCTACTATCggggc > 3:8176/1‑134 (MQ=255) aGATAAAACCAAACTTTACATACAAAAACTTCCCTCGGTTTTTTTTTGGCTGGCTGCTT‑T‑TTTA‑GCAGCCCGTTTTCCCCGCCGCAGCCCGAAATA‑AACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa < 4:293952/139‑1 (MQ=255) aTACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑TTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGA‑AAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaa > 3:136817/1‑138 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTGTTTTTGGCAGGCCGC‑CTT‑TTTA‑GCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGA‑AAAGACGGCCGAGTTCGTTTGTCGCCAACTAGCGGCTCGAGGaaaa > 2:134569/1‑138 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑TTTA‑TCAGCCCGTTTTCTCCGCCGCAGCCCGCAAGA‑AAAGACGGCCGATTTCTTTTCTCGCAAACTATCGGCTCTAGCAAAAc > 1:42814/1‑139 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑TTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGA‑AAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa > 2:176775/1‑138 (MQ=255) gaaaaCTTTTAATAAAAAAACTCCCCACGGATTCTTTTTGGCCGGCCGCTT‑T‑TTTC‑CCCGCCGGTTTGCCCCGCCCCCGCCGCAAAGA‑AACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt < 4:320349/137‑1 (MQ=255) ccAAACGTGCCATCCAAAACCATCCCTCCGCTGCTTTTTGGCCGGCAGC‑CTT‑TTTA‑GCCGCCGGTGTTCTCCGCCGCCCCCCCCAACA‑AACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt < 8:188762/138‑1 (MQ=255) cAAACGTTTCAAACAAAAACACCCCACCGCGTCTTTTTGGCAGGCACC‑CTT‑TTTA‑CCCCCCCGGTTTCCCCCCCCCCCCCCCCAACACAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt < 7:79951/139‑1 (MQ=255) aaccttttCATATAAAAACCTCGCGCCGTTTTTTTTGGGCCGGCCTCAT‑T‑TTTC‑CCAGCCCTTTTTCCCCCCCGCCCCCCGAAAAA‑AACCACGGCCGTTTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa < 8:136109/134‑1 (MQ=255) ccGTTTCATACAAAAACATCCCACCGTTTCTTTTTGGCTGGCAGC‑CTT‑TTTA‑GCAGCCCGTTTTCTCACCCGCCGCCCGCAATA‑AAAGACGGCCGATTTCTTTTCTCGAAAACTAGCGGCTCTAGCAAAACTCATGAct < 7:21023/138‑1 (MQ=255) aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑TTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGA‑AAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct > 1:279150/1‑139 (MQ=255) aCGTTTCATACAAAAACACCCCTCCGCTTCTTTTGGGCTGGCCCC‑A‑T‑TTTACGCCGCCGTTTTTCTCCGCCGCCGCCCCCAAAA‑AACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct < 7:141976/139‑1 (MQ=255) gAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑TTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGA‑AAAGACGGCCGATTTCTTTTCTCGCa < 2:223499/100‑1 (MQ=255) gAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑TTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGA‑AAAGACGGCCGATTTCTTTTCTCGCa > 1:223499/1‑100 (MQ=255) caccccTCCGCTTCTTTTTGGCCGGCCGC‑CTT‑TTTA‑GCAGCCGGTTTTCTCCCCCCCAGCCCCCAATA‑AACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTgg < 1:135/136‑1 (MQ=255) gcagC‑CTT‑TTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGA‑AAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAatat > 1:247157/1‑139 (MQ=255) | TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑C‑T‑TTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCACGGGA‑AAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATAT > CP006881/3725068‑3725337 |
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: ‑ |