Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,287,220 | 0 | A | G | 40.7% | 21.3 / 26.8 | 28 | noncoding (17/171 nt) | rttR | rtT sRNA, processed from tyrT transcript |
Reads supporting (aligned to +/- strand): ref base A (4/12); new base G (7/4); total (12/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.08e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTT > NC_000913/1287091‑1287345 | cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCa > 2:714683/1‑136 (MQ=255) cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCa < 1:714683/136‑1 (MQ=255) cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAggg < 2:3959/139‑1 (MQ=255) gccGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGa > 2:709237/3‑139 (MQ=11) aaccttaatCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGaa > 1:273367/9‑84 (MQ=255) aaCCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGaa < 2:273367/84‑1 (MQ=255) aatCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGAACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCt > 1:479332/3‑139 (MQ=14) atCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTg < 2:479332/138‑1 (MQ=14) gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcgttgcctgcg < 2:101688‑M2/139‑63 (MQ=255) ccTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcgttgcctgcggcaacgctctct < 2:394493‑M2/139‑75 (MQ=255) tCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCgg > 1:29800/1‑139 (MQ=18) ccccGCATGGGCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCggg > 2:120446/1‑139 (MQ=18) ccccGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcgttgcctgcggcaacgctctctcgct < 1:100578‑M2/139‑79 (MQ=255) ccccGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcgttgcctgcggcaacgctctctcgct < 1:639048‑M2/139‑79 (MQ=255) cGCATGGGCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAGCCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAc < 1:120446/139‑1 (MQ=18) cGCATGGGCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAc < 2:566558/139‑1 (MQ=21) ggCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGAc > 2:237293/1‑94 (MQ=11) ggCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGAc < 1:237293/94‑1 (MQ=11) gCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCccacca > 2:586083/1‑139 (MQ=18) gCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcgttgcctgcggcaacgctctctcgctggcgctcgag < 1:391249‑M2/139‑89 (MQ=255) ttgaTTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAAtt > 2:461199/4‑139 (MQ=17) ttgaTTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAAtt > 2:697032/4‑139 (MQ=18) tGCGGATTCG‑TTGGGAAGTTTAGGGACTTTTGAAAGTgatggtggtgggggaaggattcctcagttctgcgcgcttcgcccttcgggtagttacctgcgacaacgctctctcgctggcgctcgagccgaaccttaatcg < 2:568327‑M2/139‑103 (MQ=255) tGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcgttgcctgcggcaacgctctctcgctggcgctcgagccgaaccttaatcg < 1:676214‑M2/139‑103 (MQ=255) cGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcgttgcctgcggcaacgctctctcgctggcgctcgagccgaaccttaatcgaa > 2:391249‑M2/1‑35 (MQ=255) cGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcgttgcctgcggcaacgctctctcgctggcgctcgagccgaaccttaatcgaa < 1:620776‑M2/139‑105 (MQ=255) tCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcg < 1:684154‑M2/83‑54 (MQ=255) tCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcg > 2:684154‑M2/1‑30 (MQ=255) | CGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTT > NC_000913/1287091‑1287345 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |