Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 353,591 | 0 | T | C | 55.2% | 16.8 / 33.2 | 29 | K271E (AAA→GAA) | cynR | transcriptional regulator CynR |
Reads supporting (aligned to +/- strand): ref base T (10/3); new base C (0/16); total (10/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.43e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGC > NZ_CP009273/353502‑353675 | cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt > 1:429917/1‑90 (MQ=255) cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt > 1:1209045/1‑90 (MQ=255) cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt > 2:625582/1‑90 (MQ=255) gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGTTGTCTGCCAGCTATTTTTccgcc > 2:329839/1‑90 (MQ=255) gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc < 1:812271/90‑1 (MQ=255) gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc > 1:1138327/1‑90 (MQ=255) gCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgc > 1:637342/1‑90 (MQ=255) cacccGCTCTTTTCTCCACCCCCATGAAAAAAAAGGCCTCCGGGGGGGTTGTCGGCAAGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTtc < 2:951807/87‑1 (MQ=255) cacccGCGCTTTTACCCACCCCCTTGTAAAAAAAGGCTTTCGGGGGGCTTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAA‑CCGCCGTtc < 1:118683/87‑1 (MQ=255) cacccGCGCTTTTACCCACCCCCATGTCAAAAAAGCCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTtc < 1:102477/87‑1 (MQ=19) cacccGCGCTTTTACCCAACCCCATGTGAAAAAATCCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAA‑CCGCCGTtc < 1:151883/87‑1 (MQ=21) cacccGCCCTTTTACCCACCCCCTTGTAAAAAAATCCCTCCGGGGGGGTTTTCTCCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTtc < 2:280635/87‑1 (MQ=255) ccGCGCTTTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTtctct < 1:90356/90‑1 (MQ=21) gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTcc > 1:794764/1‑90 (MQ=255) ccctaTTTACCCAACCCCATGTAAAAAAATCCTTTCGCGGCGGTTGTCGGCCATTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTcc < 1:697181/86‑1 (MQ=255) agaaTTTACCCACCCCCATGTAAAAAAATCCCTCCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCa < 1:1102463/87‑1 (MQ=21) accccATGTGCAAAAATGCCTTCGGGGGGGCTTTCTGCCAGTTTTTTTTCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg < 1:41075/87‑1 (MQ=37) aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggt > 2:1102463/1‑89 (MQ=255) aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg > 2:765090/1‑90 (MQ=255) aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg > 2:933182/1‑90 (MQ=255) aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCCGTGGTGGCGGGGccgga > 1:676079/1‑86 (MQ=37) aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTA‑TTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGagag > 2:379426/1‑90 (MQ=37) aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTA‑TTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGagag > 2:480905/1‑90 (MQ=37) aaaatTCCCCTCGCGGGGGTTGTCGCCCATTTTTTTTCCCCCCAAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa < 2:1131791/90‑1 (MQ=255) aaaaaTCCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa < 1:933182/90‑1 (MQ=38) aaaaaCCCCTCCGCGGGGGTTCTCGGCCGGTTTTTTTCCCCCCCCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa < 1:648260/90‑1 (MQ=255) tCGGGGCGGCTGTCTGCCAGCTTTTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc < 1:299548/90‑1 (MQ=255) tCGCGGGGGCTGTCTGCC‑GTTTTTTTTCCCCCGCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc < 1:1130163/90‑1 (MQ=255) ccggggcggcTTTCTCCCGGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc < 2:155858/86‑1 (MQ=255) gggcggtTGTCGGCCAGTTTTTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt < 1:226742/89‑1 (MQ=255) gcggcTGTCTCCCAGTTTTTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa < 1:1101121/90‑1 (MQ=255) ttGTCGGCCACTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgtt < 1:527248/89‑1 (MQ=255) tttttcCCGCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc < 1:379426/89‑1 (MQ=255) tttttcCCGCCGAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc < 1:485822/89‑1 (MQ=255) tttttcCCCCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc < 1:183814/89‑1 (MQ=255) tttttCCCCCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc < 2:472067/90‑1 (MQ=255) | CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGC > NZ_CP009273/353502‑353675 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |