| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 298,960 | 0 | A | C | 43.2% | 56.3 / 15.5 | 44 | Y659D (TAT→GAT) | paoC | PaoABC aldehyde oxidoreductase, Moco‑containing subunit |
| Reads supporting (aligned to +/- strand): ref base A (23/2); new base C (0/19); total (23/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.49e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACAC > NC_000913/298810‑299107 | cGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATa < 1:999720/151‑1 (MQ=255) ttATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATAcc > 1:275501/1‑151 (MQ=255) cacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCg < 2:609987/149‑1 (MQ=255) cacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGcc > 1:187201/1‑151 (MQ=255) ttGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGt > 2:229203/1‑151 (MQ=255) ggtgATAGTCGCGCTCACGCCGCACAGGCCCACCTCACCCACACCTTTGGCCTTCATCGGGAAGGATATGGCGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGCG‑ATATCCGCTGGAACCGACACCTCATACCCCGCCATAACGTg > 2:137301/4‑150 (MQ=255) cGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGAtt > 1:317359/1‑151 (MQ=255) ccGCGCTCACCCCGCACGGCCCCAGCTCCCCGACCCCTTTGGCTTTCATGGGGGGGGTTATGGGTTCGGTACCACCCAGAAAACCCCCCTCCTGTTTGGGT‑ATACCCGAATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTaa < 2:939734/147‑1 (MQ=255) cgcgCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACg > 2:18001/1‑148 (MQ=255) gcTCACGCCGCACAGGCCCAGCTCACAGACACATTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATGATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACTCCCCCATATCGTGATTAACGAAg > 2:547360/1‑148 (MQ=255) cACGCCGCACAGGCCCAGCTCACCGACCCCTTGGGCCTTCATCGGGGAGGATATGGGGTCGGTATCACCAAGGAAAATCCCCTCCTGTTTTGGG‑ATATCCGCATGACCCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGcc < 1:164503/150‑1 (MQ=255) ccgcAGGGGCCCAGCTCACCGACCCCTTGGGCCTTCAGGGGGGAGGATATGGGGTCGGTATCACCCAGAAAAATCCCCTCCTTTTTGGGG‑ATACCCGCAAGACCCGGCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCaaa < 2:1003612/150‑1 (MQ=255) cgcACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACg > 1:726675/1‑151 (MQ=255) caGGCCCAGCTCCCCGACCCCTTTGGCCTTCAGGGGGGAGGATGTGGGGTCGTTATCATCAAGAAAAACCCCCCCCTTTTTTGGA‑ATATCCAAATAACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc < 2:857799/150‑1 (MQ=255) caGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc > 2:903039/1‑150 (MQ=255) aGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAt > 2:706426/1‑151 (MQ=255) ggCCCAGCTCCCCGACTCCTTTGGCCTTCAGCGGGGGGGATATGGGGTCGGTATCATCCAGGAAACTCCCCTCCTTTTTTGGA‑ATATCCGCATGACCCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAt < 1:402537/150‑1 (MQ=255) gACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc > 2:812665/1‑151 (MQ=255) acccttttGCCTTTCATGGGGGAGGATATGGGGTCGGTATCACCCAGAAAAACCCCCCCCTTTTTTGGG‑ATACCCGCATGACCCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc < 1:301836/145‑1 (MQ=255) tttGGCTTCCATCGGGGAGAATATGGGGTCGGTTTCATCAAGAAAAACCCCCTCCTTTTTTGGG‑ATATCCGCCAGACCCGGCCCCTATTCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCAt < 1:644388/150‑1 (MQ=255) ccTTCAGGGGGGAGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCCCCTGTTTTGGG‑ATACCCGCATGACCCGCCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTg < 1:381859/150‑1 (MQ=255) tGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGt > 1:480146/1‑151 (MQ=255) ggggTGGTTATCCCCCAGAAAAACCCCCCCCTTTTGGGGG‑ATACCCAAATGCCCGGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTGGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat < 2:952046/147‑1 (MQ=255) gttaacatcCAAAAAAACCCCCTCTTTTTTGGGG‑ATACCGACGAGACCCGGCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTg < 2:82313/144‑1 (MQ=255) aaaaCCCCCCCCTTTTTGGGT‑ATCCCCGCATGACCCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCaa < 2:372220/142‑1 (MQ=255) aaaTCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCt > 2:832733/1‑151 (MQ=255) aaccccctcctTTTTTGGG‑ATACCCGCATAACCCGGCCCCTCTTCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTg < 1:146123/146‑1 (MQ=255) cctcctGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc > 1:451355/1‑151 (MQ=255) tGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt > 2:93862/1‑151 (MQ=255) ttttGGA‑ATACCCGCATGACCCGCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt < 1:215410/150‑1 (MQ=255) tttGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCACCAATGACCTGGCTGCGCGCAGTTTTCg > 2:95505/1‑151 (MQ=255) tttGGG‑ATACCCGCATAACCCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTc < 1:252221/150‑1 (MQ=255) ttGGGAATACCCGCAGGACCCGCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTc < 1:873095/150‑1 (MQ=255) ggg‑ATACCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg < 1:990078/149‑1 (MQ=255) atatCCTAATGAACAGGCACCTCACACTCCGCAATATCGTGATTAACGACATAGCCAAGACGGCCAGCCACCGCCAGCTCCTCCATCAGTGCCGAGACCATCCTCATAGTCATCGCGCCAATGACCTGGCTGAGCGTGGTTTTCGTATTc > 2:619465/1‑150 (MQ=255) tacccacagaacccGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGAtt < 1:547992/136‑1 (MQ=255) atCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGACGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg > 2:845702/1‑151 (MQ=255) cccGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGa < 1:436189/150‑1 (MQ=255) aaCCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt > 1:226014/1‑151 (MQ=255) cccgcccccTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt < 2:473637/143‑1 (MQ=255) ccgcccccTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTg < 2:296023/144‑1 (MQ=255) cGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc > 1:478914/1‑150 (MQ=255) aCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgcagca > 2:995846/1‑150 (MQ=255) tCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacac > 2:543059/1‑150 (MQ=255) tACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacacac > 1:223258/1‑149 (MQ=255) | CGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACAC > NC_000913/298810‑299107 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |