Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 4,589,054 | T→G | T69P (ACC→CCC) | opgB ← | phosphatidylglycerol‑‑membrane‑oligosaccharide glycerophosphotransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,589,054 | 0 | T | G | 95.5% | 71.2 / ‑5.6 | 22 | T69P (ACC→CCC) | opgB | phosphatidylglycerol‑‑membrane‑oligosaccharide glycerophosphotransferase |
Reads supporting (aligned to +/- strand): ref base T (0/0); major base G (7/14); minor base A (1/0); total (8/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.64e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTT > NZ_CP009273/4588921‑4589163 | agcagGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGACAGGc > 2:104407/1‑139 (MQ=255) tAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTa < 2:54835/139‑1 (MQ=255) tAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTa < 2:395048/139‑1 (MQ=255) tAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTa < 2:229878/139‑1 (MQ=255) tAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTa < 2:134722/139‑1 (MQ=255) aCCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAg < 2:391517/139‑1 (MQ=255) gatgatGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACcgc < 2:368642/139‑1 (MQ=255) ggCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTgg < 1:385829/114‑1 (MQ=255) ggCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTgg > 2:385829/1‑114 (MQ=255) cgcAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGc < 2:157344/139‑1 (MQ=255) agcGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGt > 2:145619/1‑136 (MQ=255) agcGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGt < 1:145619/136‑1 (MQ=255) cgcACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATc < 2:453036/120‑1 (MQ=255) cgcACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTg > 2:481671/1‑139 (MQ=255) cgcACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGCATAGAGAACCGCGTCGTTAATACCATc > 1:453036/1‑120 (MQ=255) tGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAg < 1:478916/59‑1 (MQ=255) tGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAg > 2:478916/1‑59 (MQ=255) tGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAg > 1:225680/1‑139 (MQ=255) cccGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCa < 1:481671/139‑1 (MQ=255) cAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGtt < 2:225680/139‑1 (MQ=255) aGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCa < 1:115849/113‑1 (MQ=255) aGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCa > 2:115849/1‑113 (MQ=255) | AGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTT > NZ_CP009273/4588921‑4589163 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |