| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 260,540 | 1 | . | G | 20.8% | 54.3 / 19.5 | 24 | coding (153/1104 nt) | proB | gamma‑glutamate kinase |
| Reads supporting (aligned to +/- strand): ref base . (10/9); new base G (3/2); total (13/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCATTCACGTCGG > NC_000913/260408‑260683 | ggtggtAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGAc > 4:43063/1‑149 (MQ=255) gtggtAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGAc < 3:43063/148‑1 (MQ=255) tAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATcgcgcccgac > 1:108250/1‑137 (MQ=255) tAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATcgcgcccgac < 2:108250/144‑8 (MQ=255) cACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGt > 3:115244/1‑148 (MQ=255) gtgtGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCg < 2:133038/149‑1 (MQ=255) aaCAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCAc > 4:137407/1‑130 (MQ=255) aCAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCTCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGc < 4:115244/148‑1 (MQ=255) aTCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACGtc < 5:173914/116‑2 (MQ=255) aTCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACGtc > 6:173914/1‑115 (MQ=255) cgcCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTGACCCGGAACTGCCAGCGACCATCGc > 8:78778/1‑149 (MQ=255) ccATAGCGTTGAACTTGTTCGCCAGTGCGCGCAGGTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTGACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTgg < 7:78778/149‑1 (MQ=255) ttGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGcggcgg > 6:28968/1‑147 (MQ=255) gAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGcc < 2:103202/123‑1 (MQ=255) gAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGcc > 1:103202/1‑123 (MQ=255) gCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGa > 2:67877/1‑149 (MQ=255) tGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGAt < 5:28968/149‑1 (MQ=255) aTGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACtgtg < 7:76877/145‑1 (MQ=255) aTGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACtgtg > 8:76877/1‑145 (MQ=255) gccgccGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGaa > 4:120696/1‑147 (MQ=255) ccgccgGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAg < 3:120696/149‑1 (MQ=255) ttATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATgg > 8:91103/1‑148 (MQ=255) caCGTCCGGCGGCGATCG‑CGCCCGACGTGAGCACCTCGTTTACCCGGAACTGCCAGCTACCATCGCCTCGCAACAACTGCTGGAGGCGGTAGGGCAGAGTCGAGAGATTCCACTGTGGGAACAGCTGTTTTCGATTTATGGCATTCACg > 6:41299/2‑149 (MQ=255) tCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCATTCACGTCgg < 7:91103/149‑1 (MQ=255) | GGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCATTCACGTCGG > NC_000913/260408‑260683 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |