Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,825,061 | 0 | G | A | 25.6% | 71.3 / 11.2 | 39 | intergenic (‑124/+79) | ves/spy | HutD family protein Ves/ATP‑independent periplasmic chaperone |
Reads supporting (aligned to +/- strand): ref base G (26/3); new base A (0/10); total (26/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.50e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTT > NC_000913/1824922‑1825196 | tatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaag > 1:1460960/1‑140 (MQ=255) aGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAcc > 1:1496855/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa > 1:1456943/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa > 1:880655/1‑140 (MQ=255) cATCCTCGCTCCCGGCCTTTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTGGGGTGA‑TTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAAAAGAAAAAAGAAAACCGCCGATc < 1:1359582/140‑1 (MQ=255) tccTCGCTCCCGGCCCGTTCTTTTCA‑TTTTTATCGTTCAGTTACCCGTGGATGGTTATCGCCACCTGGGGTGA‑TTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAACCCGCCGATCCt < 2:42500/138‑1 (MQ=255) aCTCGCTCCCGGCCCGTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTg > 2:924154/1‑140 (MQ=255) aCTCGCTCCCGGCCCGTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTg > 1:275082/1‑140 (MQ=255) cTCGCTCCCGGCCCGTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGt > 2:1538259/1‑140 (MQ=255) tCGCTCCCGGCCCGTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTc > 1:893618/1‑140 (MQ=255) cTCCCGGCCCGTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTGGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc > 2:1141734/1‑140 (MQ=255) cccGGCCCGTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCg > 1:888145/1‑140 (MQ=255) cccGGCCCGTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGGATGG‑TATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGc > 1:658671/1‑140 (MQ=255) cGGCCCGTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCa > 1:1115176/1‑140 (MQ=255) ggccttttCTTTTTA‑TTGTTATTGTACAGATAATGGTGTATGTTTATCGTC‑ATTGGGGTTATTTTACGTTATTGTTAAGCAGGTTTTTTTAACAAGTGTTAGAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCAt < 2:1115176/134‑1 (MQ=255) ggCCCGTTTTTTTTATTTTTTATCTTTAAGGTACCGGTGGATGTTTATCGGCACCTGGGGTGT‑TTTACGTTAGTGGTGATCAGTTTTTTTTCAAACCTTTTAGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCa < 2:315523/140‑1 (MQ=255) ggCCCGTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGGTCTTGTCCActacct > 1:864895/1‑135 (MQ=255) cccTTTTTTTTCA‑TTGTTATCGTTCAGGTAACGGTGGATGTTTATCGGAAACTGGGGTGA‑TTTACGTTAGTGGTGATCAGTTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGGTCTTGTCCActaccttg < 2:864895/140‑8 (MQ=255) cccTTTCTTTTCA‑TTGTTATCTTTCGGGTACCCGTGGATGTTTATCGGCAACTGGGGTGA‑TTTACGTTAGTGGTGATCGGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTa < 2:1333003/140‑1 (MQ=255) cgttttttttCT‑TTTTTATCGTTCAGGTAACGGGGGATGTTTATCGGCAACTGGGGTTA‑TTTACTTTAGTGGTGATCAGTTTTTTTTTACAATTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTAc < 2:977065/134‑1 (MQ=255) tcttttttATCGTTCAGGAAACGGTGGTTGGTTATCGGCACCTGGGGTGT‑TTTACTTTAGGGGTGATCGGTTTTTTTTTACATTTGTTGGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa < 2:1460960/134‑1 (MQ=255) tatttGTTTTCGTCCAGGTAACGGGGGATGTTTATCGGAATCGCGGGTGT‑TTTAGTTTGGTGGTGATCGGGTTTTTTTAACATTTTTTGAAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa < 2:880655/137‑1 (MQ=255) tATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa > 1:309363/1‑140 (MQ=255) ttCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcc > 1:1273716/1‑140 (MQ=255) ggTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCATCGCATTACTGCAAGGTAGTGGACAAGACcggcg > 2:1528205/1‑140 (MQ=255) gggTTACCGGCAACTGGGGTGT‑TTTACGTTAGGGGTGATCAGGTTTTTTTTACATCTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttt < 2:717706/139‑1 (MQ=255) aTCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGc > 1:1250327/1‑140 (MQ=255) cTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAgct > 1:717706/1‑138 (MQ=255) tttACGTTGGTGGTGATCGGTTTTTTTTAACAATTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagca < 2:275082/140‑1 (MQ=255) tttgtggtgCTCGGTTTTTTTTCACAGTTCTTTGAGAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGTATTTCTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa < 2:1456943/137‑1 (MQ=255) ttAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa < 2:1250327/140‑1 (MQ=255) tggtggtcACCAGTTTTTTTTAACAATTGTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg < 2:1496855/138‑1 (MQ=255) tAGGGGTGATCAGTTTTTTTTAAAAATTGTTGAAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg < 2:309363/140‑1 (MQ=255) gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGc > 2:1112609/1‑140 (MQ=255) gtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTTGCTGAAAGATTATTTAGCAGTTGCCGGCAtt > 2:88138/1‑140 (MQ=255) gATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtttt > 1:1334472/1‑140 (MQ=255) aTCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa > 1:117666/1‑140 (MQ=255) aTCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCATAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTTCAGGCAttttt > 2:815011/1‑139 (MQ=255) tttttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg > 1:1248298/1‑140 (MQ=255) ttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCg > 2:588567/1‑140 (MQ=255) ttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCg > 1:1551825/1‑140 (MQ=255) aaaatcTGTTAGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTg < 2:1551825/135‑1 (MQ=255) tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt > 2:158454/1‑140 (MQ=255) gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt > 1:1060416/1‑128 (MQ=255) gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt < 2:1060416/128‑1 (MQ=255) gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAAATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGGCGTTCTGTCAGACGCt > 2:1540594/1‑140 (MQ=255) aaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCtt > 1:628199/1‑140 (MQ=255) | TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGA‑TTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTT > 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 < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |