| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,184,646 | 0 | T | C | 56.5% | 13.0 / 27.7 | 23 | E219G (GAA→GGA) | thiM | hydoxyethylthiazole kinase |
| Reads supporting (aligned to +/- strand): ref base T (6/4); new base C (0/13); total (6/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.08e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-05 | |||||||||||
AGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACAAAACTGCCTGGCCCCTCGCTTCTGGCGACTGCGCGTTCTCCGGCTTGTTTCATCCAGTGACAGGCAGATGCGACATTTTCCAGCGTATCGCCTGGTAA‑CGCACAGCAGGCAGC‑GACAACCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAACGGATCACCACCGTGAATACCAATGATACGATGTCCATCGGTAACATAATCCATCTCGCCAGTGA > NC_000913/2184499‑2184793 | agCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACAAAACTGCCTGGCCCCTCGCTTCTGGCGACTGCGCGTTCTCCGGCTTGTTTCATCCAGTGACAGGCAGATGCGACATTTTcc < 2:322932/150‑1 (MQ=255) gCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACAAAACTGCCTGGCCCCTCGCTTCTGGCGACTGCGCGTTCTCCGGCTTGTTTCATCCAGTGACAGGCAGATGCGACATTTTCCAg > 2:621309/1‑151 (MQ=255) ccTGCACCTCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACAAAACTGCCTGGCCCCTCGCTTCTGGCGACTGCGCGTTCTCCGGCTTGTTTCATCCAGTGACAGGCAGATGCGACATTTTCCAGCGTATCGCCTGGTAA‑CGc < 2:81105/151‑1 (MQ=255) tGCACCTCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACAAAACTGCCTGGCCCCTCGCTTCTGGCGACTGCGCGTTCTCCGGCTTGTTTCATCCAGTGACAGGCAGATGCGACATTTTCCAGCGTATCGCCTGGTAA‑CGcac < 1:594527/151‑1 (MQ=255) aCCTCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACAAAACTGCCTGGCCCCTCGCTTCTGGCGACTGCGCGTTCTCCGGCTTGTTTCATCCAGTGACAGGCAGATGCGACATTTTCCAGCGTATCGCCTGGTAA‑CGCAcagc < 2:229235/151‑1 (MQ=255) ggcgaATCAAGGAAAGGGGAAAAAAACCTGCCGGCCCCCTCGCTTCGGGCGACTGCGCGTTCCCCGGTTTTTTTCACCCAGGGACAGGCAGATGCGACATTTCCCAGCGTACCGCCGGGTAC‑CGCACAGCAGGCAGC‑GACAACCGCCGAt < 2:320539/149‑1 (MQ=255) aGGAAATGTGGAACAAAACTGCCTGGCACCTCGCTTCTGTCGACTGCGCGTTCTCCGGCTTGTTTCATCCAGTGACAGGCAGATGCGACATTTTCCAGCGTATCGCCTGGTAA‑CGCACAGCAGGCAGC‑AACAACCGCCGGTAATGCACAg > 2:552207/1‑150 (MQ=255) gAAATGTGGAACAAAACTGCCTGGCCCCTCGCTTCTGGCGACTGCGCGTTCTCCGGCTTGTTTCATCCAGTGACAGGCAGATGCGACATTTTCCAGCGTATCGCCTGGTAA‑CGCACAGCAGGCAGC‑GACAACCGCCGATAATGCACAGCCa > 1:587414/1‑151 (MQ=255) aaCAAAATTCCCTCCCCCCTTGTTTCTGGCCCCTGGGTGTTCCCCGTCTTGTTCCCCCCAGTGACAGGGAGAGGCAACATTTCCCGTGGTACCGCCTGTAAA‑CGCACAGTAGTATCC‑GACATCCGCCGATAATGCACAGCCAGTTCCTAcc < 1:434635/151‑1 (MQ=255) ccccccGTTTCGGGGGCCGGCGCTTCCCCCGGCTTGTTTCACCCAGGGACGGGAAGAGGCAAATTTTCCCGGCTTCCCGCCGGGAAC‑CCCACAGCAGGCAGC‑AACAACCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAAc < 2:574191/151‑1 (MQ=255) ccccTCGCTTCGGGCGCCGGCGCTTTCCCGGGTTTTTTTCACCCAGGGACAGGCAGAGGCGCCATTTCCCGGCGTACCGCCGGGTAC‑CGCACGGCAGGCGGC‑GACAACCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAAc < 2:172840/151‑1 (MQ=255) cTCGTTTCGGGCGCCGGCGCTTCCCCCGGTTTGTTCCACCCAGGGACGGGCAGAGGCGACTTTTCCCAGCGTACCCCCGGGTAC‑CGCACAGCGGGCAGC‑AAAAACCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAACGGa < 1:433341/151‑1 (MQ=255) cTTCTGGCGACTGCGCGTTCTCCGGCTTGTTTCATCCAGTGACAGGCAGATGCGACATTTTCCAGCGTATCGCCTGGTAA‑CGCACAGCAGGCAGC‑GACAACCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAACGGATca > 2:485870/1‑150 (MQ=255) cgggcggcggcGCGTTCCCTGGCTTTTTCCACCGGGAAACAGGAAGAGCCAACATTTCCCAGCTTA‑CCCCCGGTAACCCCCCAGCCGGCAGC‑AACAACCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAACGGATcaccac < 2:377661/149‑1 (MQ=255) cgggCGCCGGCGCTTCCCCCGGTTTTTTTCACCCAGGGACGGGCAGATGCGACTTTTCCCAGCTTATCGCCTGGAAA‑CGCACAGCAGGCAGC‑AAAAACCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAACGGATcaccac < 1:205712/149‑1 (MQ=255) ccggcgcGTTCCCCGTTTTTTTTCACCCAGGGACAGGCAGATGCGACTTTTCCCAGCGTACCGCCTGGAAC‑CGCACAGCAGGCGGC‑AACACCCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAACGGATCACCACCGTGa < 1:57145/147‑1 (MQ=255) ctttCTCGGGTTTTTTTCACCCAGTGACAGGCAGAGGCGACATTTCCCAGCTTATCGCTTGGAAC‑CGCACAGAAGCCAGC‑GACAACCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAACGGATCACCACCGTGAATACCaa < 1:621321/149‑1 (MQ=255) gTTTCATCCAGTGACAGGCAGATGCGACATTTTCCAGCGTATCGCCTGGTAA‑CGCACAGCAGGCAGC‑GACAACCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAACGGATCACCACCGTGAATACCAATGATACGATGTc > 1:530329/1‑150 (MQ=255) aTCCAGAGACAGGCAGATGCGACTTTTCCCGGCGTACCGCCGGTTAA‑CCCACAG‑AGGGAGCAGACAACCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAACGGATCACCACCGTGAATACCAATGATACGATGTCCATCgg < 1:116297/151‑1 (MQ=255) tGACGGGCAGATGCGACTTTTCCCAGCGTACCGCCTGTTAA‑CGCACAGCAGCCAGC‑AAAAACCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAACGGATCACCACCGTGAATACCAATGATACGATGTCCATCGGTAACAt < 1:418446/151‑1 (MQ=255) gACAGGCAGATGCGACTTTTCCCAGCGTACCGCCTGGTAA‑CCCACAGCAGGCAGC‑GACAACCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAACGGATCACCACCGTGAATACCAATGATACGATGTCCATCGGTAACAt < 1:537443/150‑1 (MQ=255) gCAGATGCGACATTTCCCGGCATCCCGCCGGGAAC‑CCCACAGCGGGCGGC‑AACAACCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAACGGATCACCACCGTGAATACCAATGATACGATGTCCATCGGTAACATAATcc < 2:288720/150‑1 (MQ=255) aCATTTTCCAGCGTATCGCCTGGTAA‑CGCACAGCAGGCAGC‑GACAACCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAACGGATCACCACCGTGAATACCAATGATACGATGTCCATCGGTAACATAATCCATCTCGCCAg > 1:378836/1‑151 (MQ=255) ttttCCAGCGTATCGCCTGGTAA‑CGCACAGCAGGCAGC‑GACAACCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAACGGATCACCACCGTGAATACCAATGATACGATGTCCATCGGTAACATAATCCATCTCGCCAGTGa > 1:377972/1‑151 (MQ=255) | AGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACAAAACTGCCTGGCCCCTCGCTTCTGGCGACTGCGCGTTCTCCGGCTTGTTTCATCCAGTGACAGGCAGATGCGACATTTTCCAGCGTATCGCCTGGTAA‑CGCACAGCAGGCAGC‑GACAACCGCCGATAATGCACAGCCAGTTCCTACCACTTTGGTCATTAACGGATCACCACCGTGAATACCAATGATACGATGTCCATCGGTAACATAATCCATCTCGCCAGTGA > NC_000913/2184499‑2184793 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |