Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,407,992 | 0 | A | T | 61.5% | ‑4.9 / 13.0 | 13 | Y2* (TAT→TAA) *65K (TAA→AAA) | rcbA ralR | double‑strand break reduction protein RcbA type I toxin‑antitoxin system endodeoxyribonuclease toxin RalR |
Reads supporting (aligned to +/- strand): ref base A (0/5); new base T (8/0); total (8/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.77e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAA‑TTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACC > NZ_CP009273/1407912‑1408073 | gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAA‑TTTTATACATGAac < 1:787933/90‑1 (MQ=255) ttCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAA‑TTTTATACATGAACAcc < 1:301212/90‑1 (MQ=255) gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAA‑TTTTTTTAATTAAAACCCCCCCTTTTTTTTTCAGTGGGTTGtttttttt > 1:85194/1‑88 (MQ=255) gTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAAATAAAAACC‑CCCTTTTTTTTTTCAGTGGGTTGTTTTTTTCCTCttttgtg > 1:134082/1‑86 (MQ=255) gTGCCACCTTCCTTTTCAATAGTGGCGGTAA‑TTTTTTTTATGAACACCCCCCTTTTTTTTTTCAGTGGTTTGTTTTTTTCATCtttttgg > 2:77873/1‑86 (MQ=255) gTGCCACCTTCCTTTTCAATAGTGGCGGTAA‑TTTTTTACATGAAAACCTCCCTTATTTTTTTCCGTGGGTCGTTTTTTTCATCtttttag > 2:315016/1‑86 (MQ=255) gTGCCACCTTCCTTTTCAATAGTGGCGGTAA‑TTTTTTAAAATAAAACCCCCCTTTTTTTTTTCCTTGGGTCGTTTTTTTCTTCtttttag > 1:507428/1‑86 (MQ=255) gTGCCACCTTCCTTTTCAATAGTGGCGGGAA‑ATTTTTTTATTAAAACCCCCCTTTTTTTTTTCATTGGGTTGGTTTTTTCATtttttttg > 1:559043/1‑86 (MQ=255) gCCACCTTCCTTTTCAATAGTGGCGGGGA‑ATTTTTAAATATAAACCCCCATTTTTTTTTTTATGGGGTGGTTTTTTTCATTTTTCTAGTg > 1:801235/1‑90 (MQ=255) tttCAATAGTGGCGGTAA‑TTTTTTACATGAACAcc > 2:104420/1‑35 (MQ=255) tttCAATAGTGGCGGTAA‑TTTTATACATGAACAcc < 1:104420/35‑1 (MQ=255) cAATAGTGGCGGTAA‑TTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACtt < 1:343212/90‑1 (MQ=255) gTAA‑TTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAAcc < 2:466886/90‑1 (MQ=255) | GATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAA‑TTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACC > NZ_CP009273/1407912‑1408073 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |