| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,747,168 | 1 | . | C | 20.0% | 45.2 / 15.6 | 20 | coding (85/1497 nt) | lyxK | L‑xylulose kinase |
| Reads supporting (aligned to +/- strand): ref base . (8/8); new base C (2/2); total (10/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGTGCCGCAGGGAGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTC > NC_000913/3747032‑3747305 | tgtgCCGCAGGGAGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGCTGCA‑CGCCTGc < 3:163856/149‑1 (MQ=255) gctgAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTgcgc < 8:159824/138‑1 (MQ=255) gctgAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTgcgc > 7:159824/1‑138 (MQ=255) ctgAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGCTGCA‑CGCCTGc > 7:59906/1‑127 (MQ=255) atgAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGCTGCA‑CGCCTGc < 8:59906/126‑1 (MQ=255) aaGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCa < 4:126607/149‑1 (MQ=255) cGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGGGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑GGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCt < 7:124326/149‑1 (MQ=255) tGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTgg < 3:33846/149‑1 (MQ=255) aaTAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGa > 2:134008/1‑148 (MQ=255) tAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGatat < 1:134008/149‑1 (MQ=255) ccATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATgg > 6:36662/1‑147 (MQ=255) aCTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGTCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCaa > 1:191665/1‑149 (MQ=255) ttGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCACGCG‑ACCAGCGCCTGCCCCAGTGCGCATTAAGCCCGCAGCCAGCCTGGGCAGAGCCCGAAAAGGCAGAACACTCGCAATCCTCCAtccctgac > 2:192088/1‑141 (MQ=255) gggCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTgg < 5:36662/149‑1 (MQ=255) gTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATgg < 6:146892/113‑1 (MQ=255) gTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATgg > 5:146892/1‑113 (MQ=255) ccGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCa < 7:8895/104‑1 (MQ=255) ccGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCa > 8:8895/1‑104 (MQ=255) acagcggcaggcGCTGCA‑CGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGt > 1:96831/6‑149 (MQ=255) gggAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTc > 6:91060/1‑149 (MQ=255) | TGTGCCGCAGGGAGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTC > NC_000913/3747032‑3747305 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |