Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,622,816 | 0 | G | A | 52.2% | 6.3 / 32.1 | 23 | intergenic (+75/+144) | yneM/dgcZ | inner membrane‑associated protein/diguanylate cyclase, zinc‑sensing |
Reads supporting (aligned to +/- strand): ref base G (10/1); new base A (0/12); total (10/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.61e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-04 |
GGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACT > NC_000913/1622681‑1622966 | 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) g‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCCTGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACt > 2:74302/1‑151 (MQ=255) | GGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACA‑GTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACT > NC_000913/1622681‑1622966 |
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: ‑ |