Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,301,621 | G→C | 25.1% | intergenic (+67/‑1233) | PP_3774 → / → PP_5627 | hypothetical protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,301,621 | 0 | G | C | 25.1% | 22.0 / 6.3 | 16 | intergenic (+67/‑1233) | PP_3774/PP_5627 | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (6/6); new base C (2/2); total (8/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01 |
CCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTCTGC > NC_002947/4301473‑4301765 | ccAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCgg > 2:4125/1‑150 (MQ=255) cgacacCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCTTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGgggctgggct > 1:498291/3‑150 (MQ=255) tGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCtt > 2:233374/1‑150 (MQ=255) gAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTAc < 2:119883/150‑1 (MQ=255) cTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGgc < 2:277035/150‑1 (MQ=255) tGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACATTCAACGGGAAGTTACTCCAGTGGGCGCTTGTGATaa < 1:247436/150‑1 (MQ=255) cccAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACt > 1:227101/1‑150 (MQ=255) tCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGccc > 1:477322/1‑150 (MQ=255) tgatagaaGCCCAGCCCCGCGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAATTTACTACAGTGGGCGCTTGTGATAACTATAACGCTTGTTATTGAATGCGCGCGGTGGGTCCGGCggg > 2:390237/7‑150 (MQ=255) tgatagaaGCCCAGCCCCGCGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCggg > 2:136829/7‑150 (MQ=255) gatagaaGCCCAGCCCCGCGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGc < 1:136829/145‑1 (MQ=255) gatagaaGCCCAGCCCCGCGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTATACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGc < 1:390237/145‑1 (MQ=255) aaaaaGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCt > 2:152951/1‑150 (MQ=255) aGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCcggcg < 1:494124/136‑1 (MQ=255) aGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCcggcg > 2:494124/1‑136 (MQ=255) aGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTctgc < 1:352699/150‑1 (MQ=255) aGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTctgc < 1:233374/150‑1 (MQ=255) | CCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTCTGC > NC_002947/4301473‑4301765 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |