Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,606,330 | 1 | . | G | 46.2% | ‑4.8 / 19.2 | 13 | noncoding (72/87 nt) | leuQ | tRNA‑Leu |
Reads supporting (aligned to +/- strand): ref base . (3/4); new base G (4/2); total (7/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.92e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGG‑G‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCT > NC_000913/4606106‑4606572 | cACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACAAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCACACGTCCGTAAGAACACTAACACCTGAAGCTAGCGAGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGAAATATCAAAGAGGTTTTTGGTGCG‑AGG‑G‑GGGGGGACTTTAACCCCCACGTCCGTacgg > 1:123500/1‑247 (MQ=39) gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCCTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGG‑G‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACAATAACACCTGAAGCTAGCGCGtt > 1:110980/1‑250 (MQ=255) caGTCATCTTACTTGTTTTGATAACACCTCGGTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACCAACACCTTAAGCTAGCGCGGCTACCAAGTCCGCCACCTCCGCCCGGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGG‑G‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCtt < 1:304433/251‑1 (MQ=37) tCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCACCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGG‑G‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgc > 1:205799/1‑251 (MQ=255) tCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCCCGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGG‑G‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACAACGGAGGCGCAt > 1:44743/1‑251 (MQ=255) gggttgACTCGAACCCCCACGTTAGTAAGGACAATAACACGTTAATATTGAGCGTCTACCAACTCCGACAGCTTAGCCTAGTGCGAGCAATAGCAACGTGGTTTTTGGTGCGTGGG‑G‑GGGGG‑ACTTGAACCCCCTCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATTCCATCAATTCTTAAAAATAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAg < 1:300065/246‑1 (MQ=18) aCGTCCGTAAGGGCACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGCTGCG‑AGG‑G‑GGGGG‑CCTTGAACCCCCACGTCCGTAAGGACACTAACATCTGGAGCTAGCGCGTATACCAATTACGCCACCTTCGCATACCCTCAATTGTTAAAAATAATTGCTACCACGGAGGCGCATTCTAGTGGTTGTCAGCTTTTCGACAATAGTTAAt > 2:115727/1‑251 (MQ=37) gTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTGCGAGCAATATCAACGAGGTTTTTGGTGCG‑AGG‑G‑GGGGG‑ACTTGAACCCCCCCGTCCGTCAGGACACTAACCCCTGCAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACTATCAATTCTTAAAAAGa > 1:212570/1‑194 (MQ=37) cgcgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAAGTGGTTTTTGGTGCG‑AGG‑G‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATg > 1:303667/1‑251 (MQ=255) gcgTCTACCAATTCCGCCACCTTCGCCCATTTCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGG‑GTGGGTG‑ACTTGATCCCCCACGTCCGTAAGGACACTAACTCCTTAAGCTAGCGCGTCTACCAATTCCTCCACCTTCGCATACCATCAATTCTTaaaaa < 1:302298/162‑1 (MQ=255) ccaagtctgacACCGTCGCCCAGTGCGAGCAATATCAGCGTGGTTATTGGTGCG‑GGGTG‑GGGGG‑ACTTGAGCCCCCACTTCCGTAAGTACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCTCCTTCGCATTCCATCAAGTCTTAAAAAGAATTGCTTCCACGGATGCGCAGACTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGAGAGAGGTGTAATTGCTGGAAAAATGGCCATCa < 1:30197/240‑1 (MQ=25) cAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGG‑G‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGa < 1:115727/251‑1 (MQ=255) ggTGCG‑AGG‑G‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgccagatg < 1:7452/251‑1 (MQ=255) ggg‑g‑ggggg‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCt < 1:168872/250‑1 (MQ=255) | CACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGG‑G‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCT > NC_000913/4606106‑4606572 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |