Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 2,000,857 | 1 | . | C | 30.8% | 19.9 / 15.4 | 13 | noncoding (535/2921 nt) | rrl3 | 23S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base . (7/2); new base C (4/0); total (11/2) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGCACACGGTTTCAGGTTCTATTTCACTCCCC‑T‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAG > CP000730/2000825‑2000996 | cgtacacGGTTTCAGGTTCTATTTCACTCCCC‑TCACAGGGGTTCTTTTCGCCTTTCCCTCACGGTACTGGTTCACTATCGGTCAGTCAG‑TAGTATTTAGCCTT‑GGAGGATGGTCCCCCCATATTCAGACAGGATATCACGTGTCCcgc > 1:120805/4‑147 (MQ=2) cacaCGGTTTCAGGTTCTATTTCACTCCCC‑T‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCcgtcgt > 1:317030/1‑148 (MQ=32) cacGGTTTCAGGTTCTATTTCACTCCCC‑T‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTAc > 1:205965/1‑148 (MQ=32) cacGGTTTCAGGTTCTATTTCACTCCCC‑T‑TCCGGGGGGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTAc > 2:120263/1‑148 (MQ=17) ggTTTCAGGTTCTATTTCACTCCCC‑TCACAGGGGTTCTTTTCGCCTTTCCCTCACGGTACTGGTTCACTATCGGTCAGT‑CAGTAGTATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGATATCACTTGTCCCGCCCTACTc > 1:317467/1‑148 (MQ=2) ggTTTCAGGTTCTATTTCACTCCCC‑TCACAGGGGTTCTTTTCGCCTTTCCCTCACGGTACTGGTTCACTATCGGTCAGT‑CAGTAGTATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGATATCACGTGTCCCGCCCTACTc > 1:205522/1‑148 (MQ=2) ttCAGGTTCTATTTCACTCCCC‑T‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGa < 2:292042/148‑1 (MQ=32) gTTCTATTTCACTCCCC‑T‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCAc < 1:345683/148‑1 (MQ=32) tATTTCACTCCCCAC‑ACAGGGGTTCTTTTCGCCTTTCCCTCACGGTACTGGTTCACTATCGGTCAGT‑CAGTAGTATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGCTATCACGTGTCCCGCCCTACTcgatttcactgaa > 2:116113/1‑136 (MQ=2) tATTTCACTCCCC‑TCACAGGGGTTCTTTTCGCCTTTCCCTCACGGTACTGGTTCACTATCGGTCAGT‑CAGTAGTATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGATATCACGTGTCCCGCCCTACTcgatttcactgaa > 2:224754/1‑136 (MQ=2) ttCACTCCCC‑T‑TCCGGGGTGCTTTTCACCTTTCCTTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAagag > 2:273728/1‑148 (MQ=17) cACTCCCC‑T‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAagagag > 2:344064/1‑148 (MQ=32) cACTCCCC‑T‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAagagag > 2:96601/1‑148 (MQ=32) | AGCACACGGTTTCAGGTTCTATTTCACTCCCC‑T‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAG > CP000730/2000825‑2000996 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |