Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 4,006,495 | C→G | A28G (GCG→GGG) | metE → | 5‑methyltetrahydropteroyltriglutamate‑‑ homocysteine S‑methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,006,495 | 0 | C | G | 98.2% | 192.6 / ‑6.0 | 55 | A28G (GCG→GGG) | metE | 5‑methyltetrahydropteroyltriglutamate‑‑ homocysteine S‑methyltransferase |
Reads supporting (aligned to +/- strand): ref base C (0/1); new base G (30/24); total (30/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.55e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGCGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCACTGGGATCAACAAAAGCAAGCGGGTATC > NZ_CP009273/4006411‑4006580 | aaaTGACAATATTGAATCACACCCTCGGTTTCCCTCTCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGa < 2:95779/90‑1 (MQ=255) aTGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGAGCAAGAAAGTTATTGGGGGGGGAAc < 1:661300/90‑1 (MQ=255) tGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACt > 2:871553/1‑90 (MQ=255) atatTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGAGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCAcg < 2:19173/90‑1 (MQ=255) tGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTg > 1:100179/1‑90 (MQ=255) gAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTga > 2:861448/1‑90 (MQ=255) gAATCACACCCTCGGTTTCCCTCGCGTTGGACTGCGTCGAGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTcc > 2:718040/1‑82 (MQ=255) gAATCACACCCTCGGTTTCCCTCGCGTTGGACTGCGTCGAGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTcc < 1:718040/82‑1 (MQ=255) acacCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTgaa > 2:867462/1‑86 (MQ=255) acacCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTgaa < 1:867462/86‑1 (MQ=255) acacCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAAc > 1:1020615/1‑90 (MQ=255) acacCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAAc > 2:506998/1‑90 (MQ=255) acacCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAAc > 2:468582/1‑90 (MQ=255) acacCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAAc > 2:389263/1‑90 (MQ=255) acacCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAAc > 1:618684/1‑90 (MQ=255) cTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAActgctg > 1:11259/1‑90 (MQ=255) cTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAActgctg > 1:322363/1‑90 (MQ=255) cTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAActgctg > 1:808194/1‑90 (MQ=255) cTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAActgctg > 2:855809/1‑90 (MQ=255) tttCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGt > 1:739323/1‑90 (MQ=255) ttCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTa > 2:72163/1‑90 (MQ=255) cccTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAgg > 1:343980/1‑90 (MQ=255) cTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGc < 1:506998/90‑1 (MQ=255) cgcgTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGa > 2:795903/1‑55 (MQ=255) cgcgTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGa < 1:795903/55‑1 (MQ=255) cgcgTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGGAGGGCgg > 1:114064/1‑89 (MQ=255) cgTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGc < 1:165083/78‑1 (MQ=255) cgTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGc > 2:165083/1‑78 (MQ=255) tGGCCTGCGTCGCGAGCTGAAAAAAGCGAAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGGGAAtt > 2:719937/1‑90 (MQ=255) ggCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGACCTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAAATg > 2:529585/1‑90 (MQ=255) ggCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTa < 1:645921/78‑1 (MQ=255) ggCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTa > 2:645921/1‑78 (MQ=255) ggCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTg > 2:1049053/1‑90 (MQ=255) ggCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTg > 1:262678/1‑90 (MQ=255) ccTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCg < 2:978567/90‑1 (MQ=255) gCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTa > 1:555086/1‑73 (MQ=255) gCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGCAACTCCCCGCGTGAAGACCTGCTGGCGGTa < 2:555086/73‑1 (MQ=255) gTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTc < 2:1032259/90‑1 (MQ=255) gTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTc < 1:535453/90‑1 (MQ=255) gTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTc < 2:863586/90‑1 (MQ=255) cgcgAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGt < 2:10626/90‑1 (MQ=255) gAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCAc < 1:467510/90‑1 (MQ=255) gAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCAc < 1:1072208/90‑1 (MQ=255) gAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCAc < 1:259963/90‑1 (MQ=255) aGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCACt > 1:1063006/1‑90 (MQ=255) cTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGTCGTGAATTGCGTGCTCGTCACTgg > 1:39188/1‑90 (MQ=255) cTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCACTgg < 2:343980/90‑1 (MQ=255) gAAAAAAGCGCAAGAAAGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGATCGTCACTGGGa < 2:827519/90‑1 (MQ=255) gcAAGAAAGTTATTGGGGGGGGACCTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCACTGGGATCAACaaaa < 1:529585/90‑1 (MQ=255) aaGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCACTGGGATCAACAAaagcaagc < 1:72163/90‑1 (MQ=255) aaGTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCACTGGGATCAACAAaagcaagc < 2:1020615/90‑1 (MQ=255) gTTATTGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCACTGGGATCAACAAAAGCAAGCgg > 1:469013/1‑90 (MQ=255) tGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTa > 1:657499/1‑39 (MQ=38) tGGGGGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTa < 2:657499/39‑1 (MQ=38) tGGGCGGGGAACTCCCCGCGTGAAGAACTGCTGGGGGTAGGGCGTGAATTGCGTGCTCTTCACTGGGCTCTTCAAAAGCAAGCGGGTATc < 2:579613/90‑1 (MQ=255) | AAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGCGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCACTGGGATCAACAAAAGCAAGCGGGTATC > NZ_CP009273/4006411‑4006580 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |