Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,798,235 | A→T | 23.3% | intergenic (+28/+4) | waaL → / ← waaU | O‑antigen ligase/lipopolysaccharide core biosynthesis |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,798,235 | 0 | A | T | 23.3% | 127.7 / 41.3 | 61 | intergenic (+28/+4) | waaL/waaU | O‑antigen ligase/lipopolysaccharide core biosynthesis |
Reads supporting (aligned to +/- strand): ref base A (23/23); new base T (7/7); total (31/30) | |||||||||||
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.12e-01 |
TGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATT > NC_000913/3798116‑3798378 | tgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCAct > 1:303346/1‑128 (MQ=255) tgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCAct < 2:303346/128‑1 (MQ=255) gcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa > 4:166481/1‑103 (MQ=255) gcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa < 3:166481/103‑1 (MQ=255) aTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa < 1:816891/95‑1 (MQ=255) aTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa > 2:816891/1‑95 (MQ=255) ttATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActt < 4:870564/96‑1 (MQ=255) ttATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActt > 3:870564/1‑96 (MQ=255) tataTCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCa < 4:344223/92‑1 (MQ=255) tataTCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCa > 3:344223/1‑92 (MQ=255) ccTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa > 1:551601/1‑72 (MQ=255) ccTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa < 2:551601/72‑1 (MQ=255) ctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCa < 4:620309/74‑1 (MQ=255) ctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCa > 3:620309/1‑74 (MQ=255) ctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCt < 1:645950/84‑1 (MQ=255) ctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCt > 2:645950/1‑84 (MQ=255) tCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTGAATCTCTTTACTGAAtt < 1:328577/105‑1 (MQ=255) tCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAAAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAtt > 2:328577/1‑105 (MQ=255) cGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg > 1:428833/1‑148 (MQ=255) ttAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAct > 3:949909/1‑65 (MQ=255) ttAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAct < 4:949909/65‑1 (MQ=255) tacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTat < 2:428833/148‑1 (MQ=255) acaatTAATTAAGAATAAACAAGTTTAATAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTaa < 3:662087/88‑1 (MQ=255) acaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTACAACTCACTTCTTATCTATACAACTCAATCTCTTTACTGAATTAGTTaa > 4:662087/1‑88 (MQ=255) attaattaAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTaa < 1:1055290/64‑1 (MQ=255) attaattaAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTaa > 2:1055290/1‑64 (MQ=255) ttaattaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt > 1:340489/1‑70 (MQ=255) ttaattaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt < 2:340489/70‑1 (MQ=255) taaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGt > 1:505370/1‑128 (MQ=255) taaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGt < 2:505370/128‑1 (MQ=255) aGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGtt < 4:794232/111‑1 (MQ=255) aGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGtt > 3:794232/1‑111 (MQ=255) aTAAACAATTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATCGATTTATGATGATTCGGGGAcc > 1:966738/1‑148 (MQ=255) aTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTATTTAAAGTTTc < 3:652696/80‑1 (MQ=255) aTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc > 4:652696/1‑80 (MQ=255) agTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAt > 1:776972/1‑51 (MQ=25) agTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAt < 2:776972/51‑1 (MQ=25) agTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtat < 4:608648/87‑1 (MQ=255) agTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtat > 3:608648/1‑87 (MQ=255) agTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg < 2:448317/93‑1 (MQ=255) agTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg > 1:448317/1‑93 (MQ=255) agTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTAt < 3:724985/69‑1 (MQ=255) agTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTAt > 4:724985/1‑69 (MQ=255) agTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAggg < 2:966738/148‑1 (MQ=255) gTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTa < 4:491734/115‑1 (MQ=255) gTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTa > 3:491734/1‑115 (MQ=255) tGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAg < 1:1091581/52‑1 (MQ=37) tGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAg > 2:1091581/1‑52 (MQ=37) tGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGat > 2:790135/1‑148 (MQ=255) gAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGata > 1:745675/1‑148 (MQ=255) gTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATAtt < 1:790135/148‑1 (MQ=255) ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGa < 1:251775/66‑1 (MQ=255) ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGa > 2:251775/1‑66 (MQ=255) ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAtt > 4:105109/1‑98 (MQ=255) ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAtt < 3:105109/98‑1 (MQ=255) ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTa > 2:758811/1‑110 (MQ=255) ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTa < 1:758811/110‑1 (MQ=255) tAAAACTCACTTCTTCTCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCCGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGta < 1:556344/101‑1 (MQ=255) tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCtt > 3:108751/1‑71 (MQ=255) tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCtt < 4:108751/71‑1 (MQ=255) tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGtcg < 4:84245/86‑3 (MQ=255) tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg > 3:84245/1‑86 (MQ=255) tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGta > 2:556344/1‑101 (MQ=255) tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAg < 4:733782/104‑1 (MQ=255) tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAg > 3:733782/1‑104 (MQ=255) tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTa < 1:88020/131‑1 (MQ=255) tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTa > 2:88020/1‑131 (MQ=255) aaaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgg < 2:817857/86‑1 (MQ=255) aaaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgg > 1:817857/1‑86 (MQ=255) aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc < 2:267922/55‑1 (MQ=255) aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc > 1:267922/1‑55 (MQ=255) aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCtttt < 1:869403/71‑1 (MQ=255) aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCtttt > 2:869403/1‑71 (MQ=255) aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGt < 1:612813/76‑1 (MQ=255) aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGt > 2:612813/1‑76 (MQ=255) aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGAttt < 1:760564/106‑1 (MQ=255) aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGAttt > 2:760564/1‑106 (MQ=255) | TGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATT > NC_000913/3798116‑3798378 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |