Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 536,014 | 0 | T | G | 54.2% | 17.6 / 28.3 | 24 | intergenic (+50/‑7) | allB/ybbY | allantoinase AllB/uracil/xanthine transporter |
Reads supporting (aligned to +/- strand): ref base T (3/8); new base G (13/0); total (16/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.24e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGT > NZ_CP009273/535927‑536103 | cgcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgtgt < 1:70928/90‑1 (MQ=255) gcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGGATTTTCCGGGTAAGgggtg > 2:505688/1‑90 (MQ=255) gAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGCTTTTCCGGGTTAGGGGTGttt > 1:512763/1‑90 (MQ=255) gAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGttt > 2:445877/1‑90 (MQ=255) tCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTc < 2:826095/90‑1 (MQ=255) aaaCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGGTAAGGGGTGTTTTTGTTCAATTTTTCAGTCa > 2:100649/1‑90 (MQ=255) aaaCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCa < 2:508851/90‑1 (MQ=255) aaaCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTTAGGGGGGTTTTTGTTCAATTTTTTAGTCa > 1:75261/1‑90 (MQ=255) agcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGa < 2:320675/90‑1 (MQ=255) agTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAg < 1:630656/90‑1 (MQ=255) aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGTGTTTTTTTTCCATTTTTCCGTCAGCCGCGAAAGGCt > 1:80746/1‑90 (MQ=255) ccTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAgg < 2:346823/90‑1 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCCGTTTAGGGGGGGTTTTGTTTAATTTTTCAGTCAGCCGCGGAAACCTGTTATTAGGa > 1:792725/1‑90 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGa > 2:347963/1‑90 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGa > 2:399246/1‑90 (MQ=255) gCCCGTCCTTGCGGCGGGCATTTTCCCGGTAAGGGGGGTTTTTTTTTAATTTTTCAGTCAGCCCCGGAAACCTGTTATTAGGGTTTTAGt > 1:1041951/1‑90 (MQ=255) cGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGGTTTTAGTGGt > 2:318622/1‑90 (MQ=255) cGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGt > 2:461830/1‑90 (MQ=255) cGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGGGTTTTTGTTTAATTTTGCAGTCAGCCGCGGAAAGCTGTTATCAGGGTTTTAGTGGt > 2:497328/1‑90 (MQ=255) tCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGttt < 2:977225/90‑1 (MQ=255) tCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGttt < 1:823057/90‑1 (MQ=255) gggCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTTAGTGGTtttttttttttt > 2:433622/1‑90 (MQ=255) gggCATTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGGTTTCAGTGGGTTTTTTTTAttt > 1:37282/1‑90 (MQ=255) gggCATTCTCCGGGTAAGGGGTGTTTTTGTTCAATTTTTCAGTCAGCCGCGAAAGCCTGTTATCAGGGTTTTAGTGGTtttttttttttt > 1:1063048/1‑90 (MQ=255) cGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac < 2:1084056/90‑1 (MQ=255) tAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGtt > 2:8223/1‑90 (MQ=255) tgtgtTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGt > 1:594556/1‑90 (MQ=255) | CGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGT > NZ_CP009273/535927‑536103 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |