Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 553,473 | C→A | 7.4% | M156I (ATG→ATT) | lpxH ← | UDP‑2,3‑diacylglucosamine pyrophosphohydrolase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 553,473 | 0 | C | A | 7.4% | 185.9 / 3.0 | 53 | M156I (ATG→ATT) | lpxH | UDP‑2,3‑diacylglucosamine pyrophosphohydrolase |
Reads supporting (aligned to +/- strand): ref base C (24/25); new base A (2/2); total (26/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCAT > NC_000913/553329‑553600 | aaGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCAtt < 4:411960/148‑1 (MQ=255) gggTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTgc < 4:106070/144‑1 (MQ=255) gggTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTgc > 3:106070/1‑144 (MQ=255) tATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTg < 4:790627/140‑1 (MQ=255) tATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCGGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTg > 3:790627/1‑140 (MQ=255) cGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTc < 1:1063808/117‑1 (MQ=255) cGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTc > 2:1063808/1‑117 (MQ=255) tCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGtt > 1:789717/1‑126 (MQ=255) tCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGtt < 2:789717/126‑1 (MQ=255) gCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTcgcg < 1:696946/111‑1 (MQ=255) gCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTcgcg > 2:696946/1‑111 (MQ=255) ttGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCAt < 2:934082/102‑1 (MQ=255) ttGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCAt > 1:934082/1‑102 (MQ=255) tGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGt < 2:516660/117‑1 (MQ=255) tGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGt > 1:516660/1‑117 (MQ=255) tGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATa > 3:974630/1‑148 (MQ=255) tcACCTGACGTTTTACCATCGCACTGACCACCGCGTTTTGGTTAACGTCCCTGATCGCCAGTGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCTCCCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAg < 4:702384/147‑1 (MQ=255) ccTGATGTTATTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGGTTTACTGCTGTTGGCTTCTTTGCAGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCga > 3:702384/1‑147 (MQ=255) tGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGAATgcgcgc > 3:737430/1‑101 (MQ=255) tGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGAATgcgcgc < 4:737430/101‑1 (MQ=255) ttCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAg > 3:469128/1‑148 (MQ=255) ccATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAggg < 4:469128/148‑1 (MQ=255) aTCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGcc > 1:1086089/1‑142 (MQ=255) aTCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGcc < 2:1086089/142‑1 (MQ=255) aTCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGtt > 4:346671/1‑148 (MQ=255) cGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtg < 4:974630/148‑1 (MQ=255) cGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGAATGcgcgcga > 4:239246/1‑88 (MQ=255) cGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGAATGcgcgcga < 3:239246/89‑2 (MQ=255) gcgTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGgcgc > 2:1036672/1‑148 (MQ=255) tttGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGcc > 1:255948/1‑123 (MQ=255) tttGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGcc < 2:255948/123‑1 (MQ=255) ttGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGaaaa < 3:346671/148‑1 (MQ=255) ttAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCt > 4:264894/1‑148 (MQ=255) cATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGtt > 1:492300/1‑68 (MQ=255) cATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGtt < 2:492300/68‑1 (MQ=255) cGCCAGCGATTTACTTCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCa > 2:853084/1‑105 (MQ=255) cGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCa < 1:853084/105‑1 (MQ=255) cAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGt > 3:654021/1‑93 (MQ=255) cAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGt < 4:654021/93‑1 (MQ=255) aGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCgg < 1:1036672/148‑1 (MQ=255) tgctgTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtg > 1:630855/1‑98 (MQ=255) tgctgTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtg < 2:630855/98‑1 (MQ=255) ctgTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCgtgt < 3:264894/148‑1 (MQ=255) tGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAg > 1:827994/1‑61 (MQ=255) tGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAg < 2:827994/61‑1 (MQ=255) tcttTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAAt > 1:832333/1‑66 (MQ=255) tcttTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAAt < 2:832333/66‑1 (MQ=255) tcttTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAg > 1:936492/1‑106 (MQ=255) tcttTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAg < 2:936492/106‑1 (MQ=255) tGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCAt > 4:1122735/1‑138 (MQ=255) tGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCAt < 3:1122735/138‑1 (MQ=255) gTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGt > 3:673842/1‑63 (MQ=255) gTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGt < 4:673842/63‑1 (MQ=255) cATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCcgc > 4:27950/1‑102 (MQ=255) cATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCcgc < 3:27950/102‑1 (MQ=255) cATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGt < 3:624161/111‑1 (MQ=255) cATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGt > 4:624161/1‑111 (MQ=255) | AAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCAT > NC_000913/553329‑553600 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |