| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,825,061 | 0 | G | A | 37.5% | 38.0 / 15.5 | 32 | intergenic (‑124/+79) | ves/spy | HutD family protein Ves/ATP‑independent periplasmic chaperone |
| Reads supporting (aligned to +/- strand): ref base G (17/3); new base A (0/12); total (17/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.02e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAAC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTT > NC_000913/1824922‑1825196 | tatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAAC‑TGTTTGAGAAGagaaaag > 2:581276/1‑140 (MQ=255) cccccGGCCCTTTTTTTTCTTTTTTATCGTTCAGGTACCGGTGGATGTTTATCGGCAACTGGGGTGTTTTACGTT‑AGTGGTGATCGGTTTTTTTTTACACT‑TTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCAc < 2:1133208/138‑1 (MQ=255) cTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAAC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc > 1:253919/1‑140 (MQ=255) tCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAAC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAcc > 1:1133208/1‑140 (MQ=255) ggccttttcttttcttTTTTCTCTTTCAGGTACCGGTGGATGTTTATCGGAACCTGGGGTGATTTACGTT‑GGGGTTAATCAGTTTTTTTTTACAAT‑TTTTGGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt < 2:1382117/134‑1 (MQ=255) gCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAAC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt > 1:720169/1‑140 (MQ=255) ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAAC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc > 2:670937/1‑140 (MQ=255) aaGGTAGCGGTGGATGTTTACCGGAAGCTGGGGTTATTT‑CGTTAAGGGGTGATC‑GGTTTTTTTTAAAACTTGTTAGAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg < 2:250162/139‑1 (MQ=255) aaGGTACCGGGGGATGTTTATCGGCACCTGGGGTGATTTACTTT‑AGGGGTGACCGGTTTTTTTTTACAAC‑TTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg < 2:9344/139‑1 (MQ=255) cGTGGATGGTTATCGGCAACTCGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAAC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTa > 1:1138035/1‑140 (MQ=255) aTGGTTATCGGCAACTCGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAAC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGgtt > 1:195236/1‑140 (MQ=255) cGGCAACTGGGGTGATTTACGTT‑AGTGGTGATCAGGTTTTTTTTACAAT‑TGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg < 1:581276/140‑1 (MQ=255) aaCTCGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAAC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAg > 2:738104/1‑140 (MQ=255) gATTTACGTT‑AGTGGTGATCAGGTTTTTTTTACAAC‑TGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag < 1:806980/140‑1 (MQ=255) gATTTACGTT‑AGTGGGGATCAGGTTTTTTTTACAAC‑TGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag < 2:1138035/140‑1 (MQ=255) ttACGTT‑AGTGGTGATCAGGCTTTTTTTACAAC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGGAAGATTATTcagcag > 2:746273/1‑140 (MQ=255) tt‑AGTGTTGATCAGGTTTTTTTTACAAC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa < 1:565716/140‑1 (MQ=255) gtggtgATCAGGCTTTTTTTACAAC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGAAGTTGCAGGc > 1:147431/1‑140 (MQ=255) tggtgATCAGGCTTTTTTTACAAC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCa > 1:211365/1‑140 (MQ=255) tggtgATCAGGCTTTTTTTACAAC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCa > 2:1233970/1‑140 (MQ=255) aGGCTTTTTTTACAAC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCt > 1:165696/1‑140 (MQ=255) tttttttACAAC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg > 1:1249771/1‑140 (MQ=255) tttACATT‑TTTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGc < 2:1249771/140‑1 (MQ=255) aC‑TTTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg < 2:1140227/140‑1 (MQ=255) aC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg < 1:738104/140‑1 (MQ=255) cttttttGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg < 2:165696/136‑1 (MQ=255) cttttttGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg < 2:196690/136‑1 (MQ=255) tGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt < 2:651325/140‑1 (MQ=255) tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt > 2:1030564/1‑140 (MQ=255) tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt > 1:1140227/1‑140 (MQ=255) tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt > 1:651325/1‑140 (MQ=255) tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCACTGGACGTTCt > 2:726760/1‑140 (MQ=255) aaaagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCa < 2:720169/138‑1 (MQ=255) aaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCtt > 1:362088/1‑140 (MQ=255) | TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAAC‑TGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTT > NC_000913/1824922‑1825196 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |