Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 4,519,968 | 0 | C | T | 42.1% | 4.3 / 17.0 | 28 | intergenic (+46/‑36) | purB/roxA | Adenylosuccinate lyase/50S ribosomal protein L16 3‑hydroxylase |
Reads supporting (aligned to +/- strand): ref base C (1/10); new base T (8/0); total (18/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.19e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGC‑GT‑TTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATC‑TCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACGCC > ACB143/4519841‑4520097 | aGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGC‑GT‑TTTTATTCCCGTACGaa < 8:124328/137‑1 (MQ=255) aaGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGGT‑TTTTTTTGCCGTA‑GAAAAttt > 3:155104/1‑136 (MQ=255) aaGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAAAGCCCCGCCTGGCCGGGCGGT‑TTTTTTTGCCGTA‑GAAAAttt > 2:84627/1‑136 (MQ=255) ggCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGC‑GT‑TTTTATTCCCGTACGAAAATTACGttttt < 3:354386/137‑1 (MQ=255) ccGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGC‑GT‑TTTTATTCCCGTACGAAAATTACGTTTTTTc < 6:364354/137‑1 (MQ=255) cGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGC‑GT‑TTTTTTTTCCGTAAGGAAATTAAGttttttga > 6:112183/1‑135 (MQ=255) cTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGAGTACAACGCCCGGCCAGGCCGGGC‑GT‑GTTTTTTTCCGTGAGAAAATTACGTTTTTTAAATat > 6:183737/1‑136 (MQ=255) tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGGCGGGC‑GT‑TTTTATTTCCCGACGAAAATTAAGTTTTTTTaaaaag > 7:369768/1‑133 (MQ=255) tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGC‑GT‑TTTTTTTACCCGACGAAAAATAAGGTTTTTGAATatg > 2:39760/1‑135 (MQ=255) tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGC‑GT‑TTTTATTCCCGTACGAAAATTACGTTTTTTCAATAgg < 6:204073/137‑1 (MQ=255) tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCGGGCCGGGC‑GG‑TGTTTTTTCCGGAGCGAAATTAAGGTTTTTTAATTgg > 1:258221/1‑137 (MQ=255) tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTAAAAAGCCCGGCCGGGGCGGGC‑GG‑TTTTTTTTCCGGAACAAAATTTAATTTTTTGAATatg > 7:338550/1‑135 (MQ=255) tCAAGCAGCTGACCCCCGCTACCAACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACCCCCGGGCTGGCCGGGC‑GT‑GTTTTTTAAGGGACGAAAATTAGGTTTTTTCaaacag > 6:68543/1‑133 (MQ=255) aCCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCGGGCCGGGC‑GG‑TTTTTTTTCCGGACGCAAAATACGGTTTTTTAATATGTTGGACATGaa > 7:243084/1‑137 (MQ=255) ccccGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGC‑GT‑TTTTATTACCGGACGAAAATTAAGGTTTTTAAATATGTTTAAAAAGAAta > 7:107524/1‑136 (MQ=255) cccGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGC‑GT‑TTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATcc < 3:223479/137‑1 (MQ=255) cGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGC‑GT‑TTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTg < 3:108528/137‑1 (MQ=255) cctacATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGC‑GT‑GTTTTTTTCCGGACGCAAATTTAGGTTTTTTAAAATGTTGTACAAGAATACGGATaaa > 4:223479/1‑135 (MQ=255) cGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGC‑GT‑TTTTATTACCGTTCGAAAATTAAGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGc > 2:202753/1‑137 (MQ=255) ggTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGC‑GT‑TTTTTTTGCCGGTCGAAAATTAAGGTTTTTGAATTTGTTGAACAAGAATACGGATACACCACTGcc > 4:174192/1‑136 (MQ=255) ccAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGC‑GT‑TTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCa < 8:303247/137‑1 (MQ=255) aaaCGCATCTAAGCTGAGCGTCGCGTACAACGCCCCGCCTGGCCGGGC‑GG‑TTTTTTTTCCGGACGAAAATTTAGTTTTTTGAAAAAGGTTTACATGAATCCGTTAAAACCACAGCCCCTGCTCGGCGGCAtcgaggc > 8:226636/1‑132 (MQ=255) gCTGAGCGTCGAGTACAACGCCCGGCCTGGCCGGGA‑GT‑TTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATC‑TCGGCCCTCGAATTCAt < 3:312102/137‑1 (MQ=255) gTCGCGTACAACGCCCGGCCTGGCCGGGC‑GTGTTTTTTAGCGGAACCAAAATAAGGTTTTTT‑AATTGGTTGGAAAAGAATACGGCTAAGACCCATCAGCAGCGCGGCGGCATCATCGGCCC‑CGGATTGATTcgcggc > 8:106476/1‑135 (MQ=255) gCCTGGCCGGGC‑GT‑TTTTTTTCCCGTAAGAAAATTTAGGTTTTTCAAAAGGTTGGAAATTAATCCTGCTACTCCACAGAAGCTGCTCGGCGGGATC‑ACGGCCCGCCAATTCAAGCGGGCAATCTGGCAGAATAAGcc > 6:162758/1‑137 (MQ=255) gCCTGGCCGGGC‑GT‑TTTTATTGCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGTATCCTGATACTCCACTGCAGCTGCTCGGCGGCATC‑TCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGcc > 6:335428/1‑137 (MQ=255) gCCGGGC‑GT‑TTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATC‑TCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCgctgc < 4:155104/137‑1 (MQ=255) ttttATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATC‑TCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACGcc < 1:184467/137‑1 (MQ=255) | AGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGC‑GT‑TTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATC‑TCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACGCC > ACB143/4519841‑4520097 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |