Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,412,781 | 0 | C | G | 41.5% | 63.9 / 15.8 | 41 | S132C (TCT→TGT) | yjfM | DUF1190 domain‑containing protein YjfM |
Reads supporting (aligned to +/- strand): ref base C (1/23); new base G (17/0); total (18/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.19e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCT > NC_000913/4412651‑4412914 | cgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGttttttt > 1:1941780/1‑140 (MQ=255) cgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGAtttttt > 1:314155/1‑140 (MQ=255) cgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGttttttt > 1:314091/1‑140 (MQ=255) cgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGTTTTTTg > 1:774559/1‑141 (MQ=255) cgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGGGTGGAtttttt > 1:275469/1‑140 (MQ=255) aacGCCAAAACGGCATTTTATGCCGATGTTCCGTAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGtt < 2:832338/141‑1 (MQ=255) cGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGGTTTTTTTTTATCCggggttt > 1:539647/1‑137 (MQ=255) ttATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGTTTTTTTTTTATGCGGGGTATTGGTTAAGATAGCGa > 1:1634238/1‑141 (MQ=255) ccGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGATTTTTGTTTAGCCCGGTGATTCGTAAAGAAAGCGAttatc > 2:116996/1‑137 (MQ=255) tGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCgtt < 1:1171994/141‑1 (MQ=255) tGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCgtt < 1:1808322/141‑1 (MQ=255) tGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCgtt < 1:1808312/141‑1 (MQ=255) tGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCgtt < 1:173196/141‑1 (MQ=255) tGTTCCGAAGAAGACGACTCAGCAGCACTGTCAGTCGAAGTACGAAAAGTGCTACTACGACAAGGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCgtt < 2:882655/141‑1 (MQ=255) ttCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCgtttg < 2:1210143/141‑1 (MQ=255) ttCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCgtttg < 2:1210234/141‑1 (MQ=255) gcagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGTTTTTTTTTAATACGGGGTATTAGTAAAGATAGCGATGATCCGGTTTTTTTTTAAAGcggcgggtg > 1:1829402/1‑138 (MQ=255) gcagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGGTTTTTTTTTATCCGGGTTTTTGGTTAATATAGCGATGTGCCGTTTGTTTTTAAAAGcggcgggta > 1:1829044/1‑138 (MQ=255) gcagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGGTTTTTTTTTAGCCGGGGTTTTAGGAAAGATAGCGGTTAGCCGTTTGTTTTTAACACcggcgggtg > 1:1829329/1‑138 (MQ=255) gcagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGATTTTTTTTTAGCCGGGGTTTTGGTAAAGATCGCGAGTTTGCGGGTGTTTATTAAAGcggcgggtg > 1:1683935/1‑138 (MQ=255) gTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGTTTTTTGTTTATCCCGGGTATTAGTAAATATCGCGATGAGCCGGTTGGTTTTAAAAGCGGCGGGTGCGCTTTTg > 1:805838/1‑141 (MQ=255) gTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGATTTTTGTTAAGACGGGGTATTAGTAAAGATAGCGGTGTGCCGGTTTTTTATAAAAGCGGCGGGTCCACttttt > 1:349828/1‑140 (MQ=255) gTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTc < 1:887241/141‑1 (MQ=255) gTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTc < 1:887275/141‑1 (MQ=255) gTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTGTGATTAGAACAGCGGCGGTTCCTCTTTTGCtgc < 2:494078/141‑3 (MQ=255) tCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcg < 2:1258321/141‑1 (MQ=255) tCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcg < 2:165958/141‑1 (MQ=255) tCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcg < 2:1259610/141‑1 (MQ=255) aTTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGATTTTTTTTAATCCGGGTTATTAGTAAATATAGCGTTGAGCCGTTTGTTTATAACAACGGCGGTTACTCTTTTTGTTAGCGCCCaggcaggc > 2:1842691/1‑141 (MQ=255) gCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgca < 2:1958566/141‑1 (MQ=255) gCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgca < 2:788503/141‑1 (MQ=255) gCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgca < 2:786356/141‑1 (MQ=255) gCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTACTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgca < 2:788476/141‑1 (MQ=255) tactacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgcagc < 1:1078625/141‑1 (MQ=255) acGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACtt < 2:1634238/141‑1 (MQ=255) cAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTgg < 2:414754/141‑1 (MQ=255) gTTGAACAGAGTTGGATCCCGGTCGTTTGTGGATTTTTGTTAATCCGGGTTATTCGTAAAGATAGCGATGAGCCGTTTGTTTATAAAAGCGGCGGTTACACTTTTGCGTAGCGCCCACTCTGGCGCACCACTTATGGTGAt > 2:1602057/1‑141 (MQ=255) ccGGTCGTGTGTGGTTTTTTGTTTATCCGGGTTATTCGTAAATATAGCGATGAGCCGTTTGTTTATAACAGCGGCGGGTGCTCTTTTTCTTAGCGGCCAGTCTTGGGCAGCACAGCTGGGGTTTAATACAGGCGCTCCCGt > 2:492207/1‑141 (MQ=255) ccGGTCGTGTGTGGATTTTTTTTAAACCGGGTTATTCGTAAAGATCGCGATGTGCCGTTTGGTTATTAAAACGGCCGTTACTATTTTGCATACACCCCAATATTGCCCACCACATCTGGGGATTACACCAGGCCCACCGGt > 2:1869877/1‑141 (MQ=255) gTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGTTGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCt < 2:987806/141‑1 (MQ=255) gTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCt < 2:1794194/141‑1 (MQ=255) | CGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCT > NC_000913/4412651‑4412914 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |