| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NZ_CP009273 | 1,200,135 | (A)5→4 | pseudogene (279/1137 nt) | BW25113_RS05965 → | terminase large subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,200,131 | 0 | A | . | 94.1% | 60.1 / ‑1.4 | 17 | pseudogene (275/1137 nt) | BW25113_RS05965 | terminase large subunit |
| Reads supporting (aligned to +/- strand): ref base A (0/1); new base . (10/6); total (10/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.12e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01 | |||||||||||
AACATTGAGGGGCCGTGCTACGACAAACGGCGGGAAGTCATCGAGATGCTCAACGGCTCGGTGCCAAACGATGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGC > NZ_CP009273/1199995‑1200260 | aacaTTGAGGGGCCGTGCTACGACAAACGGCGGGAAGTCATCGAGATGCTCAACGGCTCGGTGCCAAACGATGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGaaa < 1:177585/139‑1 (MQ=255) agtGAGGTGCCGTGCTACGACAAACGGCGGGAAGTCATCGAGATGCTCACCGGCTCGGTGCCAAACGATGAACTGTGCGGGATCATCTATACCGTGGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGaaaagc < 1:196997/137‑3 (MQ=255) gCCGTGCTACGACAAACGGCGGGAAGTCATCGAGATGCTCAACGGCTCGGTGCCAAACGATGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGG‑AAAAGCCAATCCa > 1:412056/1‑138 (MQ=255) tGCTCAACGGCTCGGTGCCAAACGATGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGtt < 1:288844/139‑1 (MQ=255) cTCGGTGCCAAACGATGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGG‑AAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTcagcagc > 1:163788/1‑139 (MQ=255) tGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGG‑AAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAAt > 2:187872/1‑138 (MQ=255) atcTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGG‑AAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCaa < 1:71048/138‑1 (MQ=255) tcTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGG‑AAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAAc > 1:27782/1‑139 (MQ=255) ccGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGG‑AAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCtt < 2:359732/138‑1 (MQ=255) ttGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGG‑AAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTaa > 2:389605/1‑138 (MQ=255) gaAGGTGACGACTGGACCGACCCGCAGGTGCTGG‑AAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGAGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTaaaacaa > 1:427402/1‑138 (MQ=255) aGGTGACGACTGGACCGACCCGCAGGTGCTGG‑AAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTaaaacaaaac < 1:419218/139‑1 (MQ=255) ggTGACGACTGGACCGACCCGCAGGTGCTGG‑AAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCG‑AAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAacac > 1:312053/1‑139 (MQ=255) acgacTGGACCGACCCGCAGGTGCTGG‑AAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCt < 1:389605/138‑1 (MQ=255) accCGCAGGTGCTGG‑AAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCg > 1:169579/1‑139 (MQ=255) accCGCAGGCGCTGG‑AAAAGCCAATCCAAATATTGGCGTGTCGGTTTATTGCGAATTTTTGTTGAGTCAGCAGCAGCGGGCGAAAAATAAGGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCACAATATATGGGCg > 2:427327/1‑139 (MQ=255) cAGGTGCTGG‑AAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGgc < 2:163788/139‑1 (MQ=255) gCTGG‑AAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAAc < 2:194034/96‑1 (MQ=255) gCTGG‑AAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAAc > 1:194034/1‑96 (MQ=255) | AACATTGAGGGGCCGTGCTACGACAAACGGCGGGAAGTCATCGAGATGCTCAACGGCTCGGTGCCAAACGATGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGC > NZ_CP009273/1199995‑1200260 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |