Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,381,402 | T→G | 6.9% | intergenic (+31/+56) | oprF → / ← cobA | porin F/uroporphyrinogen‑III C‑methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,381,402 | 0 | T | G | 6.9% | 172.5 / 4.6 | 58 | intergenic (+31/+56) | oprF/cobA | porin F/uroporphyrinogen‑III C‑methyltransferase |
Reads supporting (aligned to +/- strand): ref base T (28/26); new base G (2/2); total (30/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.90e-01 |
AGCCGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCAAACGCC > NC_002947/2381255‑2381548 | aGCCGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTc < 8:119563/149‑1 (MQ=255) ttGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCgg < 7:86233/148‑1 (MQ=255) cTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGtttt > 1:17201/1‑149 (MQ=255) cTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGtttt < 6:283459/149‑1 (MQ=255) gCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGc < 7:37006/148‑1 (MQ=255) tACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGcc < 5:45076/132‑1 (MQ=255) tACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGcc > 6:45076/1‑132 (MQ=255) tACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCt > 4:170473/1‑148 (MQ=255) tACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTg > 3:320835/1‑149 (MQ=255) aCGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGACAAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGCCCTCGGCCGGGTTTTTCTTTGCCTgg > 6:181358/1‑149 (MQ=255) cGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCt < 3:170473/146‑1 (MQ=255) gAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGttt > 5:268607/1‑148 (MQ=255) ccGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCg < 1:104313/149‑1 (MQ=255) ccGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCg > 2:188640/1‑149 (MQ=255) ccGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCg > 4:209308/1‑149 (MQ=255) tCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGcc < 5:297578/149‑1 (MQ=255) ccGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGccc < 6:268607/149‑1 (MQ=255) cTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAgg < 2:17201/148‑1 (MQ=255) tGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCGGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCATGACGGcc > 8:295559/1‑149 (MQ=255) aCAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGcc < 4:269771/147‑1 (MQ=255) aaCGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTc > 4:302386/1‑149 (MQ=255) cACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTc > 2:140079/1‑109 (MQ=255) cACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCg > 6:99090/1‑149 (MQ=255) cACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCAAGAAAAACCCGGCCGAGGCCGGGtttt > 1:79792/1‑107 (MQ=255) cACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCAAGAAAAACCCGGCCGAGGCCGGGtttt < 2:79792/107‑1 (MQ=255) ccGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGt < 8:45185/149‑1 (MQ=255) cGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGc > 4:290561/1‑111 (MQ=255) aagaagGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCCTGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAg < 7:75963/149‑1 (MQ=255) aagaagGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAg > 6:151804/1‑149 (MQ=255) aagGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCTTCAGGGCCggt < 7:162193/135‑3 (MQ=255) aagGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCTTCAGGGCCggt > 8:162193/1‑133 (MQ=255) ccGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCtgcga < 1:188640/149‑1 (MQ=255) tGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTc < 4:205158/126‑1 (MQ=255) tGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTc > 3:205158/1‑126 (MQ=255) tGTCAACCGTAGCGTTGAAGCTCAGGTAGAATCCAAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACACGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGAACTGAAGGcagatc > 3:205159/1‑121 (MQ=255) aaCCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATcc < 2:178472/148‑1 (MQ=255) aCCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGc < 4:163162/149‑1 (MQ=255) ccGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGcc > 4:85065/1‑149 (MQ=255) cgcgTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGCTTGCGATCCGCCCCa > 8:73585/1‑148 (MQ=255) gcgTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCa > 4:125316/1‑147 (MQ=255) gcgTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATa < 3:85065/149‑1 (MQ=255) gTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGc > 1:314482/1‑149 (MQ=255) tGAAGCTCAGGTAGAAGCCCAGGCCAAGTAAATGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGACTGGGTTTTAACCGGACCTGAAGGCAGCTCCTCCGGTCAGCCTGCGAATGCGATCCGCCCCATAGCct > 5:212187/1‑149 (MQ=255) aaGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCctgc > 4:2246/1‑149 (MQ=255) tCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCAAAGCCTgctgct > 4:237166/1‑149 (MQ=255) ggCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCGGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTccc < 4:320835/149‑1 (MQ=255) aaTTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTccc < 4:43293/141‑1 (MQ=255) aaTTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTccc > 3:43293/1‑141 (MQ=255) gTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGa < 2:314482/148‑1 (MQ=255) gAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAACCACCAGGATTGc < 5:99090/148‑1 (MQ=255) gCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCggg < 5:181358/149‑1 (MQ=255) aaagAAAAACCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTAc > 4:288020/4‑149 (MQ=255) aTGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACg < 3:125316/149‑1 (MQ=255) tGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCCATTGCGATCCGCCCCAGAGCCTGCTGCACCCCACCCACCGCCCCAAACCCCAGGCATCCCCGGCTACGc > 6:273311/1‑149 (MQ=255) tGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCAATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACg > 6:134811/1‑148 (MQ=255) aaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCCCCTCCCCAATCACCAGGATTGCCGGGCTACGcaaac > 3:174408/1‑149 (MQ=255) ccGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaaccaaa < 3:209308/149‑1 (MQ=255) ccTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCAAACGcc < 3:302386/149‑1 (MQ=255) | AGCCGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCAAACGCC > NC_002947/2381255‑2381548 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |