| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 228,181 | 0 | C | T | 54.4% | 28.7 / 65.0 | 47 | F191F (TTC→TTT) | yafD | endonuclease/exonuclease/phosphatase family protein |
| Reads supporting (aligned to +/- strand): ref base C (1/20); new base T (25/0); total (27/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.74e-12 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.02e-09 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCG > NZ_CP009273/228094‑228270 | gggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGG‑GTTTTaa > 1:1520232/1‑90 (MQ=255)gggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTCaa < 1:15030/90‑1 (MQ=255)gggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTCaa < 2:775384/90‑1 (MQ=255) ggATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTGATTATGGCGGGGG‑GTTTGAATGcgcc > 1:971806/1‑87 (MQ=255) ggATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGG‑GTTTTTATGCCTg > 2:644822/1‑90 (MQ=255) gATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTCAATGCCTgg < 2:2105669/90‑1 (MQ=255) gTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑GTTTTAAAGCCTggggg > 1:783667/1‑87 (MQ=255) gTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGc < 2:1259743/90‑1 (MQ=255) gTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGc < 1:920430/90‑1 (MQ=255) gTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGc < 1:650950/90‑1 (MQ=255) tCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGG‑GTTTTAATGCCTgggggc > 1:837633/1‑86 (MQ=255) taGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGCCGTAg < 2:1763229/90‑1 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGG‑GTTTTAATTCCTGGGGGCGGAGAAgg > 1:633545/1‑90 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGG‑GTTTTAATGCCCTGGGGCGTAGAAgg > 1:859018/1‑90 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGG‑ATTTTAATGCCTGGGGGCGGAGAAgg > 1:297279/1‑90 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTTTATGCCTTGGGGCGTAGAAgg > 1:1484070/1‑90 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTTTATGCCCGGGGGCGGAGAAgg > 1:343864/1‑90 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTTAATGCCTGGGGGCGGAGAAgg > 1:382160/1‑90 (MQ=255) cAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGCCGTAGAAGGATg < 1:371801/90‑1 (MQ=255) gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGG‑GTTTTTATGCCCTGAGCCGGAGAAGGGTGaa > 1:243164/1‑90 (MQ=255) gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGG‑GTTTTAATGCCCTGGGCCGCAGAAGGAGGaa > 2:1089854/1‑90 (MQ=255) aCTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGCCGTAGAAGGATGAAcgc < 1:1524801/90‑1 (MQ=255) aCTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGCCGTAGAAGGATGAAcgc < 2:1633126/90‑1 (MQ=255) aCTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGCCGTAGAAGGATGAAcgc < 2:392542/90‑1 (MQ=255) ttCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGt < 1:1089854/90‑1 (MQ=255) ttCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGt < 1:1399469/90‑1 (MQ=255) cTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGG‑GTTTTAATACCCGGGGGCGGAGAAGGGGGAAAGCGGtat > 2:640312/1‑90 (MQ=255) cTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGG‑GATTTAATGCCCGGGGGCGGAGAAGGAGGAAagcggttt > 2:709028/1‑86 (MQ=255) aTTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGG‑AATTTTATGCCTTGGGGCGGAGGAGGGTGGAAGCGtttttt > 2:795976/1‑86 (MQ=255) tCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGG‑AGTTTAATGCCTGGGGGCGTAGAAGGGTGAAAGGGGTTTATTGCttttc > 1:2117670/1‑88 (MQ=255) tCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTgc < 1:1282084/90‑1 (MQ=255) caGCGGCCCGGTCATTATGGCGGGGG‑AATTTAATGCCTGGGGCCGGAGAAGGAAGAAAGCGTTATATTGCTTTTCGCGGGGAAAGGCGCt > 2:895695/1‑90 (MQ=255) caGCGGCCCGGTCATTATGGCGGGAG‑ATTTTAATGCCTGGGGGCGGAGAAGGGGGAAAGCGGTATATTGCTTTTCGCGGGGAAAGTCGCt > 2:683024/1‑90 (MQ=255) caGCGGCCCGGTCATTATGGCGGGAG‑ATTTTAAAGCCTGGAGCCGGAGAAGGGTGAAAGGGGTATATCGCTTTTCGCGGGGAATGTCGCt > 1:285008/1‑90 (MQ=255) gCGGCCCGGTCATTATGGCGGGGGCATTTC‑ATACCCGGGGGCGGAGAAGGGGGAAAAAGGTTTATCTCTTTGCGGGGGGAAAGTCGcggg > 2:2249196/1‑87 (MQ=255) gCGGCCCGGTCATTATGGCGGGGG‑GTTTTAATTCCTGGGGGCGGAGAAGGGAGAAAGCGGTATTTTGCTTTTCGCGGGGAAAGTGGCTgc > 2:376942/1‑90 (MQ=255) gCGGCCCGGTCATTATGGCGGGAG‑ATTTTAATGCCTGGGGCCGTAGAAGGATGAACGCGTTATATCGCTTTTCGCGGGGAAAGTCGCTgc > 1:596907/1‑90 (MQ=255) gCGGCCCGGTCATTATGGCGGGAG‑ATTTTAATGCCTGGGGCCGGAGAAGGGTGAACGCGTTATATCGCTTTTCGCGGGGAATGTCGCTgc > 1:769452/1‑90 (MQ=255) gCGGCCCGGTCATTATGGCGGGAG‑ATTTTAATGCCTGGAGCCGGAGAAGGATGAACGCGGTATATCGCTTTTCGCGGGGAAAGTCGCTgc > 1:786340/1‑90 (MQ=255) ggCCCGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTgcgc < 2:437043/90‑1 (MQ=255) gCCCGGTCATTATGGCGGGAG‑GTTTTAATGCCTGGGGCCGTAGAAAGAAGAACGCGGTATATTGCTTTTCGCGGGGAAAGTGGCTgcggc > 1:997718/1‑88 (MQ=255) cGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAgg < 2:1484070/90‑1 (MQ=255) cGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAgg < 2:2208367/90‑1 (MQ=255) cGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAgg < 2:1178280/90‑1 (MQ=255) ggCGGGGG‑GTTTTAATGCCTTGGGCCGTAGAAGGATGAAAGCGGTTTATTGCGTTTCGCGGGGAAAGTCGCTGCGCCCGGGGCGgttttc > 2:1674458/1‑88 (MQ=255) gggAG‑ATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCga < 1:566251/90‑1 (MQ=255) gggAG‑ATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCga < 2:1225143/90‑1 (MQ=255) cAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAgcg < 1:2041939/90‑1 (MQ=255) | GGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAG‑ATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCG > NZ_CP009273/228094‑228270 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |