Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,381,315 | T→C | 18.8% | intergenic (+32/+33) | bamE → / ← PP_4732 | outer membrane protein assembly factor BamE/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,381,315 | 0 | T | C | 18.8% | 61.1 / 10.0 | 32 | intergenic (+32/+33) | bamE/PP_4732 | outer membrane protein assembly factor BamE/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (11/15); new base C (3/3); total (14/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.97e-01 |
CTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGA > NC_002947/5381174‑5381461 | cTGAGGAACAGCCGGCCAAGCCGTGTTCGACCGCAGAGTCCATCCAGCTCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGAGGGAACATCAGATACAAAAAAGCCCGGACTTGGTCCgg < 2:182193/150‑1 (MQ=255) tGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCggg < 1:143407/150‑1 (MQ=255) aaaaaCCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCtttt > 1:283733/1‑150 (MQ=255) aaaCCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGAATACCCCGGTCCCGACCCCGGCACCGCGGGAAACCCGCTAAGCGGGAACATCAGATGCAAAATAGCCCGGACTTGGTCCGGGCtttttt > 1:454480/1‑150 (MQ=255) aCCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTtttgt < 1:121691/150‑1 (MQ=255) cAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTgcgc < 2:174784/150‑1 (MQ=255) cTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACCAAGTCCGGGCtt > 2:444799/1‑133 (MQ=255) cTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACCAAGTCCGGGCtt < 1:444799/133‑1 (MQ=255) gTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGACACCGCTGGAAACCCGCTAGGCGGGAACATCAGATGCAAAAAAGCACGGACTTGGTCCGGTGTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAAt > 2:36235/1‑150 (MQ=255) ccatccaGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTg < 2:309400/150‑1 (MQ=255) atccaGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCt > 2:316009/1‑150 (MQ=255) agaTCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTg < 2:283733/150‑1 (MQ=255) gaTCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTgg < 1:248846/150‑1 (MQ=255) acacCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGaa < 1:316009/150‑1 (MQ=255) gagaCTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGa > 1:336854/1‑150 (MQ=255) tACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGc < 2:443746/150‑1 (MQ=255) ccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACCAAGTCCGGGCTTTtttgcatc > 1:41537/1‑85 (MQ=25) ccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACCAAGTCCGGGCTTTtttgcatc < 2:41537/89‑5 (MQ=25) ggTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTc < 2:222014/150‑1 (MQ=255) cTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCAtt > 2:19540/1‑88 (MQ=255) cTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCAtt < 1:19540/88‑1 (MQ=255) cTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTg < 1:225182/94‑1 (MQ=255) cTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTg > 2:225182/1‑94 (MQ=255) aCCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCa > 2:377199/1‑150 (MQ=255) aCCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCa > 2:445854/1‑150 (MQ=255) tCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCg > 2:296275/1‑150 (MQ=255) aaaaaaGCCCGGACCAAGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCgg < 2:282180/119‑1 (MQ=39) aaaaaaGCCCGGACCAAGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCgg > 1:282180/1‑119 (MQ=39) cccGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGa < 2:358745/105‑1 (MQ=255) cccGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGa > 1:358745/1‑105 (MQ=255) gACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGa < 1:445854/150‑1 (MQ=255) gACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGa < 1:377199/150‑1 (MQ=255) | CTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGA > NC_002947/5381174‑5381461 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |