Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,622,816 | 0 | G | A | 43.5% | 12.3 / 19.9 | 23 | intergenic (+75/+144) | yneM/dgcZ | inner membrane‑associated protein/diguanylate cyclase, zinc‑sensing |
Reads supporting (aligned to +/- strand): ref base G (13/0); new base A (0/10); total (13/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.74e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAA‑GGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTT > NC_000913/1622671‑1622950 | tGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAA‑GGGGCCTGATTTTTATGACGGCGaaaa > 2:55327/1‑150 (MQ=255) ataatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAA‑GGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTaa > 1:144484/1‑150 (MQ=255) ataatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAA‑GGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTaa > 1:380147/1‑150 (MQ=255) tttttttctgtttttcgggCCGCTTTTTCCAGCCAAAAATGGGATGATAAATGAACGGAGATAACCCCTCCCCTACCCGCCCCT‑TTGTTACATTTGTGTACAG‑GGGGCTTGTTTTTTAGGACGGCAAAAAAAAACCGCCAGTAAACcggc < 2:604968/147‑1 (MQ=255) ttttttgtttttcgggCCGTATTTTACAGCCAAAAAGGGGATAACTAATGAAGGGAAATAACCCCCCCCCAACCCGCCCCT‑TTGTTACTGTGGTGTAAAG‑GGGGCTTTTTTTTTATGCCGCAAAAAAAAAACCGCCAGTAAACCGGCGGt < 1:260428/150‑1 (MQ=255) ttttctggCCGCTTTTTTCAGCCAAAAAGGGGATGACAAATGAACGGAAATAACCCCTCCCCTACCCGCCCCT‑TTGTTACAGTTGTGTAAAA‑GGGGCTTTTTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAAt < 2:543073/146‑1 (MQ=255) ttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAA‑GGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt > 1:271503/1‑150 (MQ=255) tgttctgcCCGCTTTTTCCACCCAAAAAGGGAATGATAAAGGAAGGGAAATACCCCCCCCCCAACCCGCCCCT‑TTTTTACATTTGTGAACAG‑GGGGCTTTTTTTTTAGGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCt < 1:238075/147‑1 (MQ=255) ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAA‑GGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc > 1:155121/1‑150 (MQ=255) ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAA‑GGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc > 1:539197/1‑150 (MQ=255) tctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAA‑GGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa > 2:21634/1‑150 (MQ=255) ccGCGTTTTCCACCCAAAAAGGGGATAATTAATGAAGGGAAATACCCCCTCACCAACCCCCCCCT‑TTTTTACAGTGTTGAAAAG‑GGGGCGTTTTTTTTAGGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt < 1:129594/150‑1 (MQ=255) aGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAA‑GGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGttt > 2:187209/1‑151 (MQ=255) aaaGGGGATGATTAAGGAACGGAAATACCCCCTCACCTACCCGCCCCT‑TTTTTACGGTTGTGAAAAG‑GGGGCTTGATTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGttt < 1:28627/145‑1 (MQ=255) aagggAGATAACCCCCCCCCAACCCGGCCCT‑TTTTTACAGTTGTTTACAG‑GGGGCTTTTTTTTTAGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAt < 2:256653/146‑1 (MQ=255) accccccccccTACCCGCCCCT‑TTTTTACAGTTGTGTAAAG‑GGGGCTTGTTTTTTAGGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACt < 2:485158/142‑1 (MQ=255) aTCCCCCCCCCAACCCGCCCCTTTTTTACA‑TTGTGGACAAGGGGGCTTTTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTcc < 1:74584/145‑1 (MQ=255) aCCTAACCGGCCCC‑TTGTTACAGTTGTGTACAA‑GGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAAACGGCGGGGAATGCTTGCATGGGTAGATTTGTGTTTTTCTTTTACGCTAACAGGCATTTTCCTGCACTGATaa > 1:212352/1‑148 (MQ=255) ttGTTACAGTTGTGTACAG‑GGGGCTTGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAg < 1:55321/150‑1 (MQ=255) ttGTTACAGTTGTGTACAA‑GGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGt > 2:159228/1‑151 (MQ=255) aCAGTTGTGTACAA‑GGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt > 1:347463/1‑151 (MQ=255) tACAA‑GGAGCCTGATTTTTATGACGGCGAAAAAAAACCGGCAGGAAACCGGCGGTGAATGGTTGCATGGCTGGATTTAAGTTTTGCTTTTACGCCAACAGGCATTTTCCTGCACTGATAACGAATCGTTTACAAAGTAGCCTCAGTTTtc > 1:184226/1‑150 (MQ=255) tGATTTTTATGAAGGAGAAAAAAAAAGGCCAGTAAACCGGGGGTGAATGTTTGGATGGATAGATTTGTGTTTTGGTTTTACGATAAAAGGGATTTTTCTGGAATGATAACGAATCGTTGACACAGTAGCATCAGTTTTTTTAATGAATGtt > 1:66904/1‑151 (MQ=255) | TGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAA‑GGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTT > NC_000913/1622671‑1622950 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |