| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 298,960 | 0 | A | C | 54.4% | 46.7 / 71.6 | 47 | Y659D (TAT→GAT) | paoC | PaoABC aldehyde oxidoreductase, Moco‑containing subunit |
| Reads supporting (aligned to +/- strand): ref base A (20/1); new base C (0/25); total (20/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.74e-12 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.81e-07 | |||||||||||
CGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTG > NC_000913/298810‑299097 | cGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATa < 1:279125/151‑1 (MQ=255) gTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATAc > 2:918014/1‑151 (MQ=255) acacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGc > 1:898327/1‑151 (MQ=255) acCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGcc > 2:581830/1‑150 (MQ=255) ttGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCg > 2:838451/1‑150 (MQ=255) aTAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTa > 1:823174/1‑150 (MQ=255) aGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACg > 2:268234/1‑151 (MQ=255) gcTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGt > 1:323657/1‑149 (MQ=255) tCACGCCGCACAGGCCCAGCTCACCGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGAAAAATCCCCTCCTTTTTTGGG‑ATATCCGCATGAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAg < 1:318914/149‑1 (MQ=255) tCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGcc > 2:698654/1‑151 (MQ=255) ccgcACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc > 1:727263/1‑151 (MQ=255) cgcACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGTTATCATCCAGGAAAACCCCCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACg < 1:789912/151‑1 (MQ=255) gcACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACgg > 2:862306/1‑151 (MQ=255) ggCCCAGCTCACCGACCCTTTTGGCCTTCATCGGGGAGGATAGGGGGTCGGTATCACCCAGGAAAATCCCCTCCTTTTTTGGG‑ATATCCGCATGAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCa < 1:598340/149‑1 (MQ=255) ggCCCAGCTCACCGACCCCTTTGGCCTTCATGGGGGAGGATAGGGGGTCGGTATCACCCAGGAAAATCCCCCCCTTTTTTGGG‑ATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATc < 1:383921/151‑1 (MQ=255) acccGACCCCTTTGCCTTTCAGCGGGGAGTATATGGGGCCGCTACCACCCAGGAAATTC‑CCCCCTTTTTTGGGTATACCCGCAAAACCCCGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCa < 2:795566/148‑1 (MQ=255) ccccTTTGCCCTTCATCGGGGAGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTGTTTGGGG‑ATTTCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc < 1:652396/147‑1 (MQ=255) cacCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcctc > 2:457371/1‑151 (MQ=255) tttGGCCTTCATGGGGGAGGATATGGGGTCGGTATCACCCAGAAAAACCCCCTCCTGTTTGGGG‑ATACCCGCATGACCCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCa < 2:755696/149‑1 (MQ=255) gcccTTCATGGGGGAGGATATGGGGTCGGTATCACCCAGAAAAATCCCCCCCTGTTTGGGG‑ATACCCGCATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCa < 2:531335/147‑1 (MQ=255) gCCTTCATGGGGGAGGATATGGGGTCGGTATCATCCAGAAAAACCCCCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAg < 1:581898/149‑1 (MQ=255) ctttcacgggggAGGATAGGGGGTCGGTATAACCCAGGAAAA‑CCCCCCCTTTTTTGGGAAAATCCGAATGAACCGGCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCa < 2:157748/139‑1 (MQ=255) ccTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc > 2:616479/1‑151 (MQ=255) tCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgc > 1:911014/1‑151 (MQ=255) ggggAGGGTTGGGGGCCGTTACCCTCCGGAAAAACCACCGCTTGTTTGGGG‑ATTGCCGGGTGATCCGCCTCCCCAGTCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGcc < 2:476411/144‑1 (MQ=255) ggaggaTATGGGGTTGGTATAATCCAGGAAAATCCCCTCCTTTTTGGGG‑ATATCCCCATAAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatg < 1:536475/151‑1 (MQ=255) ggggTCGGTATCACCCAGAAAAACCCCCTCCTTTTTTGGA‑ATATCCGCAGGAACGGCCCCCCAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcc < 1:960289/144‑1 (MQ=255) gggTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc > 2:37242/1‑150 (MQ=255) gTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCAtt > 1:959912/1‑151 (MQ=255) tATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGcc > 2:833522/1‑151 (MQ=255) cccAGGAAAACCCCCCCCTTTTTGGGG‑ATACCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCaa < 2:829115/147‑1 (MQ=255) cgagaaaaaTCCCCTCTTTTTTTGGG‑AAACCCGCATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGa < 2:689985/145‑1 (MQ=255) ccAGAAAAACCCCCTCCTTTTTGGGG‑ATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg < 1:189492/149‑1 (MQ=255) ccAGAAAAACCCCCTCCTGTTTTGGG‑ATACCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg < 1:850752/149‑1 (MQ=255) ggAAAATCCCCTCCTTTTTGGGA‑ATACCCGAATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAc < 2:175578/148‑1 (MQ=255) ggAAAACCCCCCCCTGTTTTGGG‑ATATCCGTTTAACCGGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt < 1:485367/150‑1 (MQ=255) aaacccccccctGTTTGGGG‑ATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCt < 1:578080/142‑1 (MQ=255) aCCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcg > 2:270825/1‑151 (MQ=255) cctcctGTTTTGGG‑ATATCCGAATGAACCGGCACCTCATACCCCGCCTTATCGTGATTAACGAAGCAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATCCCCATAGTCATTGCGCCACTGCCCTGGCTgcgcgc > 2:151046/1‑151 (MQ=255) ttttGGG‑ATTCCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAg < 2:485243/146‑1 (MQ=255) tttgGGG‑ATACCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt < 2:656140/149‑1 (MQ=255) ggg‑ATATCCGCATGACCCGCCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg < 1:15518/149‑1 (MQ=255) tacccGCAGGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCa < 2:685769/147‑1 (MQ=255) gCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc > 1:32748/1‑151 (MQ=255) cATAAACCGGCACCTCATCCTCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTAATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGAtt < 1:423154/142‑1 (MQ=255) aCCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTg > 1:763926/1‑151 (MQ=255) ccgaccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCATCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGCCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATg < 2:823083/137‑1 (MQ=255) ccGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt < 1:436175/149‑1 (MQ=255) | CGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTG > NC_000913/298810‑299097 |
| 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: ‑ |