| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | BS168_reference_genome | 1,338,342 | A→C | 15.3% | N1126T (AAC→ACC) | xkdO → | phage PBSX; putative tape measure protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | BS168_reference_genome | 1,338,342 | 0 | A | C | 15.3% | 78.5 / 12.0 | 39 | N1126T (AAC→ACC) | xkdO | phage PBSX; putative tape measure protein |
| Reads supporting (aligned to +/- strand): ref base A (20/13); new base C (3/3); total (23/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.74e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01 | |||||||||||
CTGCTTATTCCGGAATGGGTGCAGTAGGGGGAAATGTGAAACAGTGGCTTATGGCAGCAATAATGGCCACAAAGACACCATTGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGCTCATAAAGCACCGGGCATGGGTAACATTAGAAATCCGATTCA > BS168_reference_genome/1338210‑1338474 | ctgctTATTCCGGAATGGGTGCAGTAGGGGGAAATGGGAAACAGTGGCTTATGGCAGCAATAATGGCCACAAAGACACCATTGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGaa < 5:153763/139‑1 (MQ=255) ctTATTCCGGAATGGGTGCAGTAGGGGGAAATGTGAAACAGTGGCTTATGGCAGCAATAATGGCCACAAAGACACCATTGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCACCCCGAAc < 5:213009/137‑1 (MQ=255) ccGGAATGGGTGCAGTAGGGGGAAATGTGAAACAGTGGCTTATGGCAGCAATAATGGCCACAAAGACACCATTGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGAACTCTATTa < 4:88697/138‑1 (MQ=255) ggAATGGGTGCAGTAGGGTGAAATGTGAAACAGTGGCTTATGGCAGCAATAATGGCCACAAAGACACCATTGAGCTGGCGCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGAACTCTATTAAcc < 3:186684/139‑1 (MQ=255) ggAATGGGTGCAGTAGGGGGAAATGTGAAACAGTGGCTTATGGCAGCAATAATGGCCACAAAGACACCATTGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGAACTCTATTAAcc < 5:297243/139‑1 (MQ=255) aaTGGGTGCAGTAGGGGGAAATGTGAAACAGTGGCTTATGGCAGCAATAATGGCCACAAAGACACCATTGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGAActc < 1:2561/129‑1 (MQ=255) aaTGGGTGCAGTAGGGGGAAATGTGAAACAGTGGCTTATGGCAGCAATAATGGCCACAAAGACACCATTGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGAActc > 2:2561/1‑129 (MQ=255) aaTGGGTGCAGTAGGGGGAAATGTGAAACAGTGGCTTATGGCAGCAATAATGGCCACAAAGACACCATTGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGAACTCTATTAACCtg > 8:273663/1‑139 (MQ=255) aaaTGTGAAACAGTGGCTTATGGCAGCAATAATGGCCACAAAGACACCATTGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGaa < 1:38263/108‑1 (MQ=255) aaaTGTGAAACAGTGGCTTATGGCAGCAATAATGGCCACAAAGACACCATTGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGaa > 2:38263/1‑108 (MQ=255) tGGCTTATGGCAGCAATAATGGCCACAAAGACACCATTGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGaa < 2:58993/95‑1 (MQ=255) tGGCTTATGGCAGCAATAATGGCCACAAAGACACCATTGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGaa > 1:58993/1‑95 (MQ=255) cagcaATAATGGCCACAAAGACACCATTGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCg > 6:51852/1‑83 (MQ=255) cagcaATAATGGCCACAAAGACACCATTGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCg < 5:51852/83‑1 (MQ=255) cacCATTGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCACCCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACca > 8:49788/1‑138 (MQ=255) ttGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACAttt > 4:20370/1‑138 (MQ=255) ttGAGCTGGCTATCGGGATTAATGACGATTGCACAGTATGAAACAGGGGGCAACCCGAACTCTATTAACCTGTCGGATAGCAACGCTAAGGCCGCAAATCCAACACAAGCACTCATGAAGACAGTACCAACCACAtta > 6:175278/1‑137 (MQ=255) aGCTGGCTCGCGGGAATAATGACGATTTCACAGTATGAATCAGGGCGCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAAt < 8:166388/138‑1 (MQ=255) gCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGAACTCTATTAACCTGTGGGATCGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGCTCa > 5:224176/1‑138 (MQ=255) gCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGa > 1:1246/1‑113 (MQ=255) gCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGa < 2:1246/113‑1 (MQ=255) gCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGc < 8:104064/135‑1 (MQ=255) gCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGc > 7:104064/1‑135 (MQ=255) cGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGACCTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGCTCATAAAg > 5:247156/1‑138 (MQ=255) ggATTAATGACGATTGCACAGTATGAATCAGGGGGCACCCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGCTCATAAAGCAc < 3:240432/139‑1 (MQ=255) gATTGCACAGTATGAATCAGGGGGCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGCTCATAAAGCACCGGGCATggg > 4:257362/1‑138 (MQ=255) tATGAATCAGGGGGCACCCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGCTCATAAAGCACCGGGCATGGGTAACATTAGa > 3:220779/1‑138 (MQ=255) tATGAATCAGGGGGCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGc > 5:266337/1‑107 (MQ=255) tATGAATCAGGGGGCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGc < 6:266337/107‑1 (MQ=255) tATGAATCAGGGGGCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGCTCATAAAGCACCGGGCATGGGTAACATTAGaa < 7:273663/139‑1 (MQ=255) aTCAGGGGGCAACCCGAACTCTATTAACATGTGGGATAGCAACGCGAAGGCGGGATATCCATCACAAGGACTCATGCAGACAGTTACATCCACATTAAATGCTCATAAAGCACCGGGCATGGGTAACATTTGAAATcc > 7:84873/1‑138 (MQ=255) ggggCACCCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATcc > 5:231682/1‑55 (MQ=255) ggggCACCCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATcc < 6:231682/55‑1 (MQ=255) ggggCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGCTCATAAAGCACCGGGCATGGGTAACATTAGAAATCCGATTca > 6:289594/1‑139 (MQ=255) ggggCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGCTCATAAAGCACCGGGCATGGGTAACATTAGAAATCCGATTca > 6:261958/1‑139 (MQ=255) ggggCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGCTCATAAAGCACCGGGCATGGGTAACATTAGAAATCCGATTca > 2:255491/1‑139 (MQ=255) ggggCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGCTCATAAAGCACCGGGCATGGGTAACATTAGAAATCCGATTca > 5:205202/1‑139 (MQ=255) ggggCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGCTCATAAAGCAACGGGCATGGGTAACATTAGAAATCCGATTc > 6:150223/1‑138 (MQ=255) ggcgCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGCTCATAAAGCACCGGGCATGGGTAACATTAGAAATCCGATTca > 1:208456/4‑139 (MQ=255) | CTGCTTATTCCGGAATGGGTGCAGTAGGGGGAAATGTGAAACAGTGGCTTATGGCAGCAATAATGGCCACAAAGACACCATTGAGCTGGCTCTCGGGATTAATGACGATTGCACAGTATGAATCAGGGGGCAACCCGAACTCTATTAACCTGTGGGATAGCAACGCGAAGGCGGGAAATCCATCACAAGGACTCATGCAGACAGTTCCAACCACATTTAATGCTCATAAAGCACCGGGCATGGGTAACATTAGAAATCCGATTCA > BS168_reference_genome/1338210‑1338474 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |