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% | 7.5 / 30.6 | 20 | A37S (GCG→TCG) | gadX | acid resistance transcriptional activator GadX |
Reads supporting (aligned to +/- strand): ref base C (8/0); new base A (0/12); total (8/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.94e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAC > NZ_CP009273/3658976‑3659138 | cATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaa > 2:631897/1‑90 (MQ=255) tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa > 2:349214/1‑90 (MQ=255) cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAc > 1:265870/1‑90 (MQ=255) cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAc > 1:588231/1‑90 (MQ=255) aatGAAAGCGTGTCCCTTGACACAAATAAAAATTTTCCAACCCATTTATCTCCTCAATTTTGGCTTCCACCCGAAAAACCCAGGTCACCg < 2:549651/90‑1 (MQ=255) tgtCCTTTGACACAAAAAAAAAATTTTCAACACACTTATCTACTCAATTTTGGCTTGCATCCGAAAAAACCAGGTCACCGCCATTAAAAt < 2:155850/90‑1 (MQ=255) cccTTGACACAAAAAAAAAATTTTCAACACCTTTATCTTCTAAATTTTGGCTCGCCCCCAAAAAAACCAGGTCACCGCCATTAAAATAGc < 2:285937/90‑1 (MQ=255) cccTTAACAAAAATAAAAATTTTCCAACCCATTTTTCTCCTCGATTTTGTCTTCCACCAGAAAAACCCAGGTCCCCGCCATTAAAATAGc < 2:351655/90‑1 (MQ=255) cTTGACACAAAAAAAAATTTTTCAACCCACTTTTCTACTAAATTTTGGCTTGCATCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGa < 2:71823/90‑1 (MQ=255) aaaaaatTTTTCAACACATTTTTCTACTCGATTTTGGCTTGCACCCCAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCAtt < 1:644104/87‑1 (MQ=255) aaaaaaaTTTCCAACACACTTATCTACTCAATTTTGGCTTCCACCCAAAAAAACCAGGCCACCGCCATTAAAATAGCGATATTCACCAtt < 2:588231/87‑1 (MQ=255) acTTATCTCCTCAATTTTGGCTTGCACCCGAAAAAACCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAatat < 1:390778/90‑1 (MQ=255) aTCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAttt > 1:132140/1‑90 (MQ=255) cGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCtt > 2:279401/1‑90 (MQ=255) cGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCtt > 1:427889/1‑90 (MQ=255) cGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCtt > 1:116800/1‑90 (MQ=255) tGGTTTCCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc < 1:201496/90‑1 (MQ=255) tGGCTTCCACCCAAAAAACCCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc < 2:140005/90‑1 (MQ=255) tGGCTTCCACCCAAAAAAACCAGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc < 1:510278/90‑1 (MQ=255) tGGCTTCCACCCAAAAAAACCAGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc < 2:427889/90‑1 (MQ=255) | CATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAC > NZ_CP009273/3658976‑3659138 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |