Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,395,630 | G→C | 21.1% | intergenic (‑195/+48) | dacB ← / ← PP_2099 | D‑alanyl‑D‑alanine carboxypeptidase/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,395,630 | 0 | G | C | 21.1% | 74.7 / 17.8 | 38 | intergenic (‑195/+48) | dacB/PP_2099 | D‑alanyl‑D‑alanine carboxypeptidase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (18/12); new base C (3/5); total (21/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.26e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.03e-01 |
GTGAATACCCTGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGACTTCGATCTTGGTCAGTTCGCTGCGGCCTTTGAGCCAGG > NC_002947/2395483‑2395772 | gTGAATACCCTGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATggg < 2:1136370/150‑1 (MQ=255) gAATACCCTGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGt < 2:452947/150‑1 (MQ=255) tACCCTGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATGAAAAGGCCACCCCATTCGGGTGGc < 2:95645/150‑1 (MQ=255) tACCCTGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCCATTCGGGTGGc > 2:224121/1‑150 (MQ=255) cTGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCCATTCGGGTGGCCttt > 1:1201251/1‑150 (MQ=255) cTGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCCATTCGGGTGGCCttt < 2:1335171/150‑1 (MQ=255) cggcgAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCCATTCGGGTGGCCttttta < 2:1201251/150‑3 (MQ=255) cgAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCCATTCGGGTGGCCtttttattg < 1:618021/150‑6 (MQ=255) aGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCCATTCGGGTGGCCtttttattggc < 1:224121/150‑8 (MQ=39) ggggCTGTAAGAGATGATGGGCTCACGACGTAAACAGTACGCAGAATGCCCCAAGCGCAGCCGCGATGCGTAAGGTTCGACGGAAAAGCGTCGCGAATAAAAAGGCCACCCGCATGGGGTGGCCTTTTCCGCAATCAACCTTTCAGCCCg > 2:633256/1‑150 (MQ=255) ggCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGaa > 2:709384/1‑150 (MQ=255) gCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAg > 1:478223/1‑150 (MQ=255) gTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGCCGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCACCCCGAAGCTg > 1:197434/1‑150 (MQ=255) gagaTGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACaa > 2:1131817/1‑150 (MQ=255) agaTGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAg > 2:961207/1‑150 (MQ=255) gatgatGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGt > 2:1003113/1‑150 (MQ=255) acgTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCATTCAGCCttgg < 2:1236053/150‑2 (MQ=255) tACGCATTATGCCCCAAGCGCAGCAGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCtt > 2:1004670/1‑150 (MQ=255) ttATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCgg < 1:709384/150‑1 (MQ=255) ttATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCgg < 1:961207/150‑1 (MQ=255) ttATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCgg < 2:478223/150‑1 (MQ=255) ttATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCgg < 1:1131817/150‑1 (MQ=255) ccAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAg > 1:665590/1‑150 (MQ=255) aaGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGc > 1:950959/1‑150 (MQ=255) tCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGAc > 2:1108328/1‑150 (MQ=255) ttCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGACTTCGa > 1:709459/1‑150 (MQ=255) tCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGACTTCGAt < 1:1126567/150‑1 (MQ=255) cGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGACTTCGATCTTg < 2:665590/150‑1 (MQ=255) tCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGACTTCGATCTTGGTCAGTTCGc > 1:204455/1‑150 (MQ=255) ccAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGACTTCGATCTTGGTCAGTTCGCTGc > 1:39208/1‑150 (MQ=255) aaTAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGACTTCGATCTTGGTCAGTTCGCTGCgg > 1:1303463/1‑150 (MQ=255) aTNAAAAGGCCACCCGAATGGGGTGGCNTTTTCAGCAATCNACCTTTNANCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGACTTCGATCTTGGTCAGTTCGCTGCGGc > 1:940816/1‑150 (MQ=255) tgaaaaGGCCACCCCATTCGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGACTTCGATCTTGGTCAGTTCGCTGCGGcc > 1:964460/3‑150 (MQ=255) aGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGACTTCGATCTTGGTCAGTTCGCTGCGGCCTTTGa < 2:39208/150‑1 (MQ=255) gCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGACTTCGATCTTGGTCAGTTc < 1:535704/134‑1 (MQ=255) gCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGACTTCGATCTTGGTCAGTTc > 2:535704/1‑134 (MQ=255) aCCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGACTTCGATCTTGGTCAGTTCGCTGCGGCCTTTGAGCCAg > 1:863674/1‑150 (MQ=255) aCCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGACTTCGATCTTGGTCAGTTCGCTGCGGCCTTTGAGCCAg < 2:270721/150‑1 (MQ=255) cccGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGACTTCGATCTTGGTCAGTTCGCTGCGGCCTTTGAGCCAgg < 1:1003113/150‑1 (MQ=255) ggggTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCgg > 1:658387/1‑78 (MQ=255) ggggTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCgg < 2:658387/78‑1 (MQ=255) | GTGAATACCCTGCGGCGAGGGAAAAGGGCTTTGGGGCTGTAAGAGATGATGGCCCCACGACGTAAACAGTACGCATTATGCCCCAAGCGCAGCCGCGATGCGTCACGTTCGACGGAAAAGCATCGCCAATAAAAAGGCCACCCGAATGGGGTGGCCTTTTCAGCAATCAACCTTTCAGCCCGAAGCTGACAAGTATCACTCAGCCTTGTTGATAGGCTTTTCCGGGTACCAGGCATCCAGCAGCGGGCTGACTTCGATCTTGGTCAGTTCGCTGCGGCCTTTGAGCCAGG > NC_002947/2395483‑2395772 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |