Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,123,943 | T→A | 16.2% | intergenic (‑138/+83) | PP_mr13 ← / ← gcvT‑I | cspA regulatory region/aminomethyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,123,943 | 0 | T | A | 16.2% | 92.1 / 11.8 | 37 | intergenic (‑138/+83) | PP_mr13/gcvT‑I | cspA regulatory region/aminomethyltransferase |
Reads supporting (aligned to +/- strand): ref base T (17/14); new base A (2/4); total (19/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.05e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01 |
GGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCTTTCCA > NC_002947/1123798‑1124092 | ggCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCgcctggcc > 1:289131/1‑150 (MQ=255) ccAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCAAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGcctggcctgg > 1:157621/1‑150 (MQ=255) cACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGTCGATCACCGACCGATTTGTCGGAAACGGCATGGCGCCAGGCCAGGCGCCAt < 1:185519/150‑1 (MQ=255) cACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCAGGCCAGGCGCCAt < 1:163359/150‑1 (MQ=255) gAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGtt > 2:11779/1‑150 (MQ=255) gACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTAt > 2:735693/1‑150 (MQ=255) aCCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTAtt < 2:252432/150‑1 (MQ=255) ccATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTgg < 1:748299/150‑1 (MQ=255) ccATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTgg > 1:600614/1‑150 (MQ=255) tGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCa > 1:644064/1‑150 (MQ=255) acacAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTc < 1:11779/150‑1 (MQ=255) tgcagtgACTCTACGCAAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAAc < 2:157621/150‑1 (MQ=255) tgACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGt < 2:313862/125‑1 (MQ=255) tgACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGt > 1:313862/1‑125 (MQ=255) tctACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATg < 2:123449/148‑1 (MQ=255) tctACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATg > 1:123449/1‑148 (MQ=255) tctACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTc > 1:59068/1‑150 (MQ=255) gCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTAtgc < 2:437911/150‑1 (MQ=255) gAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGc < 1:110986/114‑1 (MQ=255) gAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGc > 2:110986/1‑114 (MQ=255) gAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTc > 1:493407/1‑150 (MQ=255) gAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCAACTTGTCGAAAATGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTc > 2:108384/1‑150 (MQ=255) gTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATAGGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGc > 2:421776/1‑150 (MQ=255) ggCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACgtg < 2:289131/150‑1 (MQ=255) tGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCAc < 1:735693/150‑1 (MQ=255) cGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGaa > 1:121736/1‑150 (MQ=255) gATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGaaa < 2:59068/150‑1 (MQ=255) ccgaccgaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCAt < 1:428244/150‑1 (MQ=255) aTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACgtg < 1:437283/127‑1 (MQ=255) aTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACgtg > 2:437283/1‑127 (MQ=255) aaaCGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACgg > 1:648174/1‑108 (MQ=255) aaaCGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACgg < 2:648174/108‑1 (MQ=255) aacATGGCGCCAGGCCAGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGAc > 2:151936/3‑139 (MQ=255) aacATGGCGCCAGGCCAGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGAc < 1:151936/137‑1 (MQ=255) acATGGCGCCAGGCCAGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACgg < 1:674536/102‑1 (MQ=255) acATGGCGCCAGGCCAGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACgg > 2:674536/2‑103 (MQ=255) aTGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACct > 2:502618/1‑150 (MQ=255) aTGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACct > 2:697424/1‑150 (MQ=255) aTGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACct > 1:718962/1‑150 (MQ=255) ctggcctggcGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCTTTcc < 2:600614/150‑1 (MQ=255) tggcctggcGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCTTTCCa < 1:502618/150‑1 (MQ=255) tggcctggcGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCTTTCCa < 1:697424/150‑1 (MQ=255) | GGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCTTTCCA > NC_002947/1123798‑1124092 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |