Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,152,727 | 0 | T | G | 54.0% | 59.9 / 80.4 | 50 | G165G (GGT→GGG) | ycfH | metal‑dependent hydrolase |
Reads supporting (aligned to +/- strand): ref base T (0/23); new base G (27/0); total (27/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.26e-15 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.63e-11 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGT > NZ_CP009273/1152647‑1152816 | cTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTg < 1:2145665/90‑1 (MQ=255) cTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTg < 2:668225/90‑1 (MQ=255) cgcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAAGTAATGGAACTTggggtt > 1:1619130/1‑86 (MQ=255) cgcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATTGGTCTCGGGttt > 1:647264/1‑90 (MQ=255) cgcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATTGGTCTCGGGttt > 2:228186/1‑90 (MQ=255) aagaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACa < 1:1885882/90‑1 (MQ=255) gACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAAATACTTGATCTGGGTTTTTA‑ATTTCC‑CTTTcc > 1:1642792/1‑90 (MQ=255) ttGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTTTTTGGGCTCCGGTTTTTAATCTCC‑TTTTCCGTCAt > 2:2338672/1‑90 (MQ=255) gTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGa < 1:379171/90‑1 (MQ=255) ggCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGAcc < 2:677050/90‑1 (MQ=255) ggCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATTGATCTCGGGTTTTAAATCTCC‑TTTTTCTGCATTGTGGcc > 1:2446832/1‑90 (MQ=255) gCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCt < 2:1270262/90‑1 (MQ=255) gCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCt < 1:1019782/90‑1 (MQ=255) gTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTc < 1:2503051/90‑1 (MQ=255) gTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTc < 2:2986081/90‑1 (MQ=255) gTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGATCTCGGGTTTTACATCTCC‑TTTTTCGGCATTGTGACCTTc > 2:1110906/1‑90 (MQ=255) ctacACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGt < 1:1539215/90‑1 (MQ=255) ctacACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGt < 1:1507668/90‑1 (MQ=255) acTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTTATGGATCTCGGGTTTTTCATCTCC‑CTTTTCCGCCTTTTGAACTTCCCTTATg > 2:2380479/1‑90 (MQ=255) acTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATGGGTCTCGGGTTTTAAATTTCC‑TTTTTCCGCCTTTTGAACTTCCCTAAtt > 2:1363529/1‑89 (MQ=255) cTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGc < 2:1099961/90‑1 (MQ=255) cTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGc < 1:2882514/90‑1 (MQ=255) cTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGc < 1:1489310/90‑1 (MQ=255) cTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATATTTGGTCTCGGGTTTTTAAACTCC‑TTTTTCCGCATTTTGGCCTTCCCTAATGc > 2:2380291/1‑90 (MQ=255) cTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGATCTCGGGTTTTACATCTCC‑TTTTTCGGCATTGTGACCCTCCCTAAAGc > 2:1172737/1‑90 (MQ=255) gTTTTACAGAGGACAGAGAAACGGCGGGGGAAATACTTGGTCTCGGGTTTTTAATCTCC‑TTTTTCGGCGTTTTGAACTTTCGTAATGCgg > 1:537648/1‑90 (MQ=255) gTTTTACAGAGGACAGAGAAACGGCGGGGAAATTAATGGGTCTCGGGTTTTTCATCTCC‑TTTTTCCGCATTGTGGCCTTCCCTAATTCgg > 2:692381/1‑90 (MQ=255) gTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTTGGTTTCGGGTTTTTCATTTCC‑TTTTTCCGCCTTTTGGCCTTTCCTAATACgg > 2:2490360/1‑90 (MQ=255) gTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATTGATCTCTGGTTTTTAATTTCC‑CTTTTCTGCCGTTTGAGCTTCCCTAATTCgg > 1:2903211/1‑90 (MQ=255) gTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATGGATCTCGGGTTTTACATCTCC‑TTTTTCCGCACTGTGACCTTCCCGAATTCgg > 2:1662651/1‑90 (MQ=255) ttttACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGa < 1:2331644/90‑1 (MQ=255) tACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCa < 1:692381/90‑1 (MQ=255) cAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAAc < 2:1760220/90‑1 (MQ=255) gagGACAGAGAAACGGCGGGGAAATTACTGGGTCTCGGGTTTTACATCTCC‑TTTTTCGGCATTGTGGCCTTCCCTAATGCGGGGCAACTg > 2:1203016/1‑90 (MQ=255) gagGACAGAGAAACGGCGGGGAAAATACTGGATCTCGGGTTTTACATCTCC‑TTTTTCGGCATTTTGGCCTTCCCTAATGCGGGGCAAATg > 2:844388/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGGAATTACTGGGTCTCGGGTTTTAAATCTCC‑CTTTCCCGCCTTGTGGCCTTCCCTAAAGCGGGGCAAATgcg > 2:2921250/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGAAATTACTTGATCTCGGATTTTAAATCTCC‑TTTTTCGGCATTGTGACCTTCCCTAATGCGGGGGAACAgcg > 2:2000643/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGAAATTACTGGGTCTCGGGTTTTTCATCTCC‑TTTTTCGGCATTTTGACCTTTCCTAATGCGGGGGAAATgcg > 2:1147942/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGAAAATACTGGATCTCGGGTTTTTCATCTCC‑TTTTTCGGCCTTGTGGCCTTCCCTAATGCGGGGCAACTgcg > 2:2526876/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGAAAATACTGGATCTCGGGTTTTACATCTCC‑TTTTTCCGCCTTGTGGACTTCCCTAAATCGGGGCAACAgcg > 2:1902029/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGAAAATACTGGATCTCGGGTTTTAAATCTCC‑TTTTTCCGCATTTTGGCCTTCCCTAATTTGGGGGAACTgcg > 1:1033266/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGAAAATAATGGGTCTCGGGTTTTTTATCTCC‑TTTTTCCGCATTTTGACCTTCCCTAATGCGGGGGAACTgcg > 2:2331644/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGAAAATAATGGGTCTCGGGTTTTAAAACTCCTTTTTCCGGC‑TTGTGGCCTTTCCTAATGGGGGGGAACTtcg > 1:697458/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGAAAATAATGGATCTCGGGTTTTTTATCTCC‑TTTTTCCGCATTTTGAGCTTTCCTTATGTGGGGCAACTccg > 1:2962116/1‑90 (MQ=255) gACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTgcgc < 1:2526876/90‑1 (MQ=255) gACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTgcgc < 1:1147942/90‑1 (MQ=255) gaAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGAtgct < 2:237924/90‑1 (MQ=255) gcggGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGt < 1:2116724/90‑1 (MQ=255) gcggGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGt < 1:1363529/90‑1 (MQ=255) gcggGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGt < 2:2098342/90‑1 (MQ=255) tAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGt < 1:542290/90‑1 (MQ=255) tAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGt < 1:2729209/90‑1 (MQ=255) tAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGt < 2:559538/90‑1 (MQ=255) tAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGt < 2:647264/90‑1 (MQ=255) | CTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGT > NZ_CP009273/1152647‑1152816 |
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: ‑ |