Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,659,062 | 0 | C | A | 43.5% | 29.0 / 13.7 | 23 | A37S (GCG→TCG) | gadX | acid resistance transcriptional activator GadX |
Reads supporting (aligned to +/- strand): ref base C (10/3); new base A (0/10); total (10/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.64e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAC > NZ_CP009273/3658975‑3659148 | gCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaa < 1:386514/90‑1 (MQ=255) gCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaa > 1:457897/1‑90 (MQ=255) aTCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaa > 2:98679/1‑90 (MQ=255) tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa > 2:299116/1‑90 (MQ=255) taatGAAAGCGTGTCCTTTGACACAAAAAAAAAATTTTCAACCCCCTTATCTACTCGATTTTGGCTTGCATCCGAAAAACCCAGGTCAcc < 2:207556/90‑1 (MQ=255) aaaggAACGCGTGTCCCTTGACACGAATACAAAATTTTCACCACACTTATCTACCCGAATTTGGCTTGCACCCGCAAAACCCGGGTCAcc < 2:166444/86‑1 (MQ=255) tGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGcc > 1:112655/1‑90 (MQ=255) tGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGcc > 2:307331/1‑90 (MQ=255) gtgtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaaa > 1:471538/1‑90 (MQ=255) gtgtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaaa > 2:11623/1‑90 (MQ=255) cctttGACACAAATAAAAAATTTTCAACACCTTTATCTACTCGATTTTGGCTTGCATCCAAAAAAACCAGGTCACCGCCATTAAAATAGc < 2:112655/87‑1 (MQ=255) cacaAATAAAAAATTTTCAACACTTTTTTCTACAAAATTTTGGCTTCCACCCAAAAAAACCAGGCCACCGCCATTAAAATAGCGATATTc < 2:471538/90‑1 (MQ=255) acGAAAACAAAATTTTAAACACCTTTATCTCCTCAATTTTGGCTTCCATCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCa < 1:463913/90‑1 (MQ=255) aaaaTTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAAc > 1:405577/1‑90 (MQ=255) aaTTTCCAACACCCTTATCTACTCAATTTTGGCTTCCACCCAAAAAAACCAGGCCCCCGCCATTAAAATAGCGATATTCACCATTAACCa < 2:405577/90‑1 (MQ=255) ttccAACACATTTCTCTCCTCAATTTTGTCTCGCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATgg < 1:495798/87‑1 (MQ=255) tCAACACCCTTATCTACTCAAATTTGGCTTGCACCCGCAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 1:205234/90‑1 (MQ=255) tCAACACATTTATCTACTCAATTTTGGCTTGCACCCAAAAAACCCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 2:79994/90‑1 (MQ=255) ccAACACTTTTATCTCCTCAATTTTGGCTTCCACCCGAAAAACCCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 1:277590/89‑1 (MQ=255) cacTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAata > 2:459617/1‑90 (MQ=255) aCTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGt > 1:208539/1‑90 (MQ=255) tGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc < 1:98679/90‑1 (MQ=255) ccGAAAAACCCAGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc < 1:11623/90‑1 (MQ=255) | GCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAC > NZ_CP009273/3658975‑3659148 |
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: ‑ |