Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,659,062 | 0 | C | A | 60.0% | ‑1.2 / 23.3 | 21 | A37S (GCG→TCG) | gadX | acid resistance transcriptional activator GadX |
Reads supporting (aligned to +/- strand): ref base C (4/4); new base A (0/12); total (4/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.44e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGT‑CACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAC > NZ_CP009273/3658975‑3659138 | gCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGGCTTGCATCCGCaa < 2:100094/90‑1 (MQ=255) tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa > 1:567213/1‑90 (MQ=255) gggAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAcc < 1:545417/90‑1 (MQ=255) tGAAAGCGTGTCCCTTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGT‑CACCGcc > 2:225765/1‑90 (MQ=255) cTTGACACAAATAAAAATTTTTC‑AACCCACTTATCTACTCAATTTTGGCTTGCACCCGAAAAACCCAGGT‑CACCGCCATTAAAATAGCGa < 2:67669/90‑1 (MQ=255) acacAAAAAAAAAATTTAA‑AACACATTTTTCTCCTCAATTTTGCCTTGCACCCG‑AAAAACCCGGTCCACCGCCATTAAAATAGCGATAtt < 2:330299/90‑1 (MQ=255) ccaaaaTAAAAATTTTCC‑CACACCCTTGTCTCTAAGATTTTGGCTGGCACCCGCAAAACCCCGGA‑CGCCGCCATTAAAATAGCGATATTc < 2:638553/86‑1 (MQ=255) cacaAATAAAAATTTTCC‑AACACACTTATCTCCTCAAATTTGGCTTGCACCCGAAAAACCCAGGT‑CACCGCCATTAAAATAGCGATATTc < 2:98458/90‑1 (MQ=255) aaTACAAAATTTTC‑AACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACAAGGT‑CACCGCCATTAAAATAGCGATATTCACCa > 2:311314/1‑90 (MQ=255) taaaaaaTTTCC‑AACACACTTTTCTACTCGATTTTGGCTTGCACCCGAAAAACCCAGGT‑CACCGCCATTAAAATAGCGATATTCACCAtt < 1:40259/87‑1 (MQ=255) aaaaaaTTTTCAAACCCCTTTATCTCCTCAATTTTGGTTTGCACCCGAAAAACCCGGGC‑CCCCGCCATTAAAATAGCGATATTCACCAtt < 1:130586/88‑1 (MQ=255) ttC‑AACACATTTATCTACTGGATTTTGGTTTGCATCCGAAAAACCCGGGC‑CACCGCCATTAAAATAGCGATATTCACCATTAACCATGGt < 2:446439/90‑1 (MQ=255) tC‑AACACCCTTATCTACTCAATTTTGGCTTCCACCCAAAAAAACCAGGC‑CACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 1:637611/90‑1 (MQ=255) tC‑AACACACTTATTTTCTAGATTTTGGTTTGCACCCCAAAAACCCAGAT‑CACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 2:377033/90‑1 (MQ=255) cc‑AACACCTTTATCTACAAGATTTTGGCTTCCACCCAAAAAAACCAGGT‑CACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 1:330571/89‑1 (MQ=255) aaCCCACTTATCTACTAAATTTTGGTTTCCACCCAAAAAACCCAGGC‑CACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTgag < 1:400272/90‑1 (MQ=255) acatTTATCTACTCGATTTTGGCTTGCACCCGAAAAAACCAGGT‑CACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat < 2:633236/90‑1 (MQ=255) tACTCAATTTTGGCTTGCACCCAAAAAACCCAGGT‑CACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATg < 2:337643/90‑1 (MQ=255) tGGCTTGCACCCAAAAAACCCAGGT‑CACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc < 2:65507/90‑1 (MQ=255) ggCTTGCATCCGCAAAAACCAGGT‑CACCGCCATTAAAATAGCGACATTCACCATTAACCATGGTgag > 1:4209/1‑67 (MQ=255) ggCTTGCATCCGCAAAAACCAGGT‑CACCGCCATTAAAATAGCGACATTCACCATTAACCATGGTgag < 2:4209/67‑1 (MQ=255) | GCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGT‑CACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAC > NZ_CP009273/3658975‑3659138 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |