Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,429,984 | 0 | G | A | 61.1% | ‑6.4 / 18.4 | 18 | intergenic (‑108/+33) | uspF/ompN | universal stress protein UspF/porin OmpN |
Reads supporting (aligned to +/- strand): ref base G (7/0); new base A (0/11); total (7/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.14e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCC > NZ_CP009273/1429904‑1430067 | aaGTGTACCGCAAAATGACTTTTTTTGTGAAAGGACGGAGTTCACTTTTTTAATTTAGATAAAAGGAGGAAAAAAAAGAAAAAAAAGccc < 1:295343/90‑1 (MQ=255) aaGTGAAACGCCAGATGACCTTTTTTGTGAAATGAGGGAGTTCACATTTTTAATTTAGATAAAAGGAGAAAAAAAAAGCAAAAAAAGccc < 2:41681/90‑1 (MQ=255) ccAGATGACTTTTTTTGTGAAAGGAGGGATTTCACTTTTTTAATTTAAAAAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAg < 1:233199/90‑1 (MQ=255) ttttttGTGTAAGGAGGGATTCCATTTTTTTAATTAAAAAAAAAGGAGAAAAAAAAAGAAAAAAACCCCCGCCATAACAGCGGGCaggag < 1:100480/89‑1 (MQ=255) ttttttGTATAATGACGAATTCCATTTTTTAAATTAAAAAAAAAGGAGAAAAAAAAAAAAAAAAAACCCCGCCATAACAGCGGGCaggag < 1:776935/89‑1 (MQ=255) ttttttGTAAAAGGAGGAATTTCAATTTTTAAATTAAAAAAAAAGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCaggag < 1:358907/89‑1 (MQ=255) ttttttGGGTAATGGCGGTGTCCATATTTTTAATAAAAAAAAAAGGAGAAAAAAAAAAAAAAAAAAGCCCGCCATAACAGCGGGCaggag < 1:594285/89‑1 (MQ=255) ttGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAtt > 2:610437/1‑90 (MQ=255) tAATGCGGGAGTTCATTTTTTTATTTAAGATAAAAGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGaa < 2:850770/90‑1 (MQ=255) tAATGACGGATTTCACTTTTTTAATTTAGATAAAAGGAGGAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGaa < 2:332424/90‑1 (MQ=255) tAATGACGGAGTTCACATTTTTAATTTAAAAAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGaa < 2:27063/90‑1 (MQ=255) tAAGGACGGAGTTCCCTTTTTAAATTTAAAAAAAAGGAGGAAAAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGaa < 2:51612/90‑1 (MQ=255) ggAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATaa > 1:300015/1‑90 (MQ=255) ggAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATaa > 1:332424/1‑90 (MQ=255) tCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAg > 2:772396/1‑90 (MQ=255) ttttAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaa > 2:9948/1‑90 (MQ=255) ggaAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTag > 1:187086/1‑90 (MQ=255) tAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGcc > 1:600349/1‑90 (MQ=255) | AAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCC > NZ_CP009273/1429904‑1430067 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |