Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,529,938 | G→A | pseudogene (2037/2037 nt) R6K (AGA→AAA) |
rhsE → ydcD → |
pseudogene, Rhs family putative immunity protein for RhsE |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,529,938 | 0 | G | A | 98.5% | 206.2 / NA | 65 | pseudogene (2037/2037 nt) R6K (AGA→AAA) | rhsE ydcD | pseudogene, Rhs family putative immunity protein for RhsE |
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (29/35); minor base T (0/1); total (29/36) | |||||||||||
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 |
CCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTT > NC_000913/1529817‑1530053 | ccAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCAttt < 2:15944/139‑1 (MQ=255) ccAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCAttt < 1:1036666/139‑1 (MQ=255) tGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCt < 2:968698/139‑1 (MQ=255) aTGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTg < 1:1925870/139‑1 (MQ=255) tGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTAtt > 2:326298/1‑139 (MQ=255) aaaaaTATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAAt < 2:318991/122‑1 (MQ=255) aaaaaTATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAAt > 1:318991/1‑122 (MQ=255) tataACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTgttgtt > 2:1160806/1‑139 (MQ=255) tGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTg < 2:1196670/138‑1 (MQ=255) tGGGAAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTg > 1:1196670/1‑138 (MQ=255) ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCAtttt < 1:298595/101‑1 (MQ=255) ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCAtttt > 2:298595/1‑101 (MQ=255) ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAAt < 1:1604493/107‑1 (MQ=255) ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAAt > 2:1604493/1‑107 (MQ=255) ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGt > 1:1979304/1‑139 (MQ=255) ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGt > 2:858737/1‑139 (MQ=255) ctctTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGAt > 1:865086/1‑139 (MQ=255) cGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGc < 1:1944880/77‑1 (MQ=255) cGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGc > 2:1944880/1‑77 (MQ=255) cGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACaa > 2:298078/1‑139 (MQ=255) gTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTg < 2:814789/139‑1 (MQ=255) aaaGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATc < 2:1979304/139‑1 (MQ=255) aaGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCg > 2:306840/1‑139 (MQ=255) aaGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCg > 2:570329/1‑139 (MQ=255) gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTAcc < 2:1638632/79‑1 (MQ=255) gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTAcc > 1:1638632/1‑79 (MQ=255) gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTg > 2:1550616/1‑84 (MQ=255) gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTg < 1:1550616/84‑1 (MQ=255) gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTgttgtt > 1:327717/1‑97 (MQ=255) gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTgttgtt < 2:327717/97‑1 (MQ=255) gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTg < 2:877421/104‑1 (MQ=255) gATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTg > 1:877421/1‑104 (MQ=255) aaTAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACa < 1:298078/139‑1 (MQ=255) gcaaagcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTg > 2:15135/1‑89 (MQ=255) gcaaagcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTg < 1:15135/89‑1 (MQ=255) gcaaagcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATa < 2:568806/134‑1 (MQ=255) gcaaagcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATa > 1:568806/1‑134 (MQ=255) gcaaagcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAaca > 2:514090/1‑139 (MQ=255) aagcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTAcc > 1:768681/1‑83 (MQ=255) aagcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTAcc < 2:768681/83‑1 (MQ=255) agcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCt > 1:312644/1‑83 (MQ=255) agcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCt < 2:312644/83‑1 (MQ=255) agcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGAt < 2:664862/93‑1 (MQ=255) agcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGAt > 1:664862/1‑93 (MQ=255) gcaaGGGAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCaaa < 1:1160806/139‑1 (MQ=255) ggAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCg < 2:1116531/139‑1 (MQ=255) gAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGt < 1:326298/139‑1 (MQ=255) gAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGt < 1:514090/139‑1 (MQ=255) gAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGt < 2:1898618/139‑1 (MQ=255) gAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGt < 2:865086/139‑1 (MQ=255) gAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGt < 1:570329/139‑1 (MQ=255) gAGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGt < 1:306840/139‑1 (MQ=255) aGTTCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGtt < 1:858737/139‑1 (MQ=255) ttCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAAtt < 2:1643074/90‑1 (MQ=255) ttCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAAtt > 1:1643074/1‑90 (MQ=255) ttCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTAtttt < 1:864549/111‑1 (MQ=255) ttCAGATGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTAtttt > 2:864549/1‑111 (MQ=255) aTGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGa > 2:1194396/1‑111 (MQ=255) aTGCTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCCTATTTTATTGa < 1:1194396/111‑1 (MQ=255) cTTCAAATTATTATAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGt < 2:433650/128‑1 (MQ=255) cTTCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGt > 1:433650/1‑128 (MQ=255) tCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTgttgtt < 1:201255/58‑1 (MQ=255) tCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTgttgtt > 2:201255/1‑58 (MQ=255) tCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGt < 2:69690/126‑1 (MQ=255) tCAAATTATTAAAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGt > 1:69690/1‑126 (MQ=255) | CCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTT > NC_000913/1529817‑1530053 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |