Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,659,062 | 0 | C | A | 62.5% | 6.1 / 22.2 | 16 | A37S (GCG→TCG) | gadX | acid resistance transcriptional activator GadX |
Reads supporting (aligned to +/- strand): ref base C (5/1); new base A (0/10); total (5/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.37e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-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: ‑ |