Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 653,356 | 0 | A | C | 57.9% | 17.2 / 21.4 | 19 | G13G (GGT→GGG) | crcB | fluoride efflux transporter CrcB |
Reads supporting (aligned to +/- strand): ref base A (8/0); new base C (0/11); total (8/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.32e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAAC > NZ_CP009273/653267‑653437 | aTCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTa > 1:449268/1‑90 (MQ=255) cgcccccAACGGAATCGCCTGATGCAGGGGTTTAAATCGCATACTTAACAGCCATCTCGCCCCCCTCCCCGTCCCACCGCCAATAAAAAc < 2:310463/87‑1 (MQ=255) ccccAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGc > 1:615255/1‑90 (MQ=255) ccccAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGGACCACCGCCAATAAAAACTGc > 1:580253/1‑90 (MQ=255) gAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAaagaa > 2:68571/1‑90 (MQ=255) gAATCGCCTGATGCAGTGGGTAAAACCGCATATTTAACAGCCATCTCGCCACGCTCCCCGTCCCACCGCCAATAAAAACTGCTAAaagaa < 2:583152/90‑1 (MQ=255) aaaTCCCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCCCCCCGCTTCCCGTCCCACCGCCAATAAAAACTGCTAAaagaa < 2:214250/89‑1 (MQ=255) cTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTaa > 2:561007/1‑90 (MQ=255) gTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCaa > 2:335115/1‑90 (MQ=255) ccATACTAAACAGCCATCTCCCCACGCTTCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTatt < 2:449268/89‑1 (MQ=255) cacccATCTCGCCCCGCTTCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTAt < 2:573549/87‑1 (MQ=255) cacccACCTCCCCCCGCTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTAt < 1:68571/87‑1 (MQ=255) gCCATCCCGCCCCTCTTCCCGTCCCCCCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGa < 1:284921/90‑1 (MQ=255) aTCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTg > 1:427218/1‑90 (MQ=255) aTCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTg > 2:598273/1‑90 (MQ=255) aTCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTg > 2:760661/1‑90 (MQ=255) cccccACCCTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATa < 1:129362/87‑1 (MQ=255) ccGCCACGCTTCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATa < 2:745585/89‑1 (MQ=255) ccGCCACGCTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATa < 2:21646/89‑1 (MQ=255) ctccccGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAAc < 2:13975/87‑1 (MQ=255) | ATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAAC > NZ_CP009273/653267‑653437 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |