Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 653,132 | 0 | A | C | 53.6% | 32.8 / 41.9 | 28 | V88G (GTG→GGG) | crcB | fluoride efflux transporter CrcB |
Reads supporting (aligned to +/- strand): ref base A (7/6); new base C (0/15); total (7/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.45e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAA > NZ_CP009273/653044‑653213 | tGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAacac < 2:256839/90‑1 (MQ=255) gCCAGTGCGGTCATGGCAAAAGCCCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAAcacc < 2:802617/90‑1 (MQ=255) gCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAAcacc < 1:945865/90‑1 (MQ=255) gTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACtt < 1:64872/90‑1 (MQ=255) ggTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGc > 1:1183097/1‑90 (MQ=255) aaaGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGtt < 2:1193174/90‑1 (MQ=255) aCCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTgaga < 1:854987/90‑1 (MQ=255) aCCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTgaga < 1:816/90‑1 (MQ=255) ttGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGaa > 1:365997/1‑69 (MQ=255) ttGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGaa < 2:365997/69‑1 (MQ=255) gACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACCCCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGAcc < 2:385823/90‑1 (MQ=255) aaaCGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACCCCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCcac < 1:393961/90‑1 (MQ=255) aCGTTAAGTAATCCCAAGCCAAAGGGGCCCTCTTGAAAAAAAAACCCCCCTTCTGCCAAAAATGTTGAGAAGGTTGTTAGACCGCCACaa < 2:1192195/90‑1 (MQ=255) ttCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGaaaa < 1:511926/59‑1 (MQ=255) ttCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGaaaa > 2:511926/1‑59 (MQ=255) tAATGCCAAGCAAAAGCGGCCCTTTTGTAAAAAAACCCCCACTTCGGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCggt < 2:542423/90‑1 (MQ=255) acccccAGCAAAGGCGGCCCTTTTGAAACAAAAACCCCCCTTCTGCCAAAAATTTTGAGAGGGTTTTAAGCCCGCCACAAAATCCggtgg < 2:1247790/87‑1 (MQ=255) aTGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACCCCCCTTCTGCCAAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCggtgg < 1:1012066/90‑1 (MQ=255) aTGCCCAGCAAAAGGGGCCCTTTTGAAAAAAAAACCCCACTTCTGCCAAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCggtgg < 2:211151/90‑1 (MQ=255) aTGCCCAGCAAAAGCGGCCCTCTTGTAACAAAAACCCCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCggtgg < 1:859663/90‑1 (MQ=255) tGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCggtggt > 1:664213/1‑90 (MQ=255) gCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTaa > 1:502995/1‑90 (MQ=255) acggCGGCCCTCTTGTAACAAAAACCCCTCTTCTGCCAAAATTTTTGAGAGGGTTTTTAGACCGCCACAAAATCCGGTGGTGATTAATAc < 2:1007219/87‑1 (MQ=255) cGGCCCTTTTGAAAAAAAACCCCCCCTCCGGCCAAAAATGTTGAGAGGTTTGTTAGACCGCCACAAAATCCGGTCGTGATTAATACTTTc < 2:492405/90‑1 (MQ=255) gccccTCTTGAAAAAAAACCCCCCCTTCTGCCAAAAATTTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTcc < 2:502995/88‑1 (MQ=255) gCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCca > 1:1097762/1‑90 (MQ=255) gCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCca > 2:1126360/1‑90 (MQ=255) cttttGAAAAAAAACCCCCTCTCCGGCCAAAAATGTGGAGAGGGTTTTAAAACCGCCACAAAATCCGGTGGTGATTAATACTTTCCacac < 2:679749/87‑1 (MQ=255) tAAAAAAAACCCCCCTTCCGCCAAAAATTTTGAAAAGTTTTTTAGCCCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATc < 1:984842/90‑1 (MQ=255) aaaCAAAACCCCCCCTTCTGCAAAAAATGTTGAAAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATc < 2:823584/89‑1 (MQ=255) aCAAAAACCCCCCTTCTGCCAAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCaa < 1:720129/90‑1 (MQ=255) | TGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAA > NZ_CP009273/653044‑653213 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |