| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,004,625 | 0 | T | G | 53.3% | 36.8 / 45.0 | 30 | L227V (TTG→GTG) | ycbF | putative periplasmic pilin chaperone |
| Reads supporting (aligned to +/- strand): ref base T (3/11); new base G (16/0); total (19/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.66e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-05 | |||||||||||
ACATTAACCTGATTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGTTGAGTGCAAAGAAGGAGCAAAATCTG‑CCCTGAAACAGGTTC‑GGAAAACGTTTGCGTTTTTTTTGCCGCAGGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTAC > NC_000913/1004479‑1004773 | aCATTAACCTGATTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGTTggg > 1:79811/1‑149 (MQ=255)aCATTAACCTGATTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGGTgtg > 2:266858/1‑149 (MQ=255) aCCTGATTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGggtgggtgggaa > 2:352247/1‑147 (MQ=255) aCCTGATTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGGTGGGGGCaaa > 1:600349/1‑151 (MQ=255) cTGATTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGTTGAGTGCAaaga < 1:24611/151‑1 (MQ=255) aTTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGGTGTGGGGAAAGAAgg > 2:471882/1‑151 (MQ=255) aTTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGGTGGGGGGGAAGAAgg > 1:271650/1‑151 (MQ=255) gTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGTTGAGTGCAAAGAAGGAGCa < 1:352252/149‑1 (MQ=255) gtgtAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGGTGAGTGCAAAAAAAGAGGaaaaa > 1:88277/1‑148 (MQ=255) gtgtAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGGTGAGGGGAGAAAAGGGGGaaaaaat > 2:545407/1‑148 (MQ=255) aTGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGGTGGGTGGAAAAAAAGAGGaaaaaattcc > 1:426129/1‑142 (MQ=255) ggTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTTGGGGGCAAAAAAGAGGAAAAAACGG‑cccgaa > 1:368622/1‑148 (MQ=255) gCGTTAAGTAATGTGGGAGTTGAGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGGTGAGTGGAAAAAAGGAGGAAAAACTG‑CCCCGAAAAAggg > 1:309031/1‑150 (MQ=255) tAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGTTGAGTGCAAAGAAGGAGCAAAATCTG‑CCCTGAAACAGGTTC‑GGaa < 2:208311/149‑1 (MQ=255) tgtgGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGGTGAGTGCGAAGAAAGGGCGAAAACTG‑CCCTGTAACAAGGTC‑gggacaagg > 2:144815/1‑144 (MQ=255) ccTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGGTGGGGGGGAAGAAGGAGCAAAAACTG‑CCCCGAAAAAAGGTC‑GGGAAAAGGTTTTGtttttttttt > 1:215851/1‑148 (MQ=255) gCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGGTGGGGGGAAAGAAAGAGGAAAAAATGACCCGAAAAAAGGGTG‑GGAAAAGATTGTGGtttttttttcgccggggaaat > 2:264119/1‑138 (MQ=255) aGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGTTGGGGGCAAAAAAAGAGCAAAAACTG‑CCCTGAAAAAGGGTC‑GGGAAAAGATTGCGTTTTTTTTTTCGCAGGGC‑AAATTCCtt > 1:605048/1‑151 (MQ=255) aGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGGTGAGGGGGAAAAAGGAGCGAAAGCAG‑CCCCTGAACAGGCGCGGGAAAAGGTTTAGGTTTTTTTTTCCGGGGATA‑ATTTCCTTTTTgt > 2:98311/1‑150 (MQ=255) cTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGTTGAGTGCAAAGAAGGAGCAAAATCTG‑CCCTGAAACAGGTTC‑GGAAAACGTTTGCGTTTTTTTTGCCGCAGGTC‑AATTCCCTTTTGGTc < 1:124845/150‑1 (MQ=255) tCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGTTGAGTGCAAAGAAGGAGCAAAACCTG‑CCCTGAAACAGGTTC‑GGAAAACGTTTGCGTTTTTTTTGCCGCAGGTC‑AATTCCCTTTTGGTCCGAACTCGCACataa < 2:215848/151‑1 (MQ=255) cGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGTTGAGTGCAAAGAAGGAGCAAAATCTG‑CCCTGAAACAGGTTC‑GGAAAACGTTTGCGTTTTTTTTGCCGCAGGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACGcc < 1:476069/151‑1 (MQ=255) gtgtCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGGTTGGTGGAAAAAAGGAGGAAAAAATT‑CCCTGAAAAAAGGTC‑GGGAAAAGGTTTCGGTTTTTTTTTCGCCGGGC‑AAATTCCTTTTTGTTCGAACTCGCCCAAAATAAGcc > 2:158488/1‑150 (MQ=255) gtgtCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGGTGAGGGGCAAGAAGGAGCAAAATCTG‑CCCCGAAAAAAGGTC‑GGGAAACGGTTGCGGTTTTTTTTTCGCAGGGC‑AATTTCCTTTTTGTCCGAACTCGCCCATAATACAcc > 1:601286/1‑150 (MQ=255) gtgtCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGGTGAGGGGAAAGAAGGGGCAAAAACTG‑CCCTGAAAAAGGGTC‑GGGAAAAGGTTGCGGTTTTTTTTTCGCCGGGCAAATTCCTTTTTGGT‑CGAACTCGCACataatatacc > 2:211542/1‑146 (MQ=255) tGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGTTGAGTGCAAAGAAGGAGCAAAATCTG‑CCCTGAAACAGGTTC‑GGAAAACGTTTGCGTTTTTTTTGCCGCAGGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGtt < 1:450936/151‑1 (MQ=255) cGCCTGGCGGGCGATAAATGATTACGGC‑GGGGTGGGTGCAAAGAAGGAGGAAAAACTG‑CCCCGAAAAAAGGTC‑GGAAAAAGTTTGCGGTTTTTTTTCCCCAGGGC‑AAATCCCTTTTTGTCCGGACTCGCACATAATACGCCCCCCGGTTg > 1:131660/1‑150 (MQ=255) tGATTACGGC‑GGGTTGAGTGCAAAGAAGGAGCAAAATCTG‑CCCTGAAACAGGTTC‑GGAAAACGTTTGCGTTTTTTTTGCCGCAGGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGgag < 2:250572/151‑1 (MQ=255) tGATTACGGC‑GGGTTGAGTGCAAAGAAGGAGCAAAATCTG‑CCCTGAAACAGGTTC‑GGAAAACGTTTGCGTTTTTTTTGCCGCAGGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGgag < 2:426120/151‑1 (MQ=255) gTTGAGTGCAAAGAAGGAGCAAAATCTG‑CCCTGAAACAGGTTC‑GGAAAACGTTTGCGTTTTTTTTGCCGCAGGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGtac < 1:512782/150‑1 (MQ=255) gTTGAGTGCAAAGAAGGAGCAAAATCTG‑CCCTGAAACAGGTTC‑GGAAAACGTTTGCGTTTTTTTTGCCGCAGGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGtac < 2:445806/150‑1 (MQ=255) | ACATTAACCTGATTCAGGTGAGTGTAAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGC‑GGGTTGAGTGCAAAGAAGGAGCAAAATCTG‑CCCTGAAACAGGTTC‑GGAAAACGTTTGCGTTTTTTTTGCCGCAGGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTAC > 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 < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |