| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 598,486 | 0 | T | C | 20.0% | 55.9 / 10.7 | 30 | intergenic (+35/+247) | tufA/yhaA | elongation factor EF1A/M20D subfamily peptidase |
| Reads supporting (aligned to +/- strand): ref base T (9/15); new base C (1/5); total (10/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.33e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTAGAATTAATCGCTCCAATCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTTAAATAAAATTGT > CP000730/598347‑598613 | gTAGAATTAATCGCTCCAATCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATc > 2:186424/1‑141 (MQ=255) aaTCGCTCCAATCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATAATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGACCCCTCGAtt < 2:342018/141‑2 (MQ=255) cctcAATCGCGTTCGAATACGGTACTCGTCTCTCAATCCGTGAAAGAAGACTTACTAAAAGACAACCCGATGTTACTGAAATCATTAAATAATTTCTAATTTCATAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAgag > 2:220619/4‑139 (MQ=255) cAATCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTc < 1:143275/141‑1 (MQ=255) cAATCGCGATTGAACACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGGGGGGTc < 2:285997/141‑1 (MQ=255) aaTCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTATTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCt < 2:128679/141‑1 (MQ=255) aaTCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCt < 2:165187/141‑1 (MQ=255) aaTCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCt < 1:92540/141‑1 (MQ=255) tCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCtt‑t < 1:84626/141‑1 (MQ=255) cgcgATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCtt‑tt > 2:180589/1‑141 (MQ=255) ttGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAAAGACCCCTCGATTGAGGGATCtt‑cttttt < 2:270908/141‑2 (MQ=38) gAAGACGGTAGGCGTTTGTCAATCCGGGAAAGTGGATATACTGTATAATCATGTATTTTTACTTAAATCATTAAATAATTTCTAATTTCTTATATTTTATAT‑AAAAA‑AAGACCCCTCGATTGAGGGATCTTCTTTTTatat < 2:315686/141‑4 (MQ=12) aGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAA‑AAGACCCCTCGATTGAGGGATCTTCTTTTTatataa < 1:138102/141‑6 (MQ=37) ggTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTaaa > 2:307679/1‑141 (MQ=255) tCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTAtt > 1:45420/1‑141 (MQ=255) ccGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTc < 1:186424/141‑1 (MQ=255) gACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTagaaga > 1:257551/1‑141 (MQ=255) gACGTACTGAAGGATCAGGCGATGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCACAATCGAGGGGTCTT‑TTTTTAAAGGGTAAACTTTGTAATGACTATTGGATTTtgaaga > 2:310719/1‑141 (MQ=255) aCGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAgaagaa > 1:219968/1‑141 (MQ=255) cGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAACAAAAATCCATAAATCGAGGGGCCTT‑TTTTTAATGTGTACATTTTGTAATAGCTATTCGATTTAGAAGAAc < 1:50520/141‑1 (MQ=255) tGTAGGATCTGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCGATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaataa < 2:137274/141‑1 (MQ=255) tGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaataa < 1:153178/141‑1 (MQ=255) gATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgat < 1:307679/141‑1 (MQ=255) gCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGACCCCTcgattg < 1:32421/71‑6 (MQ=255) gCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGACCCCTcgattg > 2:32421/1‑66 (MQ=255) cATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACtt < 1:214411/141‑1 (MQ=255) tttctaatttctTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGata < 2:256731/141‑1 (MQ=255) ttctTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGtt > 2:95003/1‑141 (MQ=255) tAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTtaaa < 1:180589/141‑1 (MQ=255) tat‑aAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTTAAATAAAAttgt < 1:115166/141‑1 (MQ=255) | GTAGAATTAATCGCTCCAATCGCGATTGAAGACGGTACTCGTTTCTCAATCCGTGAAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATAT‑AAAAAGAAGATCCCTCAATCGAGGGGTCTT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGATAATACTGTTTAAATAAAATTGT > CP000730/598347‑598613 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |