Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,852,080 | T→C | 19.2% | V343A (GTC→GCC) | PP_3399 → | minor curlin subunit CsgB |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,852,080 | 0 | T | C | 19.2% | 34.3 / 4.3 | 21 | V343A (GTC→GCC) | PP_3399 | minor curlin subunit CsgB |
Reads supporting (aligned to +/- strand): ref base T (7/10); new base C (2/2); total (9/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01 |
AAAACCTCAACCTGACCATCAACACCCAAGGCAATGGTAACAGCATCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAGGCAACCAGATGTCGACCAAGCAGGATGATGTCAGCAACGA > NC_002947/3851942‑3852176 | aaaaCCTCAACCTGACCATCAACACCCAAGGCAATGGTAACAGCATCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGAccc > 2:70112/1‑150 (MQ=255) aTCAACACCCAAGGCAATGGTAACAGCATCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCaa < 2:192541/150‑1 (MQ=255) aaGGCAATGGTAACAGCATCACCGCCAGTAACTGCAAAGGCGCCAAGATGGGTGGGCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAAGGGCAACGGCGTCAATCTCACGCAGAAGGGAGCTGACCAACTGGCCACGc < 2:57389/150‑1 (MQ=255) ggCAATGGTAACAGCTTCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGCCAATCTGa > 2:576898/1‑118 (MQ=255) ggCAATGGTAACAGCTTCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGCCAATCTGa < 1:576898/118‑1 (MQ=255) ggCAATGGTAACAGCATCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGCCAATCTGa < 1:576726/118‑1 (MQ=255) ggCAATGGTAACAGCATCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGCCAATCTGa > 2:576726/1‑118 (MQ=255) tAACAGCATCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGAc < 2:710986/144‑1 (MQ=255) tAACAGCATCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGAc > 1:710986/1‑144 (MQ=255) cGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAGGCaa < 2:258657/150‑1 (MQ=255) gCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAGGCAAc < 2:194848/150‑1 (MQ=255) gCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAGGCAAc > 1:194848/1‑150 (MQ=255) tGGAAAGGCGCCAAGATGGGTGGCCAGTTAGGCACCAACCAGACCGCAACCACCAACCAGGCTGGCAATGGCCCCGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCAAGACGGATAAAGGCAACCGGACGTcc > 2:170862/1‑149 (MQ=255) cAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCt < 2:268601/91‑1 (MQ=255) cAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCt > 1:268601/1‑91 (MQ=255) aacaTGGGTGGCCAGTTCAGCACAAACCAAAACGCCACCATCATGCTGGCTGGCAATTGCAACGGCGTCAATCTGTATCAGAAGGGCGCTAACCGACTGGCCACGCTGACCCAGACGGGTAAAAGCAACCAGATGTCGACCAAGCAGgat < 2:144464/147‑1 (MQ=255) cAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAg > 1:587047/1‑112 (MQ=255) cAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAg < 2:587047/112‑1 (MQ=255) gTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAGGCAACCAGATGTCGACCAAGCAGGATGATGTCAGCAACGa < 1:70112/150‑1 (MQ=255) caaccaGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAGGCAACCAGATGTCGa > 2:414304/1‑116 (MQ=255) caaccaGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAGGCAACCAGATGTCGa < 1:414304/116‑1 (MQ=255) | AAAACCTCAACCTGACCATCAACACCCAAGGCAATGGTAACAGCATCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAGGCAACCAGATGTCGACCAAGCAGGATGATGTCAGCAACGA > NC_002947/3851942‑3852176 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |