Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP006881 | 3,725,237 | A→G | 100% | intergenic (‑69/‑460) | 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,237 | 0 | A | G | 89.5% | 43.1 / ‑2.9 | 18 | intergenic (‑69/‑460) | U712_19205/U712_19210 | Uncharacterized protein ywcI/Minor extracellular protease vpr |
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (12/4); minor base C (1/1); total (13/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.90e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑C‑‑TTTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTT > CP006881/3725119‑3725360 | ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CT‑TTTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt > 7:180560/1‑139 (MQ=255) cccGAAAAGAAACCACCACACGTTAAATCCAAACACATCCCTCCCCTGCTTTGTGGCCGCCAGCTT‑‑TTTTA‑GCAGCCCGTTTTCGCAGCCCCAGCCCGCACCAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGc < 2:317527/139‑1 (MQ=255) agcctAAACATACAAGAAAATGTTAAAAACAAAAACATCCCCCCGCTTCTTTTTGGCTGGCAGC‑CT‑TTTTC‑GCCGCCGTTTTGCTCAGCCGCCGCCCGCAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGc < 1:366089/135‑1 (MQ=255) tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CT‑TTTTA‑GCAGCCCGTTTTATCAGCCGCAGCCCCCAAGAAAAGACCGCCGATTTATTTTATCGCAAACTAGCGGcactt > 5:136696/1‑138 (MQ=255) ctaTACAAGCAAAGCTTCAATATAACAACACCCCACGGCTTTTTTTTGCCAGCCCGC‑C‑‑TTTTACGCACCCCTTTTTCCCCGCCCCCCCCCCCACCACACGCCGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa < 4:85857/137‑1 (MQ=255) aGACACAACCAAACGTTTCATACAAAAACACCCCCCCGCTGCTTTTTGGCAGCCAGC‑CT‑TTTTA‑GCAGCCCGTTTCCCCCGCCCCCCCCCGCAAGAACAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa < 7:389399/139‑1 (MQ=255) aCAAGCAAACGTTTCATAGAAACACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CT‑TTTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCGTGTCTCGCCAACTAGCGGCGCGAGCAAAAc > 4:157821/1‑139 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CT‑TTTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACGAGCGGCTCTAGCAAAAc > 7:388881/1‑139 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CT‑TTTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTCCTTTTCTCGCAAACTAGCGGCTCTAGAAAAAc > 3:126720/1‑139 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CT‑TTTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACAACCGGCTCTAGCaaaa > 5:279602/1‑138 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CT‑TTTTA‑GCAGCCCGTTTTCTAAGCCGCCGCCCGCAAAAAAAGACGGCCGATTTAGTTTCTCGCAAACTACCCGCACAAGCaaaa > 5:128703/1‑138 (MQ=255) aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCATTTTGGCAGGCAGC‑CT‑TTTTA‑GCAGCCCGTTTTCTCAGCCGCACCCCGCAAGAAAACACGGCCGATTTATATTCACGCAAACAAGCGGCGCTAGCAAAAccg > 2:366089/1‑137 (MQ=255) aaaaCTTTTCATCTAAAAACATCCCGCGGGTTC‑TTTTGGCACGCAGC‑CCTTTTTC‑GCGGCCGGTTTTCTCACCCGCTCCCCGAAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATg < 7:75327/138‑1 (MQ=255) ggcagC‑CT‑TTTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTATCGCAAACTAGCGGCGCGAGCAAAACTCATGCCGGGATGAAAGCCCGTGGAAAGAGAATTGATTTTCACAata > 2:402160/1‑139 (MQ=255) gcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGATTGAAAAAACCa > 2:300783/1‑138 (MQ=255) gcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCaa > 3:42858/1‑139 (MQ=255) cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTATCGCTAACTAGCGGCTCTAGCAAAACTCCTGACTCTATGTTAGACCGTGGAAAGAGAATTGTTTTTCACAATATGAGTGTAAAAACCaa > 4:321692/1‑138 (MQ=255) cAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGtttt > 3:140606/1‑136 (MQ=255) | TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑C‑‑TTTTA‑GCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTT > CP006881/3725119‑3725360 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |