| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,005,978 | 0 | C | A | 42.5% | 84.3 / 15.7 | 48 | intergenic (‑192/‑284) | ygeV/ygeW | putative sigma‑54‑interacting transcriptional activator/putative carbamoyltransferase |
| Reads supporting (aligned to +/- strand): ref base C (8/19); new base A (20/0); total (29/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.76e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ATGATGATATTAACTCAGGCGAAATTGGCTTTGATAAAAACATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTC > NC_000913/3005848‑3006123 | atgatgATATTAACTCAGGCGAAATTGGCTTTGATAAAAACATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGGTTATCAAAATGAAGAGAGATAgc‑ccgccc > 2:119911/1‑151 (MQ=255) tgatgaTATTAACTCAGGCGAAATTGGCTTTGATAAAAACATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTTTTAAACTTAAGAGAGAAAGC‑CTGaccc > 1:679272/1‑151 (MQ=255) tgatgaTATTAACTCAGGCGAAATTGGCTTTGATAAAAACATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGGGGGTTTATTAAAATGAAGAAATATAGCTCTGcc > 1:476448/1‑150 (MQ=255) gatgatATTAACTCAGGCGAAATTGGCTTTGATAAAAACATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCt < 1:577669/151‑1 (MQ=255) aaCTCAGGCGAAATTGGCTTTGATAAAAACATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTa < 2:433109/150‑1 (MQ=255) cTCAGGCGAAATTGGCTTTGATAAAAACATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTTTTAAAATGAAGAGAGATAGC‑CTGCCCCCTTATCTTAttt > 2:382902/1‑151 (MQ=255) cTCAGGCGAAATTGGCTTTGATAAAAACATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGGTTATCAAAATTAAAAGAGATAAA‑CTGCCCCCTTTTCtttttt > 1:158414/1‑151 (MQ=255) aGGCGAAATTGGCTTTGATAAAAACATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGGGTTTATCAAAATGAAAAGAGATAGC‑CTGCCCCTTTAtttttttttt > 1:36882/1‑148 (MQ=255) ggCGAAATTGGCTTTGATAAAAACATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGa < 2:828932/151‑1 (MQ=255) gAAATTGGCTTTGATAAAAACATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATAc < 1:701649/151‑1 (MQ=255) ggCTTTGATAAAAACATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGGGTTTTTCAAAATGAAGAGAGATAGC‑CTGCCCCTTTATTTTATTTTTGATAAtt > 1:873104/1‑147 (MQ=255) ggCTTTGATAAAAACATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTagca < 2:464090/151‑1 (MQ=255) gATAAAAACATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGCGAGTTTATCAAACAGAAGAGAGATAGC‑CTGCCCCCTTATCTTTTTTTTTCTACTTAGCAGCcaa > 2:93741/1‑150 (MQ=255) aaaaaCATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCaaataaa < 2:879766/151‑1 (MQ=255) aaaaaCATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCaaataaa < 2:667708/151‑1 (MQ=255) aaaCATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGGTTATCAAAATGAAGAGAAAAAGC‑CTGCCCCCTTATCCTATTTCTGATAATTAGCAGCaaataaa > 2:147645/1‑149 (MQ=255) aaCATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAAATGAAGAGAGATAGC‑CTGCCCCTTTATCTTTTTTTTGATACTTTACAGCaaataaat > 2:619172/1‑149 (MQ=255) cATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGTGGCTTTCAAATCGTACAGAAATAAA‑CTGCCCCTTTATCTTATTTTTGTTACTTAGCAGTAaataaata > 1:228474/1‑148 (MQ=255) cATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAAcg < 2:714381/151‑1 (MQ=255) tttttATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTTTCAAACTTAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTTTGATACTTAACAGCAAATAAAAAACGCGATaa > 2:646159/1‑150 (MQ=255) ttttATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGTGTTTTTCAAAATGAAAAAAGATAAC‑CAGCCCCCTTTTTTTTTTTTTGATACTTAAGAGCAAaaaaataaaacgaataa > 2:506818/1‑141 (MQ=255) ttttATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATaaa < 1:628590/150‑1 (MQ=255) tATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTTTCAAACTTTAGAGAGATAAA‑CTGCCCCTTTATTTTTTTTTTTATACTTTACAGGAAAAAAATAACGCGATaaaaaa > 1:857126/1‑150 (MQ=255) tATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATTAAAATGAAGAGAGATAGC‑CTGCCCCCTTTTTTTATTTTTTATACTTAGCAGCAAATAAAAAACGCGATaaaaaa > 2:789266/1‑150 (MQ=255) tATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATaaaaaa < 1:375781/150‑1 (MQ=255) aTCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATaaaaaaa < 1:828429/150‑1 (MQ=255) aTCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGAGAGTTTTTCTCAATGAGAAGAGATAGC‑CCGCCCCCTTATCTTATTTCTTATACTTTACAGCAAATAAATAACACGATaaaaaaa > 2:240757/1‑150 (MQ=255) tCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAAATGAAAAGAGAAAGA‑CTGCCCCTTTTTTTTATTTTTGATACTTAGCAGAAAATAAATAAAGCGATaaaaaaa > 1:186370/1‑149 (MQ=255) aTTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGGGGGGTTATCAAAACGAAAAGAAAAAGC‑CTGCCCCTTTTTTTTATTTTTTTTAATTTAAAACAAAAAAATAAATCGATAAAAAAAgac > 1:356739/1‑148 (MQ=255) tttCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAAATTAAGAGAGATAAC‑CTGCCCCTTTTTTTTTTTTTTTATACTTAGCAGCAAATAAATAACGCGATAAAAAAAACCaa > 2:529106/1‑150 (MQ=255) tAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGATTTTTAAAATGAAAAGAGAAAAC‑CTGCCCCCTTATTTTTTTTTTGATACTTTACAGCAAATAAAAAAAACGATaaaaaaaaaccaaagtt > 2:180486/1‑141 (MQ=255) aaaTTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGCGGGTTTTTCAAAATTAAGAGTGGAAGC‑CCGCCCCCTTTTATTATTTCTTAAACTCACCAGGAAATAAAAAACCTGCTAAAAAAATGCAAACGgttt > 2:216389/1‑148 (MQ=255) aTGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTAtttt < 1:575905/151‑1 (MQ=255) aTGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGGGGGTTTATCAAAATGAAGAGAGAAAGA‑CCTCCCCTTTATTTTATTTTTTCTACTTAGCAGCAAAAAAAAAACGCGATAAAAAAAACCAAACGTTTTCGTAtttt > 1:631116/1‑151 (MQ=255) aaGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGGTTATCAAAATTAAAAGAGATAGC‑CTGTCCCTTTTTCTTATTTTTGTTAATTAAAAGCAAAAAAAAAAAGCGATAAAAAAAACCAAAAGTTTTCGTAtttt > 1:764312/1‑147 (MQ=255) agaTATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGGTTATTAAAATTAAGAAAGATAAC‑CCGCCCCTTTATTTTATTTTTGATACTTAGCAGCAAATAAATAACACGGTAAAAAAAACCAAAAGTTTTTGTATTTTacaaa > 2:521098/1‑149 (MQ=255) gaTATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTacaaacaa < 1:125789/151‑1 (MQ=255) atCACATTTCTATATCAATATGAGAATTACGGCGGTGAGGTTTTCAAAATGAAGAGAGATAGC‑CTGCCCCTTTTTCTTATTTTTGATACTTTGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAAcc > 2:160190/1‑150 (MQ=255) ttCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCt < 1:412698/150‑1 (MQ=255) tCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTg < 2:760231/150‑1 (MQ=255) tatCAATATGAGAATTACGGCGGGGAGTTTATTCAACTGATGAGAGAAATA‑CCGCTCCCTTTTTTTATTTTTGTTAATTAGCAGCAAAAAAATAACGAGGTTACAAAAGCCCCACATTTTTGTATTTTTCAAAAAAACAGAAGCttg > 1:193821/1‑145 (MQ=255) tCAATATGAGAATTACGGCGGTGAGGTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTTTTTTTTATACTTAGCAGCAAATAAAAAACGCGATAAAAAAAACCAAACGTTTTCGTATTTTACAAAAAACCAGAAGCTgg > 2:479507/1‑145 (MQ=255) cAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAAt < 2:606712/151‑1 (MQ=255) cAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAAt < 1:594930/151‑1 (MQ=255) aaTATGAGAATTACTGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAAtt < 1:644506/151‑1 (MQ=255) aCGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAAccc < 1:833321/151‑1 (MQ=255) ggTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCcacac > 2:874490/1‑151 (MQ=255) tCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTc > 2:20367/1‑151 (MQ=255) | ATGATGATATTAACTCAGGCGAAATTGGCTTTGATAAAAACATAAGATTTTTATCATTTTCTAATGAAATTATGGAAGAGATATCACATTTCTATATCAATATGAGAATTACGGCGGTGAGTTTATCAAACTGAAGAGAGATAGC‑CTGCCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTC > NC_000913/3005848‑3006123 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |