| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,622,802 | 0 | A | T | 43.6% | 75.8 / 28.4 | 55 | intergenic (+61/+158) | yneM/dgcZ | inner membrane‑associated protein/diguanylate cyclase, zinc‑sensing |
| Reads supporting (aligned to +/- strand): ref base A (26/5); new base T (0/24); total (26/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.81e-11 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ATGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTT > NC_000913/1622658‑1622950 | atGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGAttttt < 1:685233/150‑1 (MQ=255)atGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTa < 2:993025/151‑1 (MQ=255) gAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTAt < 1:328946/150‑1 (MQ=255) aTGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGAc > 2:98608/1‑151 (MQ=255) ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGa > 2:321201/1‑150 (MQ=255) tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGaa < 2:934799/150‑1 (MQ=255) tGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGaaaa > 1:880002/1‑150 (MQ=255) ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGaaaaa > 2:439413/1‑150 (MQ=255) aCTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGc > 1:795101/1‑151 (MQ=255) aCTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGc < 2:241339/151‑1 (MQ=255) tGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGcc > 1:914384/1‑150 (MQ=255) gAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGt > 1:14085/1‑150 (MQ=255) aatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAAcc > 2:902765/1‑151 (MQ=255) ttttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcgg > 2:430642/1‑151 (MQ=255) ttttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcgg > 1:571449/1‑151 (MQ=255) tctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGaa > 1:6318/1‑150 (MQ=255) ctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATg > 1:431990/1‑151 (MQ=255) ggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCt > 1:164161/1‑151 (MQ=255) tttttcgggCCGCTTATTCCAGCAAAAAGGGGGATTATTAATGAGCGAAAATAACCCCCCACCAACCCGCCCC‑T‑TTG‑TTACATTTGTGAAAAGGGGGCT‑TGTTTTTTAGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCt < 1:607493/149‑1 (MQ=255) ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc > 2:74998/1‑150 (MQ=255) ggCCGCTTTTTCCACCCAAAAAGGGGATGACTAATGAAGGGAAATAACCCCTCCCCAACCCGCCCC‑T‑TTT‑TTACATTTGTGTAAAGGGGGCT‑TGTTTTTTATGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa < 1:988182/147‑1 (MQ=255) gcccgattattacgcccACAAATGGTATTAATAATGAGGGGAAAAATCCCCTCACCAACCCGGCCC‑T‑TTT‑TTACATTTGTGAAAAGGGGGCG‑TTTTTTTTATGAGGGGAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATgg < 2:478673/135‑1 (MQ=255) gTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt > 1:105372/1‑151 (MQ=255) tATTTCACCCAAAAAGGGGATGACTAATGAAGGGAAATAACCCCTCCCCTACCGGCCCC‑C‑TTT‑TTAAGTTTGTGTAAAGGGGGCC‑TGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt < 1:134007/150‑1 (MQ=255) aaaaGGGGATAACTAATGAAGGGAAATACCCCCCCCCCAACCCGGCCC‑T‑TTG‑TTACGGTTGTGTAAAGGGGGCC‑TGTTTTTTAGGCGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt < 1:203570/150‑1 (MQ=255) aaaTGGGAGGATAAATGAACGAAGAT‑ATCCCCCCCCCTACCGGCCC‑CTTTT‑TTACATTTGTGAACAGGGGCCT‑TGTTTTTTATAAGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCttt < 1:433603/151‑1 (MQ=255) gggATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGTTTTTTATGACGGCGAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACg < 2:224697/151‑1 (MQ=255) gggATGACTAAGGACGGGAAATAACCCCTCACCTCCCCGCCCC‑T‑TTT‑TTACATTTTTAAACAAGGGGTT‑TTTTTTTTGTGGCGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACg < 2:940419/151‑1 (MQ=255) tgaaaaaTGAAGGGAAATACTCCCTC‑CCAAACCCGCCC‑C‑TTGTTTATAGTGGTGCAGAGGGGGCC‑TATTTTTTATGGCGGAGAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCt < 1:813292/144‑1 (MQ=255) aCTAATGAGGGGAAATATCCCCCCCCATACCCGCCCC‑T‑TTT‑TTACTTTTGTCAAAAGGGGGCT‑TGTTTTTTATGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACa < 2:819459/151‑1 (MQ=255) aCTAATGAAGGGAAAAAACCCCTCCCCTACCCGGCCC‑T‑TTG‑TTACATTTGTGTACAGGGGGCT‑TGTTTTTTATGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACa < 1:826876/151‑1 (MQ=255) aaTGAGGGGAGATACCCCCTC‑CCCAACCGGCCCTC‑TTT‑TTAGATTTGTGAACAGGGGGCT‑TTTTTTTTATGCCGGGAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAg < 1:821586/149‑1 (MQ=255) aaTAAACGGAGATACCCCCTCCCCAACCCGCCCC‑T‑TTG‑TTACTGTTGTGTAAAAGGGGCT‑TGTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAg < 1:522211/149‑1 (MQ=255) aCGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtttt > 1:386770/1‑151 (MQ=255) gaTAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc > 2:282298/1‑151 (MQ=255) ccccccccaacccGCCCC‑T‑TTG‑TTAAATTTGTTAAAAGGGGGCT‑TTTTTTTTAGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc < 1:530029/136‑1 (MQ=255) cccTCCCCTACCCGCCCC‑T‑TTT‑TTACATTTGTGTACAGGGGGCT‑TTTTTTTTATGCCGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACt < 1:36973/148‑1 (MQ=255) ccTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGGGAATGGCTGCATGGGTAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTTCTGCCCt > 2:384888/1‑147 (MQ=255) ccccTACCCGGCCC‑T‑TTG‑TTACATTTGTGTACAGGGGGCT‑TGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATaa < 1:201267/147‑1 (MQ=255) aCCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACg > 1:945470/1‑150 (MQ=255) cGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg > 1:515935/1‑151 (MQ=255) ccc‑t‑TTG‑TTACAGTTGTGTACAGGGGGCT‑TTTTTTTTATGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGa < 1:141585/149‑1 (MQ=255) cccttggttaagtttGTGTA‑AAGGGGCCTTGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg < 2:838026/135‑1 (MQ=255) ttacttttGTGAAAAGGGGGCT‑TTTTTTTTAGGCCGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTa < 1:844999/143‑1 (MQ=255) ttacgtttGTGTAAAGGGGGCT‑TGTTTTTTAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAg < 1:846791/144‑1 (MQ=255) ttacatgtgggaaCAGGGGGCC‑TTTTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTa < 1:539115/137‑1 (MQ=255) taactgTTGTGTA‑AAGGGGCCTTGTTTTTTAAGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAg < 1:902875/145‑1 (MQ=255) atACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGCGAGATTTGGGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTa > 2:185462/2‑149 (MQ=255) aCAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt > 2:132591/1‑151 (MQ=255) aCAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt > 1:298939/1‑151 (MQ=255) aCAGTTGTGAACAGGGGGCT‑TTTTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt < 1:255410/151‑1 (MQ=255) ttGTGAAAAGGGGGCC‑TTTTTTTTATGGCGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGc < 1:264236/145‑1 (MQ=255) tgtACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtttt > 2:868784/1‑151 (MQ=255) tACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct > 1:477920/1‑151 (MQ=255) tGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt > 1:850904/1‑151 (MQ=255) | ATGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCC‑C‑TTG‑TTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTT > NC_000913/1622658‑1622950 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |