| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,622,815 | 0 | C | A | 50.0% | 10.8 / 15.2 | 22 | intergenic (+74/+145) | yneM/dgcZ | inner membrane‑associated protein/diguanylate cyclase, zinc‑sensing |
| Reads supporting (aligned to +/- strand): ref base C (9/2); new base A (0/11); total (9/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.21e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAAC > NC_000913/1622681‑1622965 | gggAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGca < 1:214621/149‑2 (MQ=255) tattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACcggcgg > 1:165885/1‑151 (MQ=255) ttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGa > 2:337419/1‑150 (MQ=255) ccggattttccacccAAAAAGGGAATGAATAAGGAAGGGAGATACCCCCCCACCAACCCGCCCCTTTGTTAAT‑TTTGTGAAAAGGGGTCC‑TTTTTTTTAGGAGGACA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa < 1:534302/132‑1 (MQ=255) tATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt > 2:86700/1‑151 (MQ=255) aaaagggggATGACTAATGAAGGGAGATACCCCCCCCCCTACCCGCCCCCTTTTTACA‑GTTGTGAAAAGGGGCCC‑TTTTTTTTATGCGGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTg < 1:381340/143‑1 (MQ=255) aaaaGGGGATGATAAATGAAGGGAAATAACCCCCCCCCTACCCGCCCCTTTTTTACG‑TTTTTGTACAGGGGGCT‑TGTTTTTTAGGCCGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtgtg < 1:219053/142‑1 (MQ=255) aCTAAGAAAGGGAAAAAACCCCCCACCAACCCGCCCCTTTTTTACA‑TTGGTGAAAAGGGGGTG‑TGTTTTTAAGGCGGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCt < 1:578656/147‑1 (MQ=255) cGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTc > 2:242094/1‑151 (MQ=255) agaTAACCCCTCACCAACCCGCCCCTTTTTTACA‑GTTGTGTAAAGGGGGCG‑TGTTTTTTAGGAGGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTc < 2:80149/148‑1 (MQ=255) cccTCCCCTACCCGCCCCTTTTTTACA‑TTTGGTAACAGGGGGCT‑TGTTTTTTAGGACGAAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt < 2:203795/151‑1 (MQ=255) cccaaCCCGCCCCTTTTTAAATGTTGTG‑AAAAGGGGCCTTTTTTTTTAGGAGGGAA‑AAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATa < 2:133067/144‑1 (MQ=255) acccGCCCCTTTTTTACA‑TTTGTAAACAGGGGCCT‑TTTTTTTTATGACGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa < 1:220763/145‑1 (MQ=255) cgcccctTTGTTACT‑GTGGTGAACAAGGGGCT‑GTTTTTTTATGAGGAAA‑AAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCg < 2:225149/145‑1 (MQ=255) ggCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGtt > 2:375815/1‑149 (MQ=255) cA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc > 2:191625/1‑151 (MQ=255) tgtgTACAGGGGCCT‑TGTTTTTTAGGACGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGc < 1:470157/144‑1 (MQ=255) ggggCT‑TGTTTTTTAGGACGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa < 1:406842/148‑1 (MQ=255) cttgttttttAGGAGGGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaat < 1:604003/145‑1 (MQ=255) tGATTTTTATGACGGCGAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGt > 2:369444/1‑151 (MQ=255) ttttATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaaa > 2:143351/1‑150 (MQ=255) gACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCt > 1:557367/1‑151 (MQ=255) cG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAAc > 2:626405/1‑151 (MQ=255) | GGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAAC > NC_000913/1622681‑1622965 |
| 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: ‑ |