Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,162,570 | G→A | 18.5% | intergenic (+474/‑106) | PP_2776 → / → PP_2777 | homocysteine S‑methyltransferase family protein/acyl carrier protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,162,570 | 0 | G | A | 18.5% | 52.2 / 7.6 | 27 | intergenic (+474/‑106) | PP_2776/PP_2777 | homocysteine S‑methyltransferase family protein/acyl carrier protein |
Reads supporting (aligned to +/- strand): ref base G (11/11); new base A (3/2); total (14/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01 |
GCTTCTGGAACTATCGTGCTGTTTCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCCATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAGTA > NC_002947/3162421‑3162715 | gCTTCTGGAACTATCGTGCTGTTTCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCCATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAg < 1:405082/150‑1 (MQ=255) tGGAACTATCGTGCTGTTTCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCCATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGga < 1:273338/150‑1 (MQ=255) aTCGTGCTGTTTCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCCATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTa < 1:273014/150‑1 (MQ=255) aGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCCATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAg < 1:170727/124‑1 (MQ=255) aGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCCATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAg > 2:170727/1‑124 (MQ=255) gACCCCTCTCCAAGCCCCATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGc < 1:92129/150‑1 (MQ=255) aaGCCCCATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGcc > 1:754455/1‑150 (MQ=255) ccccATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGa > 2:662056/1‑150 (MQ=255) cTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTatt < 2:508943/150‑1 (MQ=255) tCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATg > 2:711956/1‑150 (MQ=255) ggTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCt < 1:711956/150‑1 (MQ=255) gggCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGTCTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAg < 2:754455/150‑1 (MQ=255) tttCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCttttttt < 1:295050/99‑1 (MQ=255) tttCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCttttttt > 2:295050/1‑99 (MQ=255) tcACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGaat > 2:213571/1‑150 (MQ=255) gCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTg < 1:213571/150‑1 (MQ=255) tGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAAt > 2:537066/1‑150 (MQ=255) tGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAAt > 2:680386/1‑150 (MQ=255) gATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAAtt < 1:537066/150‑1 (MQ=255) gATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTATGTAGGGAAACACTCAATGAATAATCCGTTGGAAtt < 1:662056/150‑1 (MQ=255) ggaaCGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGa > 1:798011/4‑150 (MQ=255) tctaaagtctcCCCCTAAGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGc > 2:605667/8‑150 (MQ=255) taaagtctcCCCCTAAGGGGAGAGTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCgg > 1:720207/6‑101 (MQ=255) taaagtctcCCCCTAAGGGGAGAGTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCgg < 2:720207/96‑1 (MQ=255) aagtctcCCCCTAAGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGtt < 1:748671/130‑1 (MQ=255) aagtctcCCCCTAAGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGtt > 2:748671/4‑133 (MQ=255) tctcCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGtcatca < 1:680386/150‑1 (MQ=255) ccTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAGt > 1:276984/1‑150 (MQ=255) cTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTAGCCCGATGACTTTATTCATTCAATCCGGCTTTTGTAGGAAAACACCCAATGAATAATCCGTTGGAATTGGACAAAGTCATGAGCAGTa > 2:237255/1‑150 (MQ=255) cTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGAGCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGTCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAGTa > 1:662245/1‑150 (MQ=255) | GCTTCTGGAACTATCGTGCTGTTTCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCCATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAGTA > NC_002947/3162421‑3162715 |
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: ‑ |