| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 298,960 | 0 | A | C | 42.9% | 88.0 / 19.4 | 56 | Y659D (TAT→GAT) | paoC | PaoABC aldehyde oxidoreductase, Moco‑containing subunit |
| Reads supporting (aligned to +/- strand): ref base A (26/6); new base C (0/24); total (26/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.54e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACAC > NC_000913/298811‑299107 | gTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATAc < 2:920412/151‑1 (MQ=255)gTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATAc < 2:1017665/151‑1 (MQ=255) acacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACTCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑GGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGc < 1:867707/151‑1 (MQ=255) gcgTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCatat > 2:1171051/1‑151 (MQ=255) gTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCg > 2:533001/1‑151 (MQ=255) tGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTg > 1:52585/1‑151 (MQ=255) aTAGCCGCGCTCCCGCCGCACAGGCCCAGCTCCCCGCCCCCTTTGGCCTTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGGAAAACCCCCTCC‑TGTTTGGG‑G‑ATATCCGCATGAACCGCCACCTCATCCCCCGCCATATCGTGATTa < 1:917050/150‑1 (MQ=255) tggCCGCGCTCCCGCCGCACAGCCCCACCCCACCGACCCTTTTGGCCTTCATGGGGGAGGATATGGGGTGGGTATCATCCAGAAAAACCCCCTCT‑TTTTTTGG‑G‑ATATCCGCATAACCCGGCACCTCATACCCCGCCATATCGTGAt < 2:1088387/145‑1 (MQ=255) gCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACg > 1:1076250/1‑150 (MQ=255) cTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGc > 2:667520/1‑151 (MQ=255) cACGCCGCACAGGCCCAGCTCCCCGACCCTTTTGGCTTTCATCGGGGAGGATATGGGGTGGTTATCATCCAGAAAAATCCCCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGccc < 1:641256/151‑1 (MQ=255) cgccgcACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCa > 1:379916/1‑150 (MQ=255) gccgcACGGGCCCAGCTCCCCGACCCTTTTGGCCTTCATGGGGGAGGATATGGGGTCGGTATCATCCAGGAAAACCCCCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCa < 1:50245/149‑1 (MQ=255) cgcACAGGCCCAGCTCACCGACCCTTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGAAAAATCCCCTCC‑TGTTTTGG‑G‑ATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACg < 1:162289/151‑1 (MQ=255) acaGGCCCAGCTCACCGACCCCTTGGCCCTTCATGGGGGAGGATATGGGGCCGGTATCACCCAGAAAAACCCCCCCC‑TGTTTTGG‑G‑ATATCCGCAGGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACgg < 2:349225/149‑1 (MQ=255) acaGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc > 1:820671/1‑151 (MQ=255) cccAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATcc > 2:838873/1‑150 (MQ=255) ccAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCAc > 1:537369/1‑151 (MQ=255) gACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc > 2:206235/1‑151 (MQ=255) cacCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcctc > 2:434956/1‑151 (MQ=255) ctttgggCTTTCAGGGGGGAGGATATGGGGTCGGTATCACCCAAAAAAACCCCCTCC‑TGTTTGGG‑G‑ATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAg < 2:168585/138‑1 (MQ=255) tttGCCCTTCGGCGGGGAGAATAGGGGGTCGTTATCACCCAGAAAAACCCCCCCC‑TTTTTGGG‑G‑ATACCCCCATGAACCGCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcctc < 2:655307/147‑1 (MQ=255) ggCTTTCAGGGGGGAGAATAGGGGGTCGGTATCACCCAGAAAAATCCCCTCC‑TGTTTTGG‑A‑ATACCCGCAGGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAg < 1:447825/150‑1 (MQ=255) gCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTg > 1:475256/1‑151 (MQ=255) ccTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc > 1:1127938/1‑151 (MQ=255) ccTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc > 1:475613/1‑151 (MQ=255) ttCATCGGGGAGGATATGGGGTCGTTATCATCCAGAAAAACCCCCTCC‑TTTTTTGG‑A‑ATACCCGCATGACCCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc < 1:874311/149‑1 (MQ=255) tccATGGGGGAGGTTGGGGGGTGGTTACCA‑CCAGGAAAAACCCCCCCTTGTTGGGG‑G‑ATACCCGCATGACCCGGCCCCCCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc < 2:263240/147‑1 (MQ=255) cacgggggAGGATAGGGGGTCGTTATCACCCA‑GAAAAACCCCCCC‑TTTTTTGG‑GAATCCCCGCATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcg < 1:862537/145‑1 (MQ=255) aTCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCgcgc > 1:408952/1‑151 (MQ=255) tagggggTCGGTACCCCCAAGAAAACTCCCCCCC‑TTTTTTGG‑G‑ATCCCCGCAGGACCCGGCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgccca < 2:1049646/146‑1 (MQ=255) gggTCGGTATCACCCAGAAAAATCCCCTCC‑TGTTTGGG‑A‑ATACCCGCAGAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCa < 1:637137/151‑1 (MQ=255) cGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgc > 2:451063/1‑151 (MQ=255) tATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGcc > 2:856849/1‑151 (MQ=255) atcCAGAAAAACCCCCTCC‑TTTTTTGG‑A‑ATACCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCagggagaaaat < 1:635412/150‑11 (MQ=255) ccAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAc > 2:203706/1‑151 (MQ=255) ccAGAAAAATCCCCTCT‑TTTTTTGG‑A‑ATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGa < 1:743859/150‑1 (MQ=255) aaaTCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCt > 1:763798/1‑151 (MQ=255) cccctcc‑tTTTTGGG‑G‑ATACCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgc < 1:926226/149‑1 (MQ=255) ccccgtttttGGAG‑ATACCCCAATGACCCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTg < 2:177527/136‑1 (MQ=255) cc‑tGTTTTGG‑G‑ATTTCCGCATGAACCGGCACCCAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGt < 1:1139136/151‑1 (MQ=255) tGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt > 2:571663/1‑151 (MQ=255) tGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt > 1:267296/1‑151 (MQ=255) tGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt > 2:531245/1‑151 (MQ=255) ttttGG‑G‑ATATCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt < 1:497568/150‑1 (MQ=255) tttGG‑G‑ATACCCGCATAACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTc < 1:764850/150‑1 (MQ=255) ttggg‑g‑ATACCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTc < 1:539871/147‑1 (MQ=255) ccGCAGGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTc < 1:214820/146‑1 (MQ=255) cGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATg > 2:3157/1‑151 (MQ=255) acccgcccccTCAACCCCCGCCATGTCGTGATTAACGAAGTGCCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTc < 2:13050/141‑1 (MQ=255) acccGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTc < 1:310905/147‑1 (MQ=255) gcccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgca < 1:564573/146‑1 (MQ=255) cACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgcagc > 2:415763/1‑150 (MQ=255) cTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacac > 1:862759/1‑151 (MQ=255) cATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacacac < 2:565961/151‑1 (MQ=255) tACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGAGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacacac > 2:102194/1‑149 (MQ=255) | GTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCC‑TGTTTTGG‑G‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACAC > NC_000913/298811‑299107 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |