Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,659,058 | 0 | T | C | 51.3% | 39.6 / 54.9 | 37 | D38G (GAT→GGT) | gadX | acid resistance transcriptional activator GadX |
Reads supporting (aligned to +/- strand): ref base T (15/3); new base C (0/19); total (15/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.71e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-08 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGC > NZ_CP009273/3658976‑3659140 | cATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaa > 2:62598/1‑90 (MQ=255) aTCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaa > 1:802312/1‑90 (MQ=255) aTCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaa < 1:176308/90‑1 (MQ=255) cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAc > 1:180305/1‑90 (MQ=255) cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAc > 1:461834/1‑90 (MQ=255) aaaTAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTc > 1:740435/1‑90 (MQ=255) aatAAAGGCGTCCCCTTTAACACAAATAAAAAATTTTCAACCCACTTATCTCCTCAATTTTGGTTTGCACCCGAAAAAACCAGGTCACCg < 1:116060/90‑1 (MQ=255) aaGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCGTCCGCAAAAACCAGGTCACCGCCAtt > 2:746141/1‑90 (MQ=255) aGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTa > 1:782975/1‑90 (MQ=255) aGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTa > 2:768575/1‑90 (MQ=255) cGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaa > 2:342037/1‑90 (MQ=255) cctttGACACAAAAAAAAATTTTTCAACACCTTTATCTACCAGATTTTGGCTTGCCCCCGAAAAACCCAGCTCCCCGCCATTAAAATAGc < 1:431819/87‑1 (MQ=255) cccTTGACACAAAAAAAAAATTTCCAACACACTTATCTCCTCAAATTTGGCTTCCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGc < 2:490756/90‑1 (MQ=255) gACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGatat > 2:391832/1‑90 (MQ=255) cccGAATAAAAAATTTCCAACCCCCTTTTCTCCTTTATTTTGCCTTCCTCCCGAAAAACCGGGGGCCCCGCCATTAAAATAGCGATATTc < 1:24559/88‑1 (MQ=255) cacaAATAAAAATTTTTCAACACACTTATCTCCTCAATTTTGGCTTCCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTc < 2:74517/90‑1 (MQ=255) acaaaTAAAAATTTTACAACACTCTTATCTCCTCGATTTTGGCTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCa < 2:499370/87‑1 (MQ=255) gAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAcc > 1:120632/1‑90 (MQ=255) gAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAcc > 1:281293/1‑90 (MQ=255) aaTTTTCACCACATTTATCTCCTCAATTTTGGCTTCCACCCAAAAAAACCAGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCa < 1:14803/90‑1 (MQ=255) tCAACCCACTTTTCTACTCAATTTTGGCTTCCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 1:677920/90‑1 (MQ=255) tCAACACCTTTATCTACTCGATTTTGGCTTGCACCCCAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 2:122829/90‑1 (MQ=255) tCAACACCCTTATCTCCTCAATTTTGGCTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 1:135243/90‑1 (MQ=255) tCAACACACTTATCTCCTCAAATTTGGCTTGCACCCAAAAAACCCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 2:661848/90‑1 (MQ=255) acAACCCACTTATCTACTCGAATTTGCCTTCCACCCGCAAAACCCGGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 2:562821/89‑1 (MQ=255) aaaaCCCACTTATCTCCACAAATTTGTTTTGCACCCAAAAAACCGGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 1:667353/88‑1 (MQ=255) cAACACATTTATCTCCTCAATTTTGGTTTGCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTga < 1:746141/90‑1 (MQ=255) acctttATCTCCTCGATTTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat < 2:281293/86‑1 (MQ=255) cTCGATTTGGGCTTGCACCCAAAAAACCCAGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTc < 2:461834/90‑1 (MQ=255) cTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTc > 1:172267/1‑90 (MQ=255) aTTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCa > 2:580055/1‑90 (MQ=255) tgtcTTGCACCCAAAAAAACCGGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc < 1:371221/87‑1 (MQ=255) tGGCTTCCACCCAAAAAACCCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc < 1:253563/90‑1 (MQ=255) tGGCTTCCACCCAAAAAACCCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc < 2:484393/90‑1 (MQ=255) ttttGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGc < 2:194648/88‑1 (MQ=255) gCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat < 2:43821/65‑1 (MQ=255) gCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat > 1:43821/1‑65 (MQ=255) | CATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGC > NZ_CP009273/3658976‑3659140 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |