| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 353,586 | 0 | A | T | 56.2% | 0.8 / 41.5 | 32 | N272K (AAT→AAA) | cynR | transcriptional regulator CynR |
| Reads supporting (aligned to +/- strand): ref base A (13/1); new base T (0/18); total (13/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.03e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-05 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTT > NZ_CP009273/353497‑353662 | aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa > 1:1027291/1‑90 (MQ=255)aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa > 2:732832/1‑90 (MQ=255)aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa > 1:955489/1‑90 (MQ=255)aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa > 1:32372/1‑90 (MQ=255) cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt > 2:625582/1‑90 (MQ=255) cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt > 1:1209045/1‑90 (MQ=255) cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt > 1:429917/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) | AGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTT > NZ_CP009273/353497‑353662 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |