| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,004,625 | 0 | T | G | 52.0% | 20.0 / 29.4 | 25 | L227V (TTG→GTG) | ycbF | putative periplasmic pilin chaperone |
| Reads supporting (aligned to +/- strand): ref base T (4/8); new base G (13/0); total (17/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.58e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-03 | |||||||||||
ACATTAACCTGATTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGA‑GCAAAATCTG‑CCCTGAAACA‑GGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAG‑GTC‑AATTCCCTTTTGGTCCGAACTCGC‑ACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTAC > NC_000913/1004479‑1004773 | aCATTAACCTGATTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAg > 2:276259/1‑151 (MQ=255)aCATTAACCTGATTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTggg > 1:436852/1‑149 (MQ=255) aCCTGATTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGAGGGGaaa > 2:24105/1‑151 (MQ=255) aTTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGGGTGCAAAGAAgg > 1:197578/1‑151 (MQ=255) aTTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGTGTGGAAAGAAgg > 1:13331/1‑151 (MQ=255) aTTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGAGGGGAAAGAAgg > 2:94304/1‑151 (MQ=255) ggTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGGGGGGCAAGAAaga‑ggga > 2:315637/1‑148 (MQ=255) gtgtAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGGGTGGAAAAAAAGA‑AGaaaaaat > 2:113628/1‑148 (MQ=255) gtgtAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGAGTGGAAAGAAGGA‑gagaaaaat > 1:19373/1‑143 (MQ=255) gtgtAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGAGGGCGAAGAAGGA‑GGaaaaa > 2:336288/1‑148 (MQ=255) tgtAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGA‑GCAAAATc < 2:201749/149‑1 (MQ=255) gTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGAGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGGGCGAAGAAGGAGGGAAAACCGC‑CCCGGAAAAA‑GGTTGGGaa > 1:158961/1‑151 (MQ=255) gACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGTGGGCGAAGAAAGA‑GAAAAAAATG‑CCCCCAAAAACGGTCGGGAAAAGGTTTGGGTTTTTTTTTCCCCCG‑GCC‑AATCCCCtt > 1:418702/1‑151 (MQ=255) tCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGA‑GCAAAATCTG‑CCCTGAAACA‑GGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAG‑GTC‑AATTCCCTTTTGGTCCGAACt < 2:11733/150‑1 (MQ=255) gCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGA‑GCAAAATCTG‑CCCTGAAACA‑GGTTCGGAAAACGTTTGCG‑TTTTTTTGCCGCAG‑GTC‑AATTCCCTTTTGGTCCGAACTCGc‑ac < 1:303867/150‑1 (MQ=255) gtgtCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGGGTGGAAAGAAGGG‑GCAAAAACTG‑CCCCGAAAAA‑AGGTCGGAAAAAGTTTTCGTTTTTTTTTTCGCCG‑GGC‑AAATTCCTTTTTGTCCGGACTCGC‑CCCTAAAAAGcc > 2:185588/1‑150 (MQ=255) gtgtCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGAGGGCAAAGAAAGA‑GGAAAAAATGCCCCGGAAAAA‑GGTCGGGAAAAGGTTT‑CGGTTTTTTTTTCGCAG‑GGCAAATTCCCTTTTGTCCCGAA‑TCGC‑ACataaaaaccc > 2:324288/1‑144 (MQ=255) tgtCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGA‑GCAAAATCTG‑CCCTGAAACA‑GGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAG‑GTC‑AATTCCCTTTTGGTCCGAACTCGC‑ACATAATACGcccc < 1:107898/151‑1 (MQ=255) cGCCTGGCGGGCGATAAATGATTACGGCGGGGTGGGGGGAAAGAAGGA‑GCAAAAACTT‑CCCTGAAACA‑AGGTCGGGAAAAGTTTGTGGTTTTTTTTTCGCCG‑GGC‑AAATCCCCTTTTGTTCGAACTCGC‑ACATAATACGCCCCCCGGTTg > 2:22804/1‑150 (MQ=255) ccTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGA‑GCAAAATCTG‑CCCTGAAACA‑GGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAG‑GTC‑AATTCCCTTTTGGTCCGAACTCGC‑ACATAATACGCCCCCGGTTTGcac < 2:345770/151‑1 (MQ=255) ggcgggcgATAAATGATTACGGCGGGGTTAGGGGAAAGAAGGG‑GCAAAAAATT‑TCCCGAAACA‑AGGTCGGCAAACGGTTTGGTTTTTTTTTTCGCCGCGGA‑AATTCCCTTTTGG‑CCGAACTCGCAACAAAAACTGCCCCCGGttgcaccc > 2:411202/1‑144 (MQ=255) ggcgATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGA‑GCAAAATCTG‑CCCTGAAACA‑GGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAG‑GTC‑AATTCCCTTTTGGTCCGAACTCGC‑ACATAATACGCCCCCGGTTTGCACACCggg < 1:37660/150‑1 (MQ=255) ggcgATAAATGATTACGGCGGGGTTAGTGGGAAAAAAGA‑GGAAAAAATG‑CCCTGAAAAA‑GGGGCGGGAAAAGGTTGTGTTTTTTTTTTCGCCG‑GGC‑AAATTCCTTTTTGTGCGCACTCGC‑ACATAATAAACCCCCCGTTTGTACACCcggg > 1:318973/1‑150 (MQ=255) tGATTACGGCGGGTTGAGTGCAAAGAAGGA‑GCAAAATCTG‑CCCTGAAACA‑GGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAG‑GTC‑AATTCCCTTTTGGTCCGAACTCGC‑ACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGgag < 1:150245/151‑1 (MQ=255) gTTGAGTGCAAAGAAGGA‑GCAAAATCTG‑CCCTGAAACA‑GGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAG‑GTC‑AATTCCCTTTTGGTCCGAACTCGC‑ACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGtac < 1:395101/150‑1 (MQ=255) | ACATTAACCTGATTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGA‑GCAAAATCTG‑CCCTGAAACA‑GGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAG‑GTC‑AATTCCCTTTTGGTCCGAACTCGC‑ACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTAC > NC_000913/1004479‑1004773 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |