| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 298,960 | 0 | A | C | 50.9% | 65.5 / 101.1 | 55 | Y659D (TAT→GAT) | paoC | PaoABC aldehyde oxidoreductase, Moco‑containing subunit |
| Reads supporting (aligned to +/- strand): ref base A (24/3); new base C (0/28); total (24/31) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.18e-12 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-09 | |||||||||||
GTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCC‑TTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACAC > NC_000913/298811‑299107 | gTTATACCCCGCGTTGGCGATAGCCGCGCTCACGCCGCTCAGGCCCAGCTCACCGACACCTTTGGCC‑TTCATCGGGGATGATATGGGGTCGGTATCATCCAGGAAAAT‑CTCCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCAAC‑TCATa < 1:753751/150‑1 (MQ=255)gTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCC‑TTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATAc > 1:634303/1‑151 (MQ=255) tataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCC‑TTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATAccc > 2:812027/1‑151 (MQ=255) ataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCC‑TTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATAcccc < 1:504456/151‑1 (MQ=255) ccGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCC‑TTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCat < 1:469713/151‑1 (MQ=255) cgTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCC‑TTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATc > 2:907055/1‑151 (MQ=255) gCGATACCCGCGCCCCCGCCGCACAGGCCCAGCTCCCCGCCCCTTTTGGCC‑TTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGGAAAAT‑CCCCTCCTGTTTGGGG‑ATCTCCGCATGACCCGGCCCC‑TAATCCCCCGCCATATCGTg < 1:289014/149‑1 (MQ=255) ccGCCCTCACGCCGCACGGCCCCAGCGCCCCGACCCCTTGGGCC‑TTAATCGGGGAGGATTTGGGGTCGGTATCACCCAGGAAAAT‑CCCCTCCTGTTTTGGA‑ATACCCGCATGACCGGCCACC‑TCATCCCCCGCCATATCGTGATTAACg < 1:801842/149‑1 (MQ=255) gcgcTCACGCCGCACAGGCCCAGCTCACCGACCCCTTTGGCC‑TTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGAAAAAC‑CCCCTCCTGTTTGGGG‑ATATCCGCATGAACCGGCCCC‑TAATCCCCCGCCATATCGTGATTAACGAAGt < 1:887432/151‑1 (MQ=255) ccgaACAGGCCCAGCTCCCCGACCCCTTTGGCC‑TTCAGCGGGGGGGATAGGGGGTCGTTACCATCCAGGAAAAT‑CCCCCCCTGTTTTGGG‑ATACCCGCAGAACCCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCaaa < 2:310422/150‑1 (MQ=255) cacaGGCCCAGCTCACCGACACCTTTGGCC‑TTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGt > 1:333060/1‑151 (MQ=255) acgggCCCAGCTCCCCCACCCCTTTGTCT‑TCCACGGGGGAGGATGTGGGCTGGGTAACACCAAGTAAACA‑CGCCT‑CTGTTTTGGGAATCACCGATTGAACCGGCACC‑TCATCCCCCGCCATATCGTGATTAACGAAGt < 2:67463/135‑1 (MQ=255) ggCCCAGCTCCCCGCCACCTTTGGCT‑TTCATCGGGGGGGATAGGGGGTCGGTATCCTCCAGGAAAAT‑CCCCTCCTGTTTTGGG‑ATACCCACATGACCCGGCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAt < 1:160842/150‑1 (MQ=255) cccAGCTCACCGACACCTTTGGCC‑TTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATcc > 1:318073/1‑150 (MQ=255) aGCTCACCGTCGCTTTTGGCC‑TTGATGGGGGGGAATAGGGGGTGGGTATCGCCAAGAAAAAC‑CCCCTCCTTTTGTGGA‑AGACCCCTATGAGCCGGCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCg < 2:625646/151‑1 (MQ=255) ccccGACCCCTTTGGCC‑TTCATGGGGGGGGATAGGGGGTGGGTATCACCCAGGAAAAT‑CCCCCCCTGTTTGGGG‑ATACCCGCATGACCCGGCACC‑TAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGc < 1:852604/146‑1 (MQ=255) cACCGGCCCCTTTTGCCTTTAACCGGGGAGGATAGGGGTTGGGTACCAACAAGAAAAAC‑CCCCTCTTGTTTTGGA‑ATATCCGCATGACCCGGCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCa < 2:653953/151‑1 (MQ=255) ccccTTTGCCC‑TTCAGGGGGGAGGATTGGGGGTCGTTATCATCCAGAAAAAT‑CCCCTCCTTTTTTGGA‑ATATCCGCATGAACCGGCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcctc < 1:3658/149‑1 (MQ=255) ccccTTGGCCC‑TTCAGCGGGGAGATTGGGGGGTCGTTACCACCCAGAAAACC‑CCCCTCCTGTTTGGGG‑ATCCCCGAATGACCCGCCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc < 2:191915/147‑1 (MQ=255) tttGGCC‑TTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc > 1:590210/1‑151 (MQ=255) ttGGCC‑TTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc > 1:324729/1‑150 (MQ=255) ttgcctttcaggggggAGGATAGGGGGTGGTTACCACCCAGAAAAAC‑CCCCCCCTTTTTTGGG‑ATACCCGCATAACCCGGCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc < 1:226347/138‑1 (MQ=255) aTCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCgcgc > 1:157851/1‑151 (MQ=255) ggaggaTATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatg > 1:150355/1‑151 (MQ=255) gatggggggcgggaatcaccCAGGAAAAC‑CCCCTCCTGTTTGGGG‑ATACCCCCATGACCCGCCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgccc < 1:211039/136‑1 (MQ=255) atGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat > 1:895104/1‑149 (MQ=255) tGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGt > 2:887221/1‑151 (MQ=255) tGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGt > 2:107621/1‑151 (MQ=255) ggggTGGTTACCACCAAGGAAAAC‑CCCCTCTTGTTTGGGG‑ATCCCCGCAGAACCCGCCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATa < 2:911901/148‑1 (MQ=255) ggTATCACCCAGAAAACC‑CCCCTCCTTTTTGGGG‑ATACCCGAATGACCCGCCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTg < 1:920070/149‑1 (MQ=255) ggTACCACCCA‑GAAAAA‑CCCCCCCTTTTTTGGGAATACCGCCAGGACCCGCCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgc < 2:934446/150‑1 (MQ=255) gTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcgc > 2:76662/1‑151 (MQ=255) catcCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAAt > 1:31271/1‑151 (MQ=255) cacccagAAAAA‑CCCCCCCTTTTTTGGGAATCCCCGCATGACCCGCCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAAt < 2:902944/145‑1 (MQ=255) acccAGGAAAAC‑CCCCTCCTGTTGGGGG‑ATATCCGCATGACCCGCCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg < 1:849547/149‑1 (MQ=255) acccggAAAATCCCCCTCTTTTTTTGGA‑ATACCCACATG‑ACCGGCCCCTTCTTCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcg < 2:586439/141‑1 (MQ=255) ccAGAAAAAC‑CCCCTCCTTTTTGGGA‑ATACCCGCATGACCCGCCCCC‑TCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg < 2:49047/149‑1 (MQ=255) ccAGAAAAAC‑CCCCTCCTGTTTTGGG‑ATATCCGCATGACCCGCCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAc < 1:34399/151‑1 (MQ=255) aGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt > 1:207193/1‑151 (MQ=255) aaaaaC‑CCCCTCTTGTTTGGGT‑ATACCCCAATAACCCGCCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAcc < 2:708062/147‑1 (MQ=255) ccccctTTTTGGGA‑ATATCCGCATGACCCGGCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc < 1:491299/148‑1 (MQ=255) ctcctGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc > 2:620177/1‑150 (MQ=255) tGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt > 1:677738/1‑151 (MQ=255) ttttGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTc > 1:241496/1‑151 (MQ=255) ttttGGA‑ATACCCGCATGCCCCGGCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt < 2:967559/150‑1 (MQ=255) tGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg > 1:337986/1‑150 (MQ=255) ggaatacCCGCATGACCCGCCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg < 2:265454/146‑1 (MQ=255) tatCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTc > 1:696040/1‑149 (MQ=255) ccGCAGGACCCGCCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGa < 1:746079/150‑1 (MQ=255) aTGACCCGGCCCC‑TCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc < 1:406851/149‑1 (MQ=255) aTGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGt > 2:365869/1‑151 (MQ=255) gAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTcc > 1:485183/1‑151 (MQ=255) cggccccc‑TCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc < 2:200235/145‑1 (MQ=255) ccGGCCCC‑TCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTg < 2:684764/150‑1 (MQ=255) ccGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc > 1:531854/1‑151 (MQ=255) cATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacacac > 2:964535/1‑151 (MQ=255) | GTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCC‑TTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACC‑TCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACAC > NC_000913/298811‑299107 |
| 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: ‑ |