| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 1,810,378 | C→A | 31.2% | intergenic (+97/‑72) | eno → / → ftsB | enolase/cell division protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 1,810,378 | 0 | C | A | 31.2% | 43.5 / 23.6 | 32 | intergenic (+97/‑72) | eno/ftsB | enolase/cell division protein |
| Reads supporting (aligned to +/- strand): ref base C (13/9); new base A (5/5); total (18/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.12e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01 | |||||||||||
TGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGTGG > NC_002947/1810235‑1810518 | tGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAg > 1:572448/1‑150 (MQ=255) cgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCtt < 2:36742/150‑1 (MQ=255) aaaGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCTTc < 1:467910/147‑1 (MQ=255) aaaGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCTTc > 2:467891/1‑147 (MQ=255) aaaGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCTTc < 1:467891/147‑1 (MQ=255) aaaGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCTTc > 2:467910/1‑147 (MQ=255) cgGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGcccc < 2:241285/150‑1 (MQ=255) gCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTaa > 1:433282/1‑150 (MQ=255) gCTAATGAACAGATGGTAAAAATGCAGCAGCCGGAGTTGTAGGAACGTTAGTTCTGAATTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCa < 2:156851/150‑1 (MQ=255) aaGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGcc < 2:433282/150‑1 (MQ=255) cAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCttt > 2:155879/1‑150 (MQ=255) gaTGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAAc > 1:106516/1‑150 (MQ=255) gTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTgg > 2:452029/1‑150 (MQ=255) ccGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGACAGCCTTTTTAACCTGGATAACGTGATGCGCa > 2:400279/1‑150 (MQ=255) ggAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGt > 2:320288/1‑150 (MQ=255) ttGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCt > 1:282651/1‑150 (MQ=255) tAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCa > 1:470708/1‑141 (MQ=255) tAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCa < 2:470708/141‑1 (MQ=255) tGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATgc > 2:589998/1‑123 (MQ=255) tGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATgc < 1:589998/123‑1 (MQ=255) aaCTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTg > 2:355159/1‑108 (MQ=255) aaCTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTg < 1:355159/108‑1 (MQ=255) tACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGt > 2:326832/1‑51 (MQ=255) tACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGt < 1:326832/51‑1 (MQ=255) gTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGtt < 2:106516/150‑1 (MQ=255) tAACGGGTGGCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGt > 2:277163/1‑150 (MQ=255) aaGCCTGGCCGAGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTAtt < 2:1658/99‑1 (MQ=255) aaGCCTGGCCGAGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTAtt > 1:1658/1‑99 (MQ=255) aGCCTGGCCGAGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCtgtg > 2:540657/1‑150 (MQ=255) gCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTgtgg < 1:452029/150‑1 (MQ=255) gCCTGGCCGAGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTaa < 2:257296/45‑1 (MQ=21) gCCTGGCCGAGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTaa > 1:257296/1‑45 (MQ=21) | TGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGTGG > NC_002947/1810235‑1810518 |
| 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: ‑ |