| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 2,000,908 | 0 | C | G | 20.8% | 54.4 / 11.5 | 24 | noncoding (484/2921 nt) | rrl3 | 23S ribosomal RNA |
| Reads supporting (aligned to +/- strand): ref base C (11/8); new base G (5/0); total (16/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.30e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGCACACGGTTTCAGGTTCTATTTCACTCCCC‑T‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATTACCTTCTTTGATTCATCTT > CP000730/2000825‑2001041 | 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:96601/1‑148 (MQ=32) cACTCCCC‑T‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAagagag > 2:344064/1‑148 (MQ=32) ggggTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTCGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCg < 1:267986/148‑1 (MQ=17) gTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACt > 2:353091/1‑148 (MQ=32) ttttCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAg < 2:349418/148‑1 (MQ=32) ttttCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAg < 1:251607/148‑1 (MQ=32) cACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGatt < 1:367576/148‑1 (MQ=32) cACCTTTCCCTCACGGTACTGCTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGatt < 1:120014/148‑1 (MQ=17) aCCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGatta > 1:298222/1‑148 (MQ=32) aCCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTGCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGattt > 1:287323/1‑147 (MQ=17) tCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGCTCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATTACcttct > 2:204200/1‑148 (MQ=17) cGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATTACCTTCTTTg < 1:335584/148‑1 (MQ=32) ttCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATTACCTTCTTTGATTCATCtt > 2:91070/1‑148 (MQ=32) | AGCACACGGTTTCAGGTTCTATTTCACTCCCC‑T‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATTACCTTCTTTGATTCATCTT > CP000730/2000825‑2001041 |
| 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: ‑ |