Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,029,116 | 0 | A | C | 60.0% | 7.8 / 16.2 | 15 | G175G (GGT→GGG) | prfB | peptide chain release factor 2 |
Reads supporting (aligned to +/- strand): ref base A (6/0); new base C (0/9); total (6/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.00e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCT > NZ_CP009273/3029039‑3029205 | ggggggTGACCGCGGTTTTCTGTCCGCACCCAGCCGTAAGCGTAATCGCCGGAGTTTTTGATCGTCACGGATTTAATCCCCGCCACTTCa < 2:198944/86‑1 (MQ=255) aGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCa > 2:459783/1‑90 (MQ=255) tGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTcc > 2:592479/1‑90 (MQ=255) tttCTGACCCCACCCAGCCGAAAGGGTACTCGCGGGATTTTTTGATCGCCACGTTTTTAATCCCCGCCACTTCACCTTCCGACTCTTCga < 2:337436/90‑1 (MQ=255) tCTGTCCGCACCCAGCCGAAAGCGTAATCGCGGGAGTTTTTGATCGTCAGGGTTTTAATCCCCGCCACTTCACCTTCCGACTCTTCgatg < 2:19582/90‑1 (MQ=255) tAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCa > 1:207760/1‑90 (MQ=255) tAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCa > 2:350304/1‑90 (MQ=255) gCGTAACCGCGGGAGTTTTTGATCGTCAGGTTTTTAACCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCAcgc < 1:489409/90‑1 (MQ=255) gTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGa > 1:49047/1‑90 (MQ=255) gTAACCGCGGGATTTTTTGATCGCCAGGTTTTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGa < 2:49047/90‑1 (MQ=255) aTCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTc > 1:86279/1‑90 (MQ=255) attttttGATCGTCAGTTTTTTAACCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAgc < 1:345362/87‑1 (MQ=255) tttttGACCGTCCCGTTTTTAACCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAgcgc < 2:520146/89‑1 (MQ=255) tttGATCGTCACGGATTTAATCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAg < 1:532738/90‑1 (MQ=255) tCACGAATTAAACCCCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGatacatac < 2:350387/90‑1 (MQ=255) aCCCGCCCCTTCACCTTCCGAATCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCt < 2:176846/90‑1 (MQ=255) | AGGCGGTGAACGCCGGTTTCTGTACGCAGCCAGCCGTAAGCGTAATCGCCGGAGATTTTGATCGTCACGGATTTAATACCCGCCACTTCACCTTCCGACTCTTCGATGATTTCAGTTTTGAAACCACGCGATTCTGCCCAGCGCAGATACATACGCTCAAGCATGCT > NZ_CP009273/3029039‑3029205 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |