Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,405,364 | T→C | 14.3% | M81V (ATG→GTG) | PP_3886 ← | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,405,364 | 0 | T | C | 14.3% | 82.0 / 6.6 | 35 | M81V (ATG→GTG) | PP_3886 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (17/13); new base C (2/3); total (19/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.42e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01 |
CCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGG > NC_002947/4405223‑4405506 | ccAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAg < 2:581378/150‑1 (MQ=255) cGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCGTTTAACGAAAAAGCCCCGGCATGTGCCGAGGCtt < 1:959618/150‑1 (MQ=255) cTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCGGCACATGCCGGGGCtttttc < 2:884396/150‑5 (MQ=255) tGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAGAGAGCCCGGCATGTGCCGAGGCTTGGAAt > 1:732/1‑150 (MQ=255) tGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGt < 2:10557/150‑1 (MQ=255) tGGGAACAGGGCGACCCGAGGTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCACAAATAAAAAAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGt < 2:83877/150‑1 (MQ=255) aGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCGGCACATGCCGGGGCtttttc > 1:1264163/1‑134 (MQ=255) aGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCGGCACATGCCGGGGCtttttc < 2:1264163/138‑5 (MQ=255) aGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATg < 1:586730/150‑1 (MQ=255) cccGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATggg < 2:1111685/135‑1 (MQ=255) cccGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATggg > 1:1111685/1‑135 (MQ=255) cccGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAAccc > 1:976087/1‑150 (MQ=255) ttCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCGGCACATGCCGGGGCtttttc > 1:700999/1‑121 (MQ=255) ttCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCGGCACATGCCGGGGCtttttc < 2:700999/125‑5 (MQ=255) ccgccgATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCCGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGt < 1:316942/150‑1 (MQ=255) gTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGgcg > 2:840662/1‑150 (MQ=255) gggCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCt < 2:976087/150‑1 (MQ=255) ggCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGt < 1:934569/124‑1 (MQ=255) ggCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGt > 2:934569/1‑124 (MQ=255) gCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTCTCGTGGCGCTTCCCTCCAGTCCCCACGCTgc > 2:300387/1‑149 (MQ=255) tCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGTAATGGGTGGAGATGGCAACCCTGTCAACTAATCACCGTGGCGGTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGa > 1:931340/1‑150 (MQ=255) tCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGa < 1:840662/150‑1 (MQ=255) cATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGatc > 1:431587/1‑150 (MQ=255) tGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGatcat > 1:457960/1‑150 (MQ=255) tGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCa < 2:471326/96‑1 (MQ=255) tGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCa > 1:471326/1‑96 (MQ=255) aTTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGc > 2:839831/1‑115 (MQ=255) aTTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGc < 1:839831/115‑1 (MQ=255) aTTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGccgccg > 2:697854/1‑127 (MQ=255) aTTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGccgccg < 1:697854/127‑1 (MQ=255) aTTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGAACATCTCCAAAAACg > 2:32313/1‑150 (MQ=255) gAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGc > 2:1192030/1‑150 (MQ=255) ccccGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACTCCTGCACGTTGCCGCCGGGCTTTGTTCATCTCCAAAAACGACAAAGCCCAGCTg > 1:211931/1‑150 (MQ=255) cccGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTgg > 2:30612/1‑150 (MQ=255) cccGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTgg > 1:301375/1‑150 (MQ=255) | CCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGG > NC_002947/4405223‑4405506 |
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: ‑ |