Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,659,062 | 0 | C | A | 45.0% | 28.8 / 11.0 | 20 | A37S (GCG→TCG) | gadX | acid resistance transcriptional activator GadX |
Reads supporting (aligned to +/- strand): ref base C (10/1); new base A (0/9); total (10/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.19e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAAT > NZ_CP009273/3658975‑3659143 | gCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaa < 2:193025/90‑1 (MQ=255) cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAc > 1:128218/1‑90 (MQ=255) gggAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAcc > 1:90397/1‑90 (MQ=255) aaGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCAtt > 2:448400/1‑90 (MQ=255) gtgtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaaa > 2:245229/1‑90 (MQ=255) cacaAATAAAAATTTTCCAACACCCTTATCTACTCGATTTTGGCTTGCATCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTc < 1:245229/90‑1 (MQ=255) cacaAATAAAAAATTTTCAACCCCCTTATCTACTAGATTTTGTCTTCCACCCAAAAAAACCAGGTCCCCGCCATTAAAATAGCGATATTc < 1:448400/90‑1 (MQ=255) caaaaaTAAAAAATTTTAAACACACTTATCTACTCAAATTTGGCTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTc < 2:128218/86‑1 (MQ=255) acGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCa > 2:67222/1‑90 (MQ=255) aTACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAccct > 1:309191/1‑88 (MQ=255) aCAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTa > 2:388508/1‑90 (MQ=255) ttaaaaCCCATTTATCTACTAAATTTTGGCTTGCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATgg < 1:28658/86‑1 (MQ=255) ccAACACACTTACCACCCCAATTTTGGTTTCCCCCCAAAAAAACCAGGCCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 1:41463/89‑1 (MQ=255) acacTTATCTCCTCAAATTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat < 2:90397/90‑1 (MQ=255) catttATCTACTCAATTTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAata < 2:237106/87‑1 (MQ=255) ttATCTCCTCAATTTTGGCTTGCATCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAatatat < 2:443301/90‑1 (MQ=255) cTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTc > 2:69863/1‑90 (MQ=255) cTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTc > 1:15619/1‑90 (MQ=255) tgttttGCACCCAAAAAAACCGGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc < 1:557049/86‑1 (MQ=255) tGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAAt > 2:499289/1‑90 (MQ=255) | GCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAAT > NZ_CP009273/3658975‑3659143 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |