Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,385,149 | G→T | 13.7% | intergenic (‑119/+42) | proS ← / ← oprD | proline‑‑tRNA ligase/basic amino acid specific porin OprD |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,385,149 | 0 | G | T | 13.7% | 137.1 / 13.2 | 51 | intergenic (‑119/+42) | proS/oprD | proline‑‑tRNA ligase/basic amino acid specific porin OprD |
Reads supporting (aligned to +/- strand): ref base G (21/23); new base T (4/3); total (25/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.03e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01 |
AAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGA > NC_002947/1385011‑1385287 | aaaTATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCggg < 7:81232/148‑1 (MQ=255) ttGACTGGTGCGCAGGGTTAGCCGTGTCGATTGCCTCAGACGTAGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACGCCGGGCttttt < 4:43662/149‑1 (MQ=255) cgcATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATCGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCgg < 2:41300/148‑1 (MQ=255) aTGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGc < 5:215072/131‑1 (MQ=255) aTGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGc > 6:215072/1‑131 (MQ=255) gTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGTGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACGCCGGGCTTTTTc < 5:106642/134‑1 (MQ=38) gTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACGCCGGGCTTTTTc > 6:106642/1‑134 (MQ=39) aGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGAAGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGc > 1:88991/1‑148 (MQ=255) ctagtagATTTCTTAAGACGTTTAAATATCCCTGAATTTTACGGGGT‑GTAAACGAAGGCGTACAGGATGCCCCTTGGCGTGATACTTGTCAGCCAAGAAAAAGCTCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGcc < 4:43055/143‑1 (MQ=255) gtgtCGATTGCCTAAGACGTTGAATTAGGCGCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGGCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCa < 2:285317/149‑1 (MQ=255) aTTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTa > 3:32694/1‑149 (MQ=255) gCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACa < 4:222778/148‑1 (MQ=255) cTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGa > 6:187633/1‑149 (MQ=255) gTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAg > 5:236500/1‑149 (MQ=255) ttGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGc < 6:236500/149‑1 (MQ=255) aaaTAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCgg < 5:199762/148‑1 (MQ=255) ccGCATTGTACGGGGT‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCtttt < 6:333440/95‑1 (MQ=255) ccGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTTACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCAc < 6:50976/149‑1 (MQ=255) ccGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCtttt > 5:333440/1‑95 (MQ=255) aTTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATc < 6:111946/149‑1 (MQ=255) ggAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGc > 2:354454/1‑122 (MQ=255) ggAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGc < 1:354454/122‑1 (MQ=255) ggAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTtcg > 3:350306/1‑149 (MQ=255) cGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGtcgat < 4:350306/149‑1 (MQ=255) aaGGCGTACAGGATGCACCAGCACGTGACACTTGTCAGCCAAGAAAAAGCCCAGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCAATGCAAGCCAGGCTTACAGGAAGCTCAGCGGGAACTCCACCATCAGGCGCAATTCGtagatc > 2:82609/1‑148 (MQ=255) ggCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCGGGCTTACAGGATGCTCAGCGGGTACTCCAc < 6:48713/124‑1 (MQ=255) ggCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCGGGCTTACAGGATGCTCAGCGGGTACTCCAc > 5:48713/1‑124 (MQ=255) tACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATggg > 4:302010/1‑148 (MQ=255) gATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCg > 6:353806/1‑149 (MQ=255) ccccAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCgtagt < 1:82609/149‑1 (MQ=255) cccAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGCCGATCGATGGGGAATCgtagta < 5:187633/149‑1 (MQ=255) ccAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCgtag > 5:324288/1‑146 (MQ=255) cAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGa < 5:353806/149‑1 (MQ=255) gTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATc > 7:194481/1‑106 (MQ=255) gTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATc < 8:194481/106‑1 (MQ=255) acTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCTGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCgg > 1:253051/1‑147 (MQ=255) acTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCTGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCgg > 1:118786/1‑147 (MQ=255) acTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCgg > 5:33770/1‑147 (MQ=255) cTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCgg > 5:276163/1‑146 (MQ=255) cTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCgg > 7:245507/1‑146 (MQ=255) ttGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGgaag > 5:157748/1‑148 (MQ=255) aaGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCgg > 2:189225/1‑77 (MQ=255) aaGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCgg < 1:189225/77‑1 (MQ=255) gAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTgg < 6:276163/149‑1 (MQ=255) gAAAAAGCCCGGCAATCACGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAgg < 2:299453/91‑1 (MQ=37) gAAAAAGCCCGGCAATCACGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAgg > 1:299453/1‑91 (MQ=37) gAAAAAGCCCGGCAATCACGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTg > 7:136036/1‑148 (MQ=255) gAAAAAGCCCGGCAATCACGCCGGGCGTTTTCAGGTGCGGGGCCCTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGGTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTgg > 7:140750/1‑149 (MQ=255) taacaGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGgcg < 6:338109/88‑1 (MQ=255) aaaGCCCGGCGTGATTGCCGGGCTTTTTCTGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACACCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGCGGCCt < 2:118786/149‑1 (MQ=255) aGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGa < 2:88991/149‑1 (MQ=255) | AAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGA > NC_002947/1385011‑1385287 |
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: ‑ |