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 | 81.4% | 85.9 / 22.9 | 43 | intergenic (‑37/‑492) | U712_19205/U712_19210 | Uncharacterized protein ywcI/Minor extracellular protease vpr |
Reads supporting (aligned to +/- strand): ref base . (1/7); new base T (20/9); total (27/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.18e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-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. |
GAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑C‑TTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGT > CP006881/3725069‑3725329 | gAGCCGAATAGACCTCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCCCACCCGGAAAGAGACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCGTTTTGGCAGGCAG‑C‑C‑tt > 5:190117/1‑139 (MQ=255) aaTAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CATTTT‑a‑g > 7:97905/1‑138 (MQ=255) aaTAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑C‑tttttagc < 7:100833/139‑5 (MQ=255) accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGccc > 2:101754/4‑139 (MQ=255) cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGt > 7:231783/2‑139 (MQ=255) cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGt > 3:249412/2‑139 (MQ=255) cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGGAAAGATACAAGCAAACGTTTCATAGACAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CGATTT‑A‑GCAGCCCGt > 2:73688/2‑139 (MQ=255) gCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGtt < 3:6725/139‑1 (MQ=255) ggCAGCCGCAGTTTTTCTGTCTGGCGTCACCCAGTTCATCAAAAAGCCCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAG‑C‑CGTTTT‑T‑GCAGCCCGtttt > 5:217930/1‑139 (MQ=255) ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGtttt < 4:337676/139‑1 (MQ=255) ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGttt > 4:157955/1‑138 (MQ=255) ccGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTGTTTTTGTCAGGCAG‑C‑CGTTTT‑A‑TCAGCCCGTTTTCTCAg > 5:293024/1‑139 (MQ=255) cAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCg < 4:249412/139‑1 (MQ=255) cGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGAAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGcccgcaaga > 2:362729/1‑132 (MQ=255) caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGAtttcttt > 4:36961/1‑139 (MQ=255) caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATACCTCCGCTTCTTTTTGGCAGGCAG‑C‑CATTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGatttattt > 8:27984/1‑139 (MQ=255) aGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑a‑gcag > 8:213412/1‑87 (MQ=255) aGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑a‑gcag < 7:213412/87‑1 (MQ=255) aGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCg > 4:77854/1‑139 (MQ=255) aCCATCCACACCCGAACAGATAAACGCAAACTTTTCAAACAAAACCAGCCCTCCGCTTCTTTTTGGCAGGCAG‑CGT‑TTTT‑A‑GCCGCCCGTTTCCGCCGCCGCCGCCCGAAAGAAACGACGGCCGATTTCTTTTCTCGCa < 7:243976/139‑1 (MQ=255) cATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑C‑tttt > 8:209247/1‑77 (MQ=255) cATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑C‑tttt < 7:209247/77‑1 (MQ=255) tCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTTGCAGGCAG‑C‑CGTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAACAAAAGACGGACGAATTCTTATCTCGCaa > 3:33650/1‑136 (MQ=255) tAAAGATACAAGCAAACGTTTCATAGAACAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGGCCGTTTTCTCCGCCGCAGCCCGCCAGAAAAGACGGCCGATTGCTTTTCTCGCAAACTAGCGGcttgt > 2:268618/1‑136 (MQ=255) tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCCGCCGCAGCCCGCAAGATAAGACGGCCGAGTTCTTTTCTGGCAACCTAGCGGctctc > 2:138322/1‑138 (MQ=255) tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTa > 4:330537/1‑139 (MQ=255) aGATACAAGCAATTGTTTCATATAAAAAAATCCCTCCTGTGTTTTTTGTGAGTCAGAG‑T‑TTTT‑C‑GCTGCCTGGTTTCTCCCCCGCCCCCCCAAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa < 6:254091/139‑1 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc > 2:12730/1‑139 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGCTTTCTTTTCTCGCAACATAGCGGCTCTAGCaaaa > 5:286888/1‑138 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCTCAAACTAGCGGCTCTAGCaaaa > 8:361302/1‑138 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAAATAGCGGCTCTAGCaaaa > 5:254091/1‑138 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTATTTTATCGCAAACTAGCCGCGCTAGAAAAAc > 1:187479/1‑139 (MQ=255) gCACACGTTTCATACAAAAACATCCCTCCGCTTGTTTTTGGCTGGCAG‑C‑CTTTTT‑C‑GCAGCCGGTTTTCTCCGCCGCCGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt < 6:286888/139‑1 (MQ=255) gCAACCGTTGCAAACAACAACATCCCGCCGCTGATTTTGGGCAGGCCG‑C‑CTTTTT‑C‑GCCGCCCGTTTTCTCCGCCGCAGCCCGCAAGAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt < 7:72748/139‑1 (MQ=255) caaaccGTTTCATAGAAAAACATCCATCCGCTTCTTTTTGCCAGGAAG‑CAT‑TTTT‑A‑GCACCCCTTTTTCTCAGCCGCCCCCCCCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt < 7:361302/134‑1 (MQ=255) ccGTTCCATAAAAAACCACCCAGCCGCTTCTTTTGGCCTGCCCG‑C‑CTTTTT‑A‑GCCGCCCGTTTTCCCCGCCGCAGCCCCCAACAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct < 1:64467/138‑1 (MQ=255) atttttCATATAAAAACATCGCTGCGCTTTTTTTTGCCTGTCCG‑CAT‑TTTT‑A‑GCAGCCGGTTTTCGCAGCCGCACCCCCAAAAAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct < 8:69013/136‑1 (MQ=255) aCGTTTCATAGAAAAACATCCCTCCGGTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct < 3:157955/139‑1 (MQ=255) aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑CAT‑TTTT‑A‑GCAGCCCGTTTTCTAAGCCGCAGCCCGCAAGAAAAGACGGACGAATTATTTTATAGCAAACTAGCGGCTCTAGCAAAACTCATGAct > 5:131165/1‑139 (MQ=255) aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAtcccg > 8:256491/1‑135 (MQ=255) aCGTTTAATACAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑CAT‑TTTT‑A‑GCAGCCCGTTTTCCCCGCCGCAGCCCGCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct < 3:330537/139‑1 (MQ=255) aCGTTAAATATAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCG‑C‑C‑TTTTAA‑GCACCCCGTTTTCTCCGCCCCACCCCGCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct < 3:77854/139‑1 (MQ=255) gTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCCCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGActct < 8:231783/139‑1 (MQ=255) aaaCATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGAAGGCCGATTTAGTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAATCTATGATAGACCg > 4:220986/1‑139 (MQ=255) cttcttTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTGTCGCAAACTAGCGGCTCGAACAAAACTCATTACTCTATGATATAACGTGGAAAGAGAAtt > 3:33326/1‑139 (MQ=255) cttTTTGGCAGGCAG‑C‑CTTTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCtt > 6:178733/1‑139 (MQ=255) ttGGCCGGCAG‑C‑C‑TTTT‑ATCCGGCCCTTTTGCTCAGCCGCAGCCCCCAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGt < 7:27984/138‑1 (MQ=255) | GAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑C‑TTTT‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGT > CP006881/3725069‑3725329 |
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: ‑ |