Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,659,062 | 0 | C | A | 63.6% | 1.8 / 26.0 | 22 | A37S (GCG→TCG) | gadX | acid resistance transcriptional activator GadX |
Reads supporting (aligned to +/- strand): ref base C (7/1); new base A (0/14); total (7/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.69e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTC‑GAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAC > NZ_CP009273/3658977‑3659138 | aTCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTC‑GAATTTGGCTTGCATCCGCaaaa > 2:529427/1‑90 (MQ=255) tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTC‑GAATTTGGCTTGCATCCGCaaaaa > 2:530486/1‑90 (MQ=255) aaaTAATGAAAGCGTGTCCCTTGACACAAATACAAAATTTTCAACACACTTATCTCCTC‑GATTTTGGCTTGCACCCGAAAAACCCAGGTc < 1:215907/90‑1 (MQ=255) gtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTC‑GAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATa > 1:115640/1‑90 (MQ=255) gtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTC‑GAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATa > 2:104164/1‑90 (MQ=255) cccTTAACACAAAAAAAAATTTTCCAACACCCTTATCTACTC‑AATTTTGGCTTCCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGc < 1:185741/90‑1 (MQ=255) ctggACACAAAAAAAAATTTTTCAACCCACTTATCT‑CTCAAAATTTGCCTTCCACCCGCAAAAACCAGCCCACCGCCATTAAAATAGCGa < 1:394584/87‑1 (MQ=255) cTTGACACAAAAAAAAATTTTTAAACCCATTTATCTCCTC‑GATTTTGGCTTCCACCCGAAAAACCCAGGCCACCGCCATTAAAATAGCGa < 2:439435/90‑1 (MQ=255) acacAAAAAAAAAATTTTAAACACCTTTATCTCCTC‑AATTTTGGCTTCCACCCAAAAAAACCAGCTCACCGCCATTAAAATAGCGATAtt < 2:540524/90‑1 (MQ=255) cacaAATAAAAATTTTTAAACCCCCTTATCTCCTC‑GATTTTGGTTTGCACCCAAAAAACCCAGGTCCCCGCCATTAAAATAGCGATATTc < 1:104164/90‑1 (MQ=255) cacaAATAAAAAATTTTCAACACATTTATCTACTC‑AATTTTGGCTTCCACCCGAAAAACCCAGGCCACCGCCATTAAAATAGCGATATTc < 2:490841/90‑1 (MQ=255) gAATACAAAATTTTCAACACACTTATCTACTC‑GAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAcc > 1:564961/1‑90 (MQ=255) aaTTTTCAACCCATTTATCTACTC‑AATTTTGGTTTGCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCa < 1:68541/90‑1 (MQ=255) tCAACACTTTTATCTACTC‑GATTTTGGCTTGCACCCGAAAAACCCAGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 2:115640/90‑1 (MQ=255) tCAACACCTTTATCTACTC‑AAATTTGGTTTGCACCCAAAAAAACCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 2:564961/90‑1 (MQ=255) ccAACACCTTTCTCTACAA‑GATTTTGCCTTGCACCCAAAAAACCCGGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 1:515565/89‑1 (MQ=255) ccacTTATCTCCTC‑AATTTTGGTTTCCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat < 1:157619/89‑1 (MQ=255) ctttATCTACTC‑AATTTTGGCTTCCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAatat < 1:129485/88‑1 (MQ=255) cTTATCTACTC‑GAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAatata > 1:256208/1‑90 (MQ=255) tATCTCCTC‑GATTTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAtt < 2:250669/90‑1 (MQ=255) c‑GAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCtt > 2:426860/1‑90 (MQ=255) tGGCTTACACCCCAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc < 2:154514/90‑1 (MQ=255) | ATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTC‑GAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAC > NZ_CP009273/3658977‑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: ‑ |