Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,622,816 | 0 | G | A | 42.5% | 40.9 / 36.4 | 47 | intergenic (+75/+144) | yneM/dgcZ | inner membrane‑associated protein/diguanylate cyclase, zinc‑sensing |
Reads supporting (aligned to +/- strand): ref base G (25/2); new base A (0/20); total (25/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.37e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGG > NC_000913/1622667‑1622956 | tttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCg < 2:433622/150‑1 (MQ=255) tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGaa > 1:782817/1‑150 (MQ=255) tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGaa > 1:53762/1‑150 (MQ=255) aTGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGaaaa < 1:533969/151‑1 (MQ=255) gCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGaagaaaa < 2:756372/151‑1 (MQ=255) gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAAcc > 1:590582/1‑151 (MQ=255) tACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCg > 2:606717/1‑151 (MQ=255) gggAATAATTTTATTTTTTGGTTTTGTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGcc > 2:88439/1‑149 (MQ=255) attttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcg > 2:249265/1‑151 (MQ=255) ttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGCTGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggcg > 2:795489/1‑151 (MQ=255) ttttttttggtttcctggCCGCTTTTTCCACCCAAAAAGGGGATAAAAAAGGAACGGAAAAAACCCCCCACCTACCCGCCCCT‑TTTTA‑ACAGTTGTGAAAAGGGGGCT‑TGATTTTTAAGAGGGCAAAAAAAAACCGCCAGTAAACcggcg < 2:461415/148‑1 (MQ=255) ttttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGACTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTg > 2:548640/1‑151 (MQ=255) tttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTg > 2:317989/1‑150 (MQ=255) tttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTg > 2:478238/1‑150 (MQ=255) ttttctggCCGCTTATTACAGCAAAAAGGGGGATAACTAA‑AAAGAGAAAAAAACCCCCCCCCAACCCGCCCCTTTTTT‑AAATTTGTGAAAAGGGGGCT‑TGTTTTTTAGGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt < 2:550201/150‑1 (MQ=255) ttttcgggCGGTGTATTCCAGCCACAAAGGGGATGATTAATGAACGGAAAAAACCCCCCCCCTACCCGCCCCT‑TTGTA‑AAAGTTGTATACAGGGGGCT‑TGTTTTTTAGGGGGGGAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt < 1:225691/150‑1 (MQ=255) ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc > 1:655633/1‑150 (MQ=255) ctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg > 1:725995/1‑151 (MQ=255) ggcggcGTATTTCACCCACAAAGGGGATGAAAAAGGAACGAAAAAACCCCCCCCCCTACCCGGCCCT‑TTGTT‑ACAGTTGTGTACAGGGGGCT‑TGTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGa < 2:13572/151‑1 (MQ=255) ggCCGCGTTTTACAGCCAAAAAGGGGATGACTAATGAAGGGAAATAACCCCTCCCCTACCCGGCCCT‑TTGTT‑ACAGTTGTGTAAAGGGGGCT‑TGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg < 1:611192/149‑1 (MQ=255) ccGCTTTTTTCAGCCAAAAAGGGGATAACAAATAAAGGGAAATAACCCCCCCCCAACCCGCCCCT‑TTTTT‑ACGTTTGTGTACAGGGGGCT‑TTTTTTTTAGGCCGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt < 1:637521/150‑1 (MQ=255) tttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACCAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAGAGATTtgt > 1:113236/1‑151 (MQ=255) ttccAGCCAAAAGTGGGATGAAAAATGAAGGGAAATACCCCCTCCCCTACCCGCCCCT‑TTTTT‑ACAGTTGTGAAAAGGGGGCT‑TGTTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt < 2:612800/146‑1 (MQ=255) aaaGGGAATAACAAATGAAGGAAAAAAACCCCCCCCCTACCCGCCCCT‑TTTTT‑ACATTTGTGTAAAGGGGGCT‑TTTTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTg < 2:793446/147‑1 (MQ=255) aaTGGGATGATTAATGAAGGGAAATATTCCCTCCCCGAGCCGCCCCT‑TTTTT‑AGATTGTTGCAAGGGGGGCT‑TTTTTTTTAGGCCGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt < 1:161521/151‑1 (MQ=255) ggATAACTAATGAAGGGAAATAACCCCTCCCCTACCCGCCCCT‑TTGTT‑ACAGTTGTGTACAGGGGGCT‑TGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGc < 1:422080/151‑1 (MQ=255) cTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc > 1:184233/1‑149 (MQ=255) tAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAgg > 1:49708/1‑151 (MQ=255) ggAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtttt > 1:286386/1‑149 (MQ=255) gaTAATCCCTCACCTATCCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc > 1:808933/1‑151 (MQ=255) ccccccccTACCCGCCCCT‑TTTTT‑ACAGTTGT‑TAAAAGGGGCCTTTTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc < 2:95506/139‑1 (MQ=255) acccGCCCCT‑TTTTT‑ACATTTGTGAAAAGGGGGCC‑TGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCg < 2:263828/149‑1 (MQ=255) aCCCGCCCC‑TTTTTAACTGTTGTTAAAAGGGGGCT‑TTTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAAc < 2:768841/145‑1 (MQ=255) c‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGacac > 2:582028/1‑149 (MQ=255) ttGTT‑ACAGTTGTGAACAGGGGGCT‑TGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAg < 1:548695/150‑1 (MQ=255) tGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGt > 2:240992/1‑150 (MQ=255) tt‑ACAGTTGTTTAAAGGGGGCC‑TGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAg < 1:699298/150‑1 (MQ=255) tt‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAg > 2:407824/1‑150 (MQ=255) tgtgTAAAAGGGGCT‑TGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGt < 2:782319/150‑1 (MQ=255) aCAAGGGGCG‑TGTTTTTTAGGACGGAAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc < 2:26546/151‑1 (MQ=255) aaaggggcgttgttttttAGGACGGTAAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct < 2:725924/137‑1 (MQ=255) aaGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa > 2:332238/1‑150 (MQ=255) agggggctttttttttAGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGttt < 1:36506/134‑1 (MQ=255) ggCC‑TGATTTTTATGACGGCGAGAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga > 1:651378/1‑150 (MQ=255) gCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaa > 1:786770/1‑150 (MQ=255) cttgttttttATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaa < 2:61244/144‑1 (MQ=255) tGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt > 2:507548/1‑151 (MQ=255) ttATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACgg > 1:884306/1‑151 (MQ=255) | TTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGG > NC_000913/1622667‑1622956 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |