| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,622,816 | 0 | G | A | 44.4% | 41.4 / 50.6 | 54 | intergenic (+75/+144) | yneM/dgcZ | inner membrane‑associated protein/diguanylate cyclase, zinc‑sensing |
| Reads supporting (aligned to +/- strand): ref base G (28/2); new base A (0/24); total (28/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.32e-13 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTT > NC_000913/1622667‑1622961 | tttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑g > 1:485507/1‑150 (MQ=255) ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑Ga < 1:257316/150‑1 (MQ=255) tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑Gaa > 2:998826/1‑150 (MQ=255) tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑Gaa > 2:555279/1‑150 (MQ=255) ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑Gaaaaa > 2:581613/1‑150 (MQ=255) gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAAcc > 1:702818/1‑151 (MQ=255) gggAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAg < 1:948523/151‑1 (MQ=255) aataatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTa > 1:914251/1‑150 (MQ=255) ataatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTaa > 2:363325/1‑150 (MQ=255) taatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAAc > 2:202358/1‑151 (MQ=255) ttttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACcgg > 1:425852/1‑151 (MQ=255) tttttttttgtttttcgggCCGCTTTTTCCAGCCACAAAGGGGATGACAAATGAGGGGAAAAAA‑CCCCCCACCTACCCGGCCCTTTGTT‑AAAG‑TTGTGAAAAGGGGGCT‑TGTTTTTTATGCGGGA‑AAAAAAAAACCGCCAGTAAACcggc < 1:125469/147‑1 (MQ=255) ttattttTTGTTTTTCTGGCCTCTTTTTCCCGCCAAAAGGGGGATGAATAATGAAC‑GAAATAACTCCCTCCCCTACCCGCCCCTTTGTA‑ACAT‑TGGTGAAAAGGGCGCT‑TTTTTTTTATGGCGGC‑AAAAAAAAACCGCCAGTAAACcggc < 2:245383/150‑1 (MQ=255) tattttTTGTTTTTCTGGCCGCGTTTTTCACCCAAAAAGGGAATGACAAATGAACGGAGAAAA‑CCCCCCCCCTACCCGCCCCTTTGTT‑ACAG‑TTGTGTAAAGGGGGCT‑TGTTTTTTATGAGGGC‑AAAAAAAAACCGCCAGTAAACcggcg < 1:318164/150‑1 (MQ=255) tattttTGGTTTTTCGGCCCGCTTTTTCCACCCAAAAAGGGGATGACAAAGGAA‑GGAAAAAA‑ACCCCCCCCCAACCCGCCCCTTTTTAACAT‑TTGGGAACAGGGGGCT‑TGTTTTTTAGGACGGA‑AAAAAAAAACCGCCAGTAAACcggcg < 2:229468/150‑1 (MQ=255) ttttctgtttttctggCCGTTTATTCCACCCAAAAATGGGATGACAAATGAACGGAAATAA‑CCCCCCCCCTACCCGCCCCTTTTTT‑ACAG‑TTGTGTAAAAGGGCCT‑TGTTTTTTATGACGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGt < 1:458701/150‑1 (MQ=255) ttttgttttgctggCCGTTTATTCCACCCAAAAATGGGTTAACAAAGGAAGGGAAATAT‑CCCCCCACCAACCCGGCCCCTTTTA‑AAAG‑TTGGGTCAAAGGGGCC‑TGTTGTTTA‑GCGGGCAAAAAAAAAACCGCCAGTAAACCGGCGGTGa < 2:242058/144‑1 (MQ=255) ttttgtggCCGCGTATTCCACCCAAAAAGGGGATGACAAATAAAGGGAAAAAA‑CCCCCCACCAACCCGGCCCTTTTTA‑ACGG‑TTGGGTAAAGGGGGCC‑TGTTTTTTATGACGGC‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt < 2:326744/150‑1 (MQ=255) tttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc > 1:154883/1‑151 (MQ=255) tctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCAt > 1:156412/1‑151 (MQ=255) tggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATgg > 1:370288/1‑151 (MQ=255) ggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg > 2:382429/1‑149 (MQ=255) ccGGGTATTCCACCAAAAATTGGGATGACAAATGACGGGAAATAA‑CCCCTC‑CCCAACCCGCCCCTTTTT‑TAAGTTTGTGTAAAAGGGGCT‑TATTTTTTATGCGGGG‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt < 1:919226/150‑1 (MQ=255) cgcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt > 1:266132/1‑149 (MQ=255) aTTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt > 2:861273/1‑150 (MQ=255) aTGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGaa > 2:20196/1‑118 (MQ=255) ggATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGc > 2:65470/1‑151 (MQ=255) aCTAATGACCGGAGATAA‑CCCCCCCCCTACCCGGCCCCTTGTT‑ACAG‑TTGTGTACAGGGGGCC‑TTTTTTTTATGACGAA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTa < 1:864538/148‑1 (MQ=255) aaggAAGGGAAATAA‑CCCCCCACC‑AACCCGCCCCTTTTT‑TACG‑TTTTGGAAAAAGGGCCTTGTTTTTTATGAGGAA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACa < 2:235655/145‑1 (MQ=255) aaTGAAGGGAAATAA‑CCCCTCCCCAACCCGGCCCTTTTTT‑ACGG‑TTGTGAAAAAGGGGCT‑TGTTTTTTAGGCCGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTa < 1:705039/145‑1 (MQ=255) aagAAGGGAAATAC‑CCCCTCCCCAACCCGCCCCTTTTTT‑ATGT‑TTGTGTAAAGGGGGCT‑TTTTTTTTAGGGGGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACa < 1:495470/145‑1 (MQ=255) aaCGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAttt > 2:709705/1‑151 (MQ=255) aCGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtttt > 1:908685/1‑151 (MQ=255) agaTAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTg > 2:8858/1‑151 (MQ=255) cccccccctccccGGCCCTTTTTT‑AATT‑TTGTGAGCAGGGGTCC‑TTTTTTTTAGGGCGAA‑AAAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc < 1:6775/134‑1 (MQ=255) ccccccccaacccGCCCCTTTTTT‑ACAG‑TTGTGTAAAGGGGGCT‑TGTTTTTTATGACGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTg < 1:803407/138‑1 (MQ=255) cccTCACCAACCCGGCCCTTTTTT‑AGGG‑TTGTGAACAGGGGGCT‑TTTTTTTTATGACGGA‑AAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt < 1:675811/151‑1 (MQ=255) ccTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt > 2:617673/1‑150 (MQ=255) aCCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TTATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACCCTAACAGGCATTTTCCTGCACTGATAACg > 2:961599/1‑150 (MQ=255) ccTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa > 1:892502/1‑150 (MQ=255) acccGCCCCTTTTTT‑ACAG‑TTGTGAAAAGGGGGCT‑TGTTTTTTAGGACGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATc < 2:154871/148‑1 (MQ=255) gcccctTTTTT‑ACGG‑TGGTGAAAAGGGGGCT‑TTTTTTTTAGGGGGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATc < 2:27798/144‑1 (MQ=255) t‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGc > 1:957025/1‑150 (MQ=255) gtgcACAGGGGCCC‑TGATTTTTGGGACGGA‑AAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAg < 1:202373/148‑1 (MQ=255) cAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc > 1:934447/1‑150 (MQ=255) aaagggggctttttttttAGGACGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc < 1:91381/138‑1 (MQ=255) aaGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa > 1:305504/1‑150 (MQ=255) aaGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa > 1:373259/1‑150 (MQ=255) ggggcCC‑TGTTTTTTAGGACGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaa < 1:751209/149‑1 (MQ=255) ggggctttttttttAGGACGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa < 1:657587/139‑1 (MQ=255) gggtgtgtttttaAGGACGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga < 1:495973/143‑1 (MQ=255) ggcctttttttttATGAGGAA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga < 2:914159/143‑1 (MQ=255) c‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGt > 1:327141/1‑151 (MQ=255) aTGACGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGa < 1:91897/150‑1 (MQ=255) cGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCtt > 2:768127/1‑150 (MQ=255) | TTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAA‑TCCCTCACCTAACCGGCCCCTTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTT > NC_000913/1622667‑1622961 |
| 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: ‑ |