Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,659,058 | 0 | T | C | 58.8% | 9.0 / 17.6 | 17 | D38G (GAT→GGT) | gadX | acid resistance transcriptional activator GadX |
Reads supporting (aligned to +/- strand): ref base T (6/1); new base C (0/10); total (6/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.66e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAC > NZ_CP009273/3658973‑3659138 | gAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGc > 1:77798/1‑90 (MQ=255) gCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaa < 2:214999/90‑1 (MQ=255) tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa > 1:304129/1‑90 (MQ=255) tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa > 1:407726/1‑90 (MQ=255) taataaAAGCGTGTCCCTTGACACAAATAAAAAATTTTCAACACCCTTATCTACTCAAATTTGGCTTGCACCCGAAAAAACCAGGTCAcc < 2:215443/90‑1 (MQ=255) aaGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCAtt > 2:487900/1‑90 (MQ=255) cctttGACACGAAAAAAAATTTTTCAACACATTTATCTCCTCGAATTTGGCTTCCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGc < 1:26978/87‑1 (MQ=255) acacAAAAAAAAATTTTTCAACACATTTATCTCCTCGATTTTGGCTTGCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATAtt < 1:381268/90‑1 (MQ=255) cacaAATACAAAATTTTCAACCCATTTATTTTCCCGATTTTGGTTTGCACCCGAAAAACCCGGGTCACCGCCATTAAAATAGCGATATTc < 2:491996/90‑1 (MQ=255) cacaAATAAAAAATTTTCAACACACTTATCTCCTCGAATTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTc < 1:176270/90‑1 (MQ=255) gAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAcc > 2:293367/1‑90 (MQ=255) ttttCACCACACTTACCTACTCAATTTTGGTTTGCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATg < 2:254003/90‑1 (MQ=255) tCAACACATTTATCTCCCCAATTTTGGCTTCCACCCGAAAAAACCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 1:298021/90‑1 (MQ=255) caaaCACACTTACCTCCAAAAATTTGGCTTCCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg < 2:11115/88‑1 (MQ=255) cTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTc > 2:398815/1‑90 (MQ=255) tGGCTTCCACCCGAAAAACCCGGGCCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc < 1:293367/90‑1 (MQ=255) tGGCTTCCACCCAAAAAACCCGGGCCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc < 2:23909/90‑1 (MQ=255) | GAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAC > NZ_CP009273/3658973‑3659138 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |