| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,029,116 | 0 | A | C | 52.3% | 44.4 / 67.7 | 44 | G175G (GGT→GGG) | prfB | peptide chain release factor 2 |
| Reads supporting (aligned to +/- strand): ref base A (20/1); new base C (0/23); total (20/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.19e-11 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.33e-08 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAA‑TACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCT > NZ_CP009273/3029035‑3029205 | cACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAA‑TACCCGCCAc > 1:110068/1‑90 (MQ=255) ccAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAA‑TACCCGCCACtt > 2:616448/1‑90 (MQ=255) aGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAA‑TACCCGCCACTTCa > 1:1706236/1‑90 (MQ=255) aGGCGGTGAACGCCGGGTTCTGTACGGAGCAAGCCGTAAGCGTAATTTCCGGAGGTTTTGATCGTCACGGATTTAA‑TACACGCCACTGCa > 1:525505/1‑90 (MQ=255) ggTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAA‑TACCCGCCACTTCACCtt > 1:344335/1‑90 (MQ=255) aaCGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGTTTTAA‑TACCCGCCACTTCACCTTCCGa > 2:807947/1‑90 (MQ=255) aCGCCGTTTTCTGTACGCAGCCACCCGAAAGCGTAATCGCCGGATTTTTTGATCGTCAGGGATTTAA‑TCCCCGCCACTTCACCTTCCGAc < 2:285565/90‑1 (MQ=255) cGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAA‑TACCCGCCACTTCACCTTCCGAct > 1:539471/1‑90 (MQ=255) gtttttCTGCCCCCACCCAGCGGAAAGCGTACCCGCGGGAGTTTTTGATCGCCAGGGATTAAA‑TCCCCGCCACTTCACCTTCCGACTCtt < 1:567322/87‑1 (MQ=255) gTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAA‑TACCCGCCACTTCACCTTCCGACTCTTCg > 2:1010282/1‑90 (MQ=255) ttttTGTCCGCACCCAGCCGTAAGCGTAATCGCGGGATTTTTTGATCGTCACGGTTTTAA‑TCCCCGCCACTTCACCTTCCGACTCTTCga < 2:344335/90‑1 (MQ=255) tGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAA‑TACCCGCCACTTCACCTTCCGACTCTTCgatgat > 1:1564507/1‑90 (MQ=255) gCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAA‑TACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAg > 1:1263989/1‑90 (MQ=255) ttttAATCGCGGGATTTTTTGACCGCCACGTTTTTAA‑TCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCAcgc < 2:1426449/87‑1 (MQ=255) gCGTACTCGCGGGATTTTTTGATCGTCCCGTTTTTAA‑TCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCAcgc < 2:1098507/90‑1 (MQ=255) gCGTACCCGCGGGATTTTTTGTTCGTCAGGTTTTTAA‑TCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCAcgc < 1:616448/90‑1 (MQ=255) gCGTACCCGCGGGATTTTTTGACCGTCACGTATTTAA‑CCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCAcgc < 1:166393/90‑1 (MQ=255) gCGTACCCGCGGGAGTTTTTTCTCGTCAGGTTTTTAA‑CCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCAcgc < 1:1459158/90‑1 (MQ=255) gCGTAACCGCGGGATTTTTTGATCGTCACGTTTTTAA‑TCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCAcgc < 2:1128067/90‑1 (MQ=255) tAATCGCCGGAGATTTTGATCGTCACGGATTTAA‑TACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGAt > 1:1426449/1‑90 (MQ=255) tAATCGCCGGAGATTTTGATCGTCACGGATTTAA‑TACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGAt > 1:1344209/1‑90 (MQ=255) tAATCGCCGGAGATTTTGATCGTCACGGATTTAA‑TACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGAt > 2:84689/1‑90 (MQ=255) ctCGCGGGAGTTTTTGTTCGTCAGGTTTTAAA‑CCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTc < 1:1423947/89‑1 (MQ=255) gCCGGAGATTTTGATCGTCACGGATTTAA‑TACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGc > 1:115212/1‑90 (MQ=255) gggAGTTTTTGATCGCCACGGATTTAA‑TCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGccc < 2:140047/89‑1 (MQ=255) gggAGTTTTTGATCGCCACGATTTTAA‑TCCCCGCCACTTCGCCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGccc < 2:496091/89‑1 (MQ=255) cGGAGATTTTGATCGTCACGGATTTAA‑TACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGccc > 1:403057/1‑90 (MQ=255) aattttttGTTCGTCA‑GGATTTAACCCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAg < 1:896350/86‑1 (MQ=255) tttttGCTCGTACCGTATTTAA‑CCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAgcgc < 1:23022/89‑1 (MQ=255) tttttGATCTTCAGTTTTTTAA‑CCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAgcgc < 1:1163609/89‑1 (MQ=255) tttttGATCTTCAGGGTTTTAA‑TCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAgcgc < 2:347832/89‑1 (MQ=255) tttttGATCGTCACGATTTAAA‑TCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAgcgc < 2:716499/89‑1 (MQ=255) tttttAATCGTCAGGGTTTAAA‑TCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAgcgc < 2:372076/89‑1 (MQ=255) ttttgGATCGTCACGTTTTTAA‑CCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAgcgc < 1:886709/89‑1 (MQ=255) ttttGATCGTCACGAATTAAA‑CCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCa < 1:282417/90‑1 (MQ=255) tttGATCGCCACGGTTTTAA‑TCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAg < 2:122466/90‑1 (MQ=255) tGATCGCCACGTATTTAA‑CCCCCGCCACTTCACCTTCCGACTCTTCGATGTTTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGat < 2:1442275/90‑1 (MQ=255) gTCACGGATTTAA‑TACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGatacata > 2:977995/1‑90 (MQ=255) gTCACAGATTTAA‑TACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGatacata > 2:77705/1‑90 (MQ=255) ttAA‑TCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGc < 2:290899/90‑1 (MQ=255) aa‑TACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCAcgc > 1:666903/1‑55 (MQ=255) aa‑TACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCAcgc < 2:666903/55‑1 (MQ=255) aa‑TACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCAt > 2:170292/1‑90 (MQ=255) aa‑TACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCAt > 1:1526958/1‑90 (MQ=255) aCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCt < 2:1344209/90‑1 (MQ=255) | CACCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAA‑TACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCT > NZ_CP009273/3029035‑3029205 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |