Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,103,745 | (T)7→8 | intergenic (‑549/‑206) | csgD ← / → csgB | DNA‑binding transcriptional dual regulator CsgD/curlin, minor subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,103,738 | 1 | . | T | 91.4% | 106.5 / 3.2 | 36 | intergenic (‑542/‑213) | csgD/csgB | DNA‑binding transcriptional dual regulator CsgD/curlin, minor subunit |
Reads supporting (aligned to +/- strand): ref base . (2/1); new base T (10/22); total (12/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.66e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑A‑TTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATA > NC_000913/1103614‑1103849 | cAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑A‑tttttttt < 2:27975/133‑2 (MQ=255) cAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑A‑tttttttt < 1:318392/133‑2 (MQ=255) aaaaGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑A‑ttttttttgc < 1:358225/133‑4 (MQ=255) cGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTACAAGAAAAGAAATTAA‑TTTTTTA‑TCGTAGAt > 1:157500/1‑133 (MQ=255) gTTTGTCACGTCTTGGGCGCAAGACATATCGCCGCAATCAGGGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACa < 2:164803/133‑1 (MQ=255) gTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACa < 1:474284/133‑1 (MQ=255) tGTCACGGGTTGTGCTCATGACCTATCGCACCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCg < 1:210817/133‑1 (MQ=255) tGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGGGCTTTTTGATAGCGGAAAACGGAGAGTTAAAAGATAACCAAAT‑A‑TTTTTTTCGCGTAGATAACAGCg > 2:143138/1‑133 (MQ=255) cTTGTGCGTAAGACATATTGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACg < 1:163799/133‑1 (MQ=255) cTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACg < 1:57596/133‑1 (MQ=255) cTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACg < 2:329966/133‑1 (MQ=255) tgtgCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTg < 1:267504/133‑1 (MQ=255) tgtgCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTg < 1:281624/133‑1 (MQ=255) tgtgCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTg < 1:511985/133‑1 (MQ=255) gcAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGttt < 1:354130/133‑1 (MQ=255) aGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATCACAGCGTAGTTACATGGGTTTTaa > 2:145188/1‑133 (MQ=255) atatCGCAGCTATCACCGACGGGCCAGAGGAAAGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACt < 2:157500/133‑1 (MQ=255) agcaATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGGGTATTTACGTGGGTGTTAATACTTTGGTAt > 1:169107/1‑133 (MQ=255) agcaATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTAt > 2:66875/1‑133 (MQ=255) gcaATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACATCGTATATCCGTGGGTTTTAATCGTTTGGTacg > 2:238814/1‑131 (MQ=255) cGACGGGCAAGAAGAATGCCTGTCTGGTGCGTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTaaaa < 1:172435/133‑1 (MQ=255) cGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATATCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTaaaa < 2:33497/133‑1 (MQ=255) aCGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTaaaaaa < 1:483295/133‑1 (MQ=255) ggCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGaa > 2:390535/1‑133 (MQ=255) tGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGCTATAAACTAAAAGAGAACAATACAAcgc > 2:192840/1‑133 (MQ=255) gACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAAcgcg < 1:66875/133‑1 (MQ=255) gACTGTCTCGTGCTTTTTCATACCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAAcgcg < 1:143138/133‑1 (MQ=255) tCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGt < 2:431199/133‑1 (MQ=255) ggTGCTTTTTGATAGGGGTAAAAGGAGATTGAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAg < 1:236383/133‑1 (MQ=255) gTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGt < 1:264291/133‑1 (MQ=255) tttttGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGttatt > 1:426360/1‑133 (MQ=255) tttGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGCGAGTTATTaa > 2:373112/1‑133 (MQ=255) tAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAata > 1:64311/1‑133 (MQ=255) gCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATAtt > 1:63407/1‑133 (MQ=255) ggAAAACGCAGATTTCAACGAAAACAAAATCT‑TTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTc < 1:238814/133‑1 (MQ=255) gAAAACGGAGATTTAAAAAAAAGCAAAAT‑ACTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTcc < 1:192840/133‑1 (MQ=255) aaaaCGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCg < 2:426360/133‑1 (MQ=255) aCGGAGATTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAg < 2:64311/133‑1 (MQ=255) gaTTTAAAAGAAAACAAAAT‑ATTTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATa < 2:63407/133‑1 (MQ=255) | CAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATTTAAAAGAAAACAAAAT‑A‑TTTTTTT‑GCGTAGATAACAGCGTATTTACGTGGGTTTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATA > NC_000913/1103614‑1103849 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |