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.8% | 76.0 / 15.3 | 33 | intergenic (‑37/‑492) | U712_19205/U712_19210 | Uncharacterized protein ywcI/Minor extracellular protease vpr |
Reads supporting (aligned to +/- strand): ref base . (4/2); new base T (20/4); total (27/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.71e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. |
GAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑C‑T‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATAT > CP006881/3725069‑3725337 | gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCCCCCAGTTCATCAAAAAGCCCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGGAGTTTTGGCCGGCAGc‑cgt > 1:181662/1‑137 (MQ=255) gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTTGCCGGAAGC‑C‑t‑t > 1:314649/1‑139 (MQ=255) gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGGTCATCAAATAGCCAATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACACCCCTCCGCTTGTTTTTGGCCGGCAGC‑C‑t‑t > 5:399859/1‑139 (MQ=255) aaTAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑C‑t‑t‑tttagc < 2:103669/139‑5 (MQ=255) accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCAACCACACCCGTAAAGAGACACGCAACCGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGC‑C‑T‑TGTTAGCAGccc > 3:160026/4‑139 (MQ=255) cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGCCCATCCACCCCCGTAAAGCTACAAGCAAACGTTTCATAGAAAAACATCACTCCGCTTCTGTTTGGCAGGCAGC‑CTT‑T‑TTAGGAGCCCgg > 5:230937/2‑138 (MQ=255) cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑C‑TGT‑TTAGCAGCCCGt > 7:42063/2‑139 (MQ=255) ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCAGCCCGtttt > 7:29728/1‑139 (MQ=255) ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CGT‑T‑TTAGCCGCCCGttt > 5:117754/1‑138 (MQ=255) tCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAACCGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGGAGC‑C‑TGT‑TTAGCCGCCCGTTTTCTCCGCCGCCGcccgc > 2:30452/1‑134 (MQ=255) tCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGcccgc < 5:12609/136‑3 (MQ=255) tCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGcccgc > 8:68005/1‑134 (MQ=255) ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACTGCCGAtt > 8:170490/1‑139 (MQ=255) ttCATAAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGATTCATAGAAAAACATCCCACCCCGTCTATTTGGCAGGCAGC‑C‑TGT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGcccagt > 8:400240/1‑135 (MQ=255) tcatcaAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCTGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAttt > 7:160464/1‑139 (MQ=255) aTCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCACC‑CTT‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCAAGCaaa > 7:380341/1‑138 (MQ=255) aTCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCACGCAGC‑CAT‑T‑TTAGCACCCCGGTTTCTCAGCCGCAGCCCCCCAGAAAAGCCCGCCGCATTCTAATCTCGCaaa > 6:242390/1‑138 (MQ=255) tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGcc < 1:364049/92‑1 (MQ=255) tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGcc > 2:364049/1‑92 (MQ=255) aaGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGAAGCCCGTTTTCTCAGCCGCAGCCCGCAACAAAAGACGGCCGATTTGTTTTCACGCAAACTAGCGGCTCTAg > 3:175736/1‑138 (MQ=255) aGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCACCCCGTTTTCTCAGCCGCCCCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa < 3:337730/139‑1 (MQ=255) aTACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCCGCCCGTTTTCTCAGCCGCCGCCCGCCAGAAAAGACGGCCGCTTTCTTTTCTCGCAACCTAGCGGCGCTATCaa > 8:357591/1‑138 (MQ=255) aTACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaa > 4:361695/1‑138 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa > 8:12548/1‑138 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa > 1:277184/1‑138 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc > 8:342772/1‑139 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc > 3:288943/1‑139 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTAGTTTTGGCAGGCAGC‑CAT‑T‑TTAGCAGCCCGGTTTATCAGCCGCAGCCCCCAAGAAAAGACGGCCGATGTATTTTATCGCAAAATAGCGCCACTAGCaaaa > 4:119252/1‑138 (MQ=255) accAACCGTTTCATATAAAAACCTCCCTCCGTTTCTTTTTGGCAGCCCGC‑CTT‑T‑TTAGCCGCCCGTTTCCGCCGCCCCCCCCCTAAATAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCa < 8:29728/137‑1 (MQ=255) ccAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCCGGCAGCTT‑T‑T‑TTAGCAGCCCGTTTTCTCCCCCGCAGCCCGCAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt < 8:14399/138‑1 (MQ=255) aaaCGTTTCATACAAAAACTCCCCTCCGCTTCTTTTTGGCAGGCAGCTT‑T‑T‑TTAGCAGCCCGTTTTCTCAGCGGCAGCCCGAAAGAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa < 3:361695/139‑1 (MQ=255) cATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGATTTGATTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGGCACTaag > 3:275637/1‑136 (MQ=255) cATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGa > 2:14886/1‑139 (MQ=255) cGCTTCTTTTTGGCAGGCAGC‑CTT‑T‑TTAGCAGCCCGTTTTATCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCCAACTACCGGCTCGAGCAAAACAACAGACACTATGATAGCCCCGGGAAAGAGaa > 6:400038/1‑139 (MQ=255) ggcagC‑CTT‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCAAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCATTGTCACAata > 1:285662/1‑139 (MQ=255) gcagC‑CTT‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATTACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAatac > 7:323104/1‑138 (MQ=255) gcagC‑CTT‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACTAGAGGCTCTAGCAACACTCATGACTCGATTCTAGCCCGTGGAAAGCGAATTCTTTGTCACAatat > 6:102206/1‑139 (MQ=255) | GAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑C‑T‑T‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATAT > CP006881/3725069‑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: ‑ |