| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,042,604 | 0 | A | C | 42.2% | 60.2 / 19.5 | 45 | V182G (GTG→GGG) | yqfA | hemolysin III family HylIII inner membrane protein |
| Reads supporting (aligned to +/- strand): ref base A (18/8); new base C (0/19); total (18/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.10e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCG > NC_000913/3042456‑3042748 | gcgcGTCTTATCAGGCCTACAAAGGCATACCCCTTACGCCTGCCCAATCTACAAATAGATCGCCAGACAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGGTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CAcc > 2:882297/1‑151 (MQ=255)gcgcGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CAcc < 2:554769/151‑1 (MQ=255) aTCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGaa < 1:926178/151‑1 (MQ=255) ggCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATaaa < 2:589983/151‑1 (MQ=255) ggCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATaa > 1:219931/1‑150 (MQ=255) cTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑cc > 2:1137375/1‑150 (MQ=255) aaaGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCAcgcc < 1:489250/151‑1 (MQ=255) cATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCt > 2:385590/1‑150 (MQ=255) ccATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGcc > 2:144105/1‑151 (MQ=255) aTTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGcc > 2:207963/1‑149 (MQ=255) tACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGt > 1:863882/1‑150 (MQ=255) ccggCCAAATATACAAATAGACCGCCAGAAAGTGG‑CACACCCTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAGGCGTTGGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAC‑CCACGCCGCCTACCGCCAGTAAGGt < 1:234816/148‑1 (MQ=255) ccTGCCCAATATACAAATAGATCGCCAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGTAAgg > 2:932843/1‑150 (MQ=255) ccTGCCAAATATACAAATAGATCGCCAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGTAAGGt < 2:793248/151‑1 (MQ=255) ccaaaTATACAAATAGATCGCCAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGTAAGGTAACg < 2:495785/148‑1 (MQ=255) aCAAATAGATCGCCAGAAAGTGG‑CACCCCCTCCCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CCCCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGTAAGGTAACGCTGCccg < 2:765152/150‑1 (MQ=255) aaaTAGACCCCAAG‑AAGTGGCCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CCCCCCGAGCGAATAAA‑CCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCccgcc < 1:651759/150‑1 (MQ=255) cGCCAGAAAGGGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGAAAGGTAACGCTGCCCGCCGCGAGCtt > 1:576707/1‑150 (MQ=255) ccAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAAc > 2:197418/1‑151 (MQ=255) cAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACt > 2:138739/1‑151 (MQ=255) gggccacaccCTCCCGCCGGCCACGAAGCGGTGCAAGAGGGCAGGTTTTTATGAACTGCTTTTGCAGACGAAAAAAAC‑CCCCCCGGGAGAAAAAC‑CCACGCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCAtt < 2:14771/147‑1 (MQ=255) cactcTCCCGCGGAGCAGGAACCGGTCCCAGAGGGCTTGTTTTTAGGGAGTGCTTTGGGAGCCGAAAAAAAC‑CCCCCCGCGCAAAAAAC‑CCACGCCCCCTACCGCCAGTAGGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCAt < 1:1245450/149‑1 (MQ=255) cccccccgccgAGCAGGAACCCGTGCCAGAGGGCTTGGTTGTAGGAAATGGTTTTGCAGACGAAAAAAAC‑CCCCCCGGGCAAAAAAC‑CCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATa < 1:119509/143‑1 (MQ=255) ccccccccgAGCACGAGGCCGTGCCAGAGGGCAGGTTTGTATGAAAGGCTTTTGCAGCCGAAAAAAAC‑CCCCCCGGGCAAAAAAC‑CCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATa < 1:1250028/140‑1 (MQ=255) ccgccgAGCACGAAGCCGTGCAAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAAAAAAT‑CCCCCCGAGCGAAAAAC‑CCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAAtt < 2:921673/147‑1 (MQ=255) cccccgAGCAGGAGCCCGTCCCAGAGGGCATGTTTTTATGAAATGCGTTTGCAGAGGTAAAAAAC‑CCCCCCGAGCAAATAAC‑CCCCGCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTac < 1:865294/146‑1 (MQ=255) cgccgAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTaccac > 1:1145210/1‑151 (MQ=255) gccgAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTaccacc > 1:26411/1‑151 (MQ=255) aaGCGGTGCCAGAGGGCATGTTTGTATGAAAGGCGTTGGCAGACGTAAAAAAC‑CCCCCCGGGCAAAAAAC‑CCCCCCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCg < 2:685574/145‑1 (MQ=255) aaGCGGTGCCAGAGGGAAGGTTTGTATGAAATGCGTTTGCAGGC‑TAAAAAATCCCCCCCGGGCAAAAAAC‑CCCGGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCg < 2:1167624/145‑1 (MQ=255) tGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGccc > 2:917449/1‑150 (MQ=255) aGATGGCATGGTTGTATGAAATGCTTTTGCAGACGTAAAAATT‑CCCCCCGAGCAAAAAAC‑CCCCGCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACagccagcc < 2:1064829/145‑1 (MQ=255) ggctgggtggtATGGAGTGCGTTTCAAGCCGTAGAAAAC‑CCCCCCGAGCAAATAAC‑CCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAg < 2:1215579/145‑1 (MQ=255) ggCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGcc > 1:395827/1‑148 (MQ=255) gtttGTAGGAAATGCTTTTGCAGACGTAAAAAAC‑CCCCCCGGGCAAAAAAC‑CCCCCCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAgg < 2:941044/148‑1 (MQ=255) ggTTGTAGGAAATGCGTTTGCAGACGAAAAAAAC‑CCCCCCGAGCGAAAAAC‑CCCCGCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAgg < 2:823573/150‑1 (MQ=255) ttGTATGGAACGCGTTTGCAGAGGTAAAAAAC‑CCCCCCGGGCAAAAAAA‑CCGCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTc < 2:549216/150‑1 (MQ=255) tgttATGGAAGGCGTTGGCAGACGTAAAAAAC‑CCCCCCG‑GCGAAAAAACCCCCCCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATa < 1:11877/143‑1 (MQ=255) tttgaaacgCGTTCGCAGATGTAAAAACC‑CACCCCGGGAGAAAAAA‑CCCCCCCCCCTACCGCCATTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAcc < 1:674822/142‑1 (MQ=255) gaaaTGCTTTGGCAGACGTAAAAAAC‑CCCCCCGAGCAAAAAAC‑CCCCGCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCagag < 2:1243010/149‑1 (MQ=255) gTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAAt > 1:1073426/1‑151 (MQ=255) ttGCAGAGGTAAAAAAC‑CCCCCCGGGAAAAAAAA‑CCCCCCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAAtt < 2:1155071/150‑1 (MQ=255) gTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCgg > 2:836683/1‑151 (MQ=255) gAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTgc > 1:197775/1‑151 (MQ=255) aaaaC‑CCCCCGGAGCAAAAAAC‑CCCCCCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTgcg < 1:507141/151‑1 (MQ=255) | GCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGG‑CACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAAT‑CACCCCGAGCGAATAAA‑CCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCG > NC_000913/3042456‑3042748 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |