| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 781,010 | 0 | C | . | 70.6% | 37.5 / 57.6 | 34 | intergenic (+92/‑58) | lysW/valZ | tRNA‑Lys/tRNA‑Val |
| Reads supporting (aligned to +/- strand): ref base C (5/5); new base . (13/11); total (18/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAG > NC_000913/780939‑781089 | tGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGA‑CCACCACTaa < 2:152454/81‑1 (MQ=255) gTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGA‑CCACCACTAAt > 2:152496/1‑81 (MQ=255) aaGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGA‑CCACCACTAATGAc > 1:152509/1‑81 (MQ=255) gtgtAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGA‑CCACCACTAATGACGGTg > 2:152505/1‑81 (MQ=255) tgtAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTg > 1:152494/1‑81 (MQ=255) gtAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATTCTGCACGA‑CCACCACTAATGACGGTggg > 1:152497/1‑81 (MQ=255) tAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGt < 1:152455/81‑1 (MQ=255) gATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTc < 1:152456/81‑1 (MQ=255) gATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCg > 2:152500/1‑81 (MQ=255) tAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGt < 2:152457/81‑1 (MQ=255) aCGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGTg < 2:152458/80‑1 (MQ=255) gTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGTgg > 1:152511/1‑79 (MQ=255) tCAGCAACGGCCCGTAGGGCGAGTGAAGCGAGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGTGGAAGTAGtt < 2:152459/81‑1 (MQ=255) cAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGtt > 1:152516/1‑81 (MQ=255) cAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGTGGAAGTAGttt > 2:152508/1‑81 (MQ=255) aGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTg < 1:152460/81‑1 (MQ=255) aaCGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTgtag > 1:152499/1‑81 (MQ=255) aCGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTgtagt < 1:152461/81‑1 (MQ=255) gCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTAt < 1:152462/81‑1 (MQ=255) cccGTAGGGCGAGCGAAGCGAGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATcc < 2:152463/81‑1 (MQ=255) gTAGGGCGAGCGAAGCGAGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAgc > 1:152515/1‑81 (MQ=255) cgagcgaAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGt < 1:152464/81‑1 (MQ=255) cgagcgaAGCGAGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTa < 1:152465/81‑1 (MQ=255) gagcgaAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTa > 1:152517/1‑81 (MQ=255) gagcgaAGCGAGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTAt > 2:152513/1‑81 (MQ=255) gaagcgaGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCggg > 2:152506/1‑81 (MQ=255) aagcgaGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGt > 2:152489/1‑81 (MQ=255) agcgaGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGt > 2:152527/1‑80 (MQ=255) agcgaGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTg < 1:152466/81‑1 (MQ=255) gcgaGTCATCCTGCACGA‑CCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGa < 1:152467/81‑1 (MQ=255) aTCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGc < 1:152468/80‑1 (MQ=255) aTCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCt > 1:152526/1‑81 (MQ=255) cACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTg > 1:152503/1‑81 (MQ=255) aCCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGgag < 2:152470/81‑1 (MQ=255) | TGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAG > NC_000913/780939‑781089 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |