Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,429,972 | 0 | G | A | 57.9% | 4.5 / 18.9 | 19 | intergenic (‑96/+45) | uspF/ompN | universal stress protein UspF/porin OmpN |
Reads supporting (aligned to +/- strand): ref base G (8/0); new base A (0/11); total (8/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.32e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAG > NZ_CP009273/1429883‑1430060 | tcctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAaggagg < 1:131545/90‑1 (MQ=255) tcctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAaggagg > 1:486523/1‑90 (MQ=255) tATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCa > 1:99849/1‑90 (MQ=255) ccGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGaa > 1:425169/1‑90 (MQ=255) cGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGaaa > 2:708609/1‑90 (MQ=255) aaGTGTAACCCAAGATGCCTTTTTTTGTGTAAGGAGGGAGTTCACATTTTTAATTTAGAAAAAAGGAGAAAAAAAAAGAAAAAAAAGccc < 2:744650/90‑1 (MQ=255) tgtAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGcc > 2:546161/1‑90 (MQ=255) tAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCAt > 1:18903/1‑90 (MQ=255) aaCGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATa > 1:225027/1‑90 (MQ=255) ccAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAg > 2:808889/1‑90 (MQ=255) caaGATGACTTTTTTTGTGAAAGGAGGGATTTCATTTTTTTAATTAAAAAAAAGGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAg < 1:322787/88‑1 (MQ=255) aGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCg > 1:259029/1‑90 (MQ=255) tttttttGTAAAAGGAGGGATTTAAATTTTTAAATTAAAAAAAAAGGAGAAAAAAAAAAAAAAAAAACCCCGCCATAACAGCGGGCagga < 1:546161/88‑1 (MQ=255) ttttttGTGTAATGACGGAGTTCACTTTTTTAATTTAGATAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCaggag < 2:239896/89‑1 (MQ=255) ttttttGTGAAAGGAGGAATTTCAATTTTTAAATTAAGAAAAAAGGAGAAAAAAAAAAAAAAAAAACCCCGCCATAACAGCGGGCaggag < 1:489832/89‑1 (MQ=255) aaaggACGGATTTCATTTTTTTAATTAAGAAAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGaa < 1:258759/86‑1 (MQ=255) aCGGAGTCCACATTTTTAATTTAAAAAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGAt < 2:138938/90‑1 (MQ=255) cGGATTTCAATTTTTTAATTAAAAAAAAAGGAGAAAAAAAAAAAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATa < 2:225027/90‑1 (MQ=255) tttAATTAAGAAAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaaa < 2:594955/90‑1 (MQ=255) ttaaaTTTAAACAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaaa < 2:333048/87‑1 (MQ=255) aaGGAGAAAAAAAAAAAAAAAAACCCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatcatc < 1:808889/90‑1 (MQ=255) ggaAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTag > 2:892926/1‑90 (MQ=255) | TCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAG > NZ_CP009273/1429883‑1430060 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |