Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,029,116 | 0 | A | C | 57.1% | 6.8 / 16.4 | 14 | G175G (GGT→GGG) | prfB | peptide chain release factor 2 |
Reads supporting (aligned to +/- strand): ref base A (6/0); new base C (0/8); total (6/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.33e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGC > NZ_CP009273/3029037‑3029200 | ccAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACtt > 1:48068/1‑90 (MQ=255) aGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCa > 1:317497/1‑90 (MQ=255) aGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCAcg > 2:726757/1‑90 (MQ=255) gCGTACGCGCGGGATTTTTTTATCGTCAGGTTTTTAACCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCAcgc < 1:499049/90‑1 (MQ=255) gCGAAATCGCGGGAGTTTTTGCCCGTCACGGTTTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCAcgc < 1:165691/90‑1 (MQ=255) tAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGAt > 1:760623/1‑90 (MQ=255) tAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGAt > 1:795753/1‑90 (MQ=255) aaTCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGAtt > 2:645864/1‑90 (MQ=255) cGGAGTTTTTGATCGTCACGGTTTTAACCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGccc < 1:630595/90‑1 (MQ=255) agtttttGACCGTCAGGTTTTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAgc < 2:760623/87‑1 (MQ=255) tttttGATCGTCAGGGATTAAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAgcgc < 2:48068/89‑1 (MQ=255) aTTTTGATCGCCGCGAATTAAACCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAgcgc < 2:67351/90‑1 (MQ=255) cGGATTAAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCt < 2:175342/90‑1 (MQ=255) taaaTCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGc < 2:611284/88‑1 (MQ=255) | CCAGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGC > NZ_CP009273/3029037‑3029200 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |