Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 92,505 | C→G | A6G (GCC→GGC) | mraY → | phospho‑N‑acetylmuramoyl‑pentapeptide‑ transferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 92,505 | 0 | C | G | 85.0% | 117.0 / 14.2 | 40 | A6G (GCC→GGC) | mraY | phospho‑N‑acetylmuramoyl‑pentapeptide‑ transferase |
Reads supporting (aligned to +/- strand): ref base C (4/2); new base G (20/14); total (24/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-01 |
CGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCT > NZ_CP009273/92419‑92589 | cGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGa < 2:569770/90‑1 (MQ=255) cGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGa < 1:196398/90‑1 (MQ=255) ttAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACAtt < 1:901344/90‑1 (MQ=255) tAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACAttt < 1:819652/90‑1 (MQ=255) tAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACAttt < 1:302702/90‑1 (MQ=255) tAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACAttt < 2:390733/90‑1 (MQ=255) aaGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCa < 2:117380/90‑1 (MQ=255) tAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTcc < 2:705500/90‑1 (MQ=255) tAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTcc < 1:948899/90‑1 (MQ=255) aGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCg > 1:264720/1‑90 (MQ=255) aGTGCCGCCATGGAAGAGGCAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCg > 2:453718/1‑90 (MQ=255) gTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTCAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCgg > 2:811342/1‑90 (MQ=255) tGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGc > 1:676404/1‑90 (MQ=255) ccgccATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCtt < 2:664132/90‑1 (MQ=255) ccATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTaa < 2:382628/90‑1 (MQ=255) ccATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTaa < 1:707442/90‑1 (MQ=255) tagtaCGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTa < 1:941512/90‑1 (MQ=255) aCGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTg > 2:987590/1‑90 (MQ=255) gcgcTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGAc > 2:574695/1‑90 (MQ=255) ttACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTc > 2:617745/1‑90 (MQ=255) gaATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATc > 2:48231/1‑90 (MQ=255) gaATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATc > 2:452345/1‑90 (MQ=255) aTGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGt > 2:631976/1‑90 (MQ=255) aTGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGt > 1:437551/1‑90 (MQ=255) aTGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGt > 1:369574/1‑90 (MQ=255) aTGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGt > 2:91883/1‑90 (MQ=255) aTGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCCGACGTTTCGCGCCATCGt > 1:670094/1‑90 (MQ=255) aTGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGt > 2:760404/1‑90 (MQ=255) aTGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGt > 1:95898/1‑90 (MQ=255) aTGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGt > 1:665274/1‑90 (MQ=255) aTGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGt > 1:646298/1‑90 (MQ=255) gggACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCa < 2:439090/90‑1 (MQ=255) ggACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAg > 1:556662/1‑90 (MQ=255) cATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCct > 1:225593/1‑90 (MQ=255) tGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCctgc > 1:991218/1‑90 (MQ=255) tGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCctgc > 2:704137/1‑90 (MQ=255) gTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGAcc > 2:450639/1‑90 (MQ=255) ttGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACcgc < 2:556662/90‑1 (MQ=255) ttGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACcgc < 1:574695/90‑1 (MQ=255) ttGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACcgc > 1:424158/1‑90 (MQ=255) tggctgggCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACcgcg > 1:343223/1‑90 (MQ=255) ggctgggCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCgcgc < 1:720373/90‑1 (MQ=255) ggctgggCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCgcgc < 1:631976/90‑1 (MQ=255) gctgggCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCt > 1:760222/1‑90 (MQ=255) | CGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCT > NZ_CP009273/92419‑92589 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |