Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,530,025 | T→A | 29.4% | intergenic (+31/‑137) | PP_4018 → / → topB | acetyltransferase/DNA topoisomerase III |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,530,025 | 0 | T | A | 29.4% | 58.2 / 28.1 | 44 | intergenic (+31/‑137) | PP_4018/topB | acetyltransferase/DNA topoisomerase III |
Reads supporting (aligned to +/- strand): ref base T (19/12); new base A (5/8); total (24/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.99e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01 |
GGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGT > NC_002947/4529886‑4530171 | ggTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCcgcg > 8:162306/1‑148 (MQ=255) gggCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa < 6:213470/147‑1 (MQ=255) gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCGCCTGAGCCGGCCTCTTCGCGGGTAAACCCGCGaa < 5:1147/149‑1 (MQ=39) gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa < 5:158298/149‑1 (MQ=255) gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa < 6:201168/149‑1 (MQ=255) gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGaa < 7:162306/149‑1 (MQ=255) cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGGGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga < 6:245299/148‑1 (MQ=255) cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAagag < 2:186821/149‑1 (MQ=255) cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAagag < 1:158927/149‑1 (MQ=255) cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga < 1:311347/148‑1 (MQ=255) cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga < 1:180515/148‑1 (MQ=255) cgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGGGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGc < 5:194584/149‑1 (MQ=255) gccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga < 6:185702/145‑1 (MQ=255) gccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga > 5:185702/1‑145 (MQ=255) gccacgcaGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCACGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGcc < 6:413585/149‑1 (MQ=255) gccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga < 3:393293/140‑1 (MQ=255) gccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga > 4:393293/1‑140 (MQ=255) gcttcGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCt > 7:273853/5‑149 (MQ=255) gCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGTGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTa > 1:106328/1‑149 (MQ=255) gCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCTAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTa < 2:289548/149‑1 (MQ=255) gCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTa > 3:222618/1‑149 (MQ=255) aGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTAATCGCGGGGTTAACCCCGTATAGGCCGG‑CTCAGGCGATATACAGCCCCGCGGACAAAAAACCTGTAca > 8:152587/1‑148 (MQ=255) cTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc > 2:100480/1‑146 (MQ=255) cTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGGTTACCCGCGCAGAGGCCGG‑CTCAGGCGATATAAATCCCCGCTGACAAAAAACCTGGAcatccctc > 4:355149/1‑146 (MQ=255) tGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACt > 1:50946/1‑149 (MQ=255) cATCCGGCATGTTTAACTGGTGCTCGACTCGCGCGCCTGTTCCATAACCTGACGTGGACTCTTCGCGGTTTTACCGGCGATGAGGCCTG‑CTCAGGCGATATACATCCCAGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGc < 8:273853/149‑1 (MQ=255) cgcgACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCAt > 4:169849/1‑149 (MQ=255) cgACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATcc > 7:64089/1‑149 (MQ=255) cTGATCCATCACCTGCACTGGCCTCTTCGCCGGTTTACCCGCGAAGACGCCGG‑CTCAGGCGCTATACATACCCGCAGACAAAAAACGTGGACAACCCACCACAAAA‑AAAATGCCTCTGCGCCCCCCCATTCCCATCCAATGGCccccca > 8:170893/1‑149 (MQ=255) tGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCccacccc > 6:295806/1‑147 (MQ=255) gATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCaccacc > 6:369505/1‑149 (MQ=255) aTCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCAccc > 7:300575/1‑149 (MQ=255) tCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGAt > 1:76607/1‑148 (MQ=255) aCCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAACAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATg < 1:100480/148‑1 (MQ=255) aCCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGa < 2:76607/148‑1 (MQ=255) gCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCgg > 6:101981/1‑149 (MQ=255) gCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCgg > 3:340597/1‑149 (MQ=255) gCCTCTTCGCGGGTAAACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCCTCCTATGGCCCACCACCCCGATGCAGCGTCCCCgg > 6:411457/1‑149 (MQ=255) ctctTCGCGGGTAAACCCGCGAAGAGGCCGGTCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGc > 1:247862/1‑149 (MQ=255) ttCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATg > 2:373554/1‑149 (MQ=255) ttCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATg > 3:167100/1‑149 (MQ=255) cgGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGcc < 3:169849/149‑1 (MQ=255) ggttACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑ACCTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCCCCACCCCGAGGAAGCGTCCCCGGCCAACCGccagc > 5:314375/3‑146 (MQ=255) gTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAg < 7:290498/128‑1 (MQ=255) gTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAg > 8:290498/1‑128 (MQ=255) | GGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGG‑CTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAA‑AACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGT > NC_002947/4529886‑4530171 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |