Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,861,127 | A→G | 12.9% | intergenic (‑78/+159) | ligA ← / ← PP_4275 | DNA ligase/cell division protein ZipA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,861,127 | 0 | A | G | 12.9% | 72.5 / 3.1 | 31 | intergenic (‑78/+159) | ligA/PP_4275 | DNA ligase/cell division protein ZipA |
Reads supporting (aligned to +/- strand): ref base A (10/17); new base G (2/2); total (12/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.30e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01 |
CACGTAGTAGCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCG > NC_002947/4860981‑4861247 | cACGTAGTAGCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGc > 1:269725/1‑149 (MQ=255) cACGTAGTAGCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGc > 7:161327/1‑149 (MQ=255) cACGTAGTAGCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACCGTGCCGGc > 1:25584/1‑149 (MQ=255) cGTAGTAGCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCct > 5:228291/1‑149 (MQ=255) gCGATAATTGTGCTGATCGGGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAATAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTcgc < 8:131247/148‑1 (MQ=255) gCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTcgcg < 2:269725/149‑1 (MQ=255) gATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCggg > 4:259141/1‑149 (MQ=255) aTAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGGCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGt < 5:194809/149‑1 (MQ=255) aTTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTaaa < 5:149727/149‑1 (MQ=255) ttGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAc < 6:228291/149‑1 (MQ=255) tgCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACCGTGCCGGCCTCTTCGCGGGTAAAccc < 2:25584/148‑1 (MQ=255) tGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTc > 8:250508/1‑149 (MQ=255) aTCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTcc > 3:131377/1‑148 (MQ=255) tCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑ca < 2:74269/149‑1 (MQ=255) gCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAg < 2:114126/149‑1 (MQ=255) tCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAggg > 2:267750/1‑149 (MQ=255) cACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAGGGCGTc < 4:131377/149‑1 (MQ=255) gTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAGGGCGTCACATGc < 4:233249/149‑1 (MQ=255) ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCCGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAGGGCGTCACATGCCCTAGGGCGGTGa < 5:279890/149‑1 (MQ=255) cggcggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAGGGCGTCACATGCCCTAGGGCGGTGATGTa < 3:259141/148‑1 (MQ=255) tCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAg < 7:289752/149‑1 (MQ=255) ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAGGGCGTCACATGCCCTAGGGCGGTg > 1:71087/1‑119 (MQ=255) ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAGGGCGTCACATGCCCTAGGGCGGTg < 2:71087/118‑1 (MQ=255) gcaaagcaaaAGAGCAGCTTGGGGTGGTGCGTGTTGTTCCA‑GGTGTCTGTTTACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAGGGCGTCACATGCCCTAGGTCGTTGATGTACCTGTAGTAGCGGGTTTACCCGCGa < 6:312223/146‑1 (MQ=39) caaagcaaaAGAGCAGCTTGGGGTGGTGCGGGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCTCACAAGGGCGTCACATGCCCTAGGGCGGTTATGTACCTGTAGGAGCGGGTTTACCCGCGaa < 1:267750/148‑1 (MQ=39) aagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAGGGCGACACATGCCCTAGGGCGGGGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGc < 5:242382/149‑1 (MQ=255) cagcTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCt < 2:39081/110‑1 (MQ=255) cagcTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCt > 1:39081/1‑110 (MQ=255) ggggTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCg > 2:255178/1‑149 (MQ=255) gtcattGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa > 7:229208/5‑101 (MQ=255) gtcattGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa < 8:229208/97‑1 (MQ=255) | CACGTAGTAGCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCA‑GGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCC‑‑CACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCG > NC_002947/4860981‑4861247 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |