Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,462,880 | 1 | . | C | 72.7% | 56.7 / 17.1 | 22 | intergenic (+234/‑132) | fadL/yfdF | long‑chain fatty acid outer membrane transporter/uncharacterized protein |
Reads supporting (aligned to +/- strand): ref base . (2/4); new base C (7/9); total (9/13) | |||||||||||
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. |
GCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CT‑‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAA > NC_000913/2462763‑2462993 | gCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CT‑‑‑CCCCCCCCTGGCTATTGTgcg < 1:243375/139‑1 (MQ=255) cTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CTCCCCCCCCCCCTGGCTATtgtg < 1:38334/139‑1 (MQ=255) aTAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CTC‑CCCCCCCCCTGGCTAtt < 2:308012/118‑1 (MQ=255) aTAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CTC‑CCCCCCCCCTGGCTAtt > 1:308012/1‑118 (MQ=255) cGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CTC‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGa < 2:475649/139‑1 (MQ=255) ggTAGCCACAGACAAAGCTGAAGGTATAGGCATATTTTCAAGCTTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCTGGATACCC‑‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGAt < 1:511240/139‑1 (MQ=255) ccACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CTC‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGtt < 2:499328/139‑1 (MQ=255) aGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CT‑‑ccccccc > 1:353012/1‑92 (MQ=255) aGACAAACCTGAAGGGAAAAGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCCGGATA‑CT‑‑ccccccc < 2:353012/92‑1 (MQ=255) aGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CT‑‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGAt < 1:496743/117‑1 (MQ=255) aGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CT‑‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGAt > 2:496743/1‑117 (MQ=255) aaaGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CTC‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGttttt < 2:372669/139‑1 (MQ=255) aaaGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CT‑‑ccccccccc < 1:219853/78‑1 (MQ=255) aaaGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CT‑‑ccccccccc > 2:219853/1‑78 (MQ=255) aGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CTCCCCCCCCCCCTGGCTATtgtg < 1:310602/89‑1 (MQ=255) aGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CTCCCCCCCCCCCTGGCTATtgtg > 2:310602/1‑89 (MQ=255) ttCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CTC‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCa < 1:312649/135‑1 (MQ=255) ttCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CTC‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCa > 2:312649/1‑135 (MQ=255) gTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATG‑CT‑‑‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGAt < 2:367447/85‑1 (MQ=255) gTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CT‑‑‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAAGTAAAGAt > 1:367447/1‑85 (MQ=255) tctcAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CTC‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATAttt > 1:107461/1‑139 (MQ=255) tACTTTTCTGCATTTCCAGGATA‑CTC‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGAt > 2:170333/1‑70 (MQ=255) tACTTTTCTGCATTTCCAGGATA‑CTC‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGAt < 1:170333/70‑1 (MQ=255) tACTTTTCTGCATTTCCAGGATA‑CTC‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAaatcaa > 1:18513/1‑139 (MQ=255) tttCCAGGATA‑CTC‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGAt < 2:335565/58‑1 (MQ=255) tttCCAGGATA‑CTC‑CCCCCCCCCTGGCTACTGTGCGCTCATACACGCAAATTAAAGAt > 1:335565/1‑58 (MQ=38) | GCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATA‑CT‑‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAA > NC_000913/2462763‑2462993 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |