| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 298,960 | 0 | A | C | 42.2% | 67.3 / 17.2 | 45 | Y659D (TAT→GAT) | paoC | PaoABC aldehyde oxidoreductase, Moco‑containing subunit |
| Reads supporting (aligned to +/- strand): ref base A (22/4); new base C (0/19); total (22/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.26e-09 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGC > NC_000913/298810‑299110 | cGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATa < 2:570818/151‑1 (MQ=255) acCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGcc > 2:747358/1‑150 (MQ=255) gTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCg > 2:57977/1‑151 (MQ=255) gagaGCCCCCCTCACGCCGCACAGGCCCCGCTCCCCGCCACCTTTGCCCTTCATCGGGGAGGATTGGGGGTCGGTTTCTCCCAGAAAAACCCCCTCCTGTTTGGGTATAT‑CCGCGTGAACCGGCACCTCACACCCCGCCATATCGTGATTa < 1:503493/148‑1 (MQ=255) aGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACg > 2:338187/1‑151 (MQ=255) gcgcTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGt > 2:506175/1‑151 (MQ=255) cgcTCCCGCCGCACGGCCCCAGCTCCCCGACACTTTTGCCTTTCATGGGGGAGGATAGGGGGTGGGTATCACCCAGGAAAATCCCCGCCTTTTTGGGGATAT‑CCGCATGACCCGCCCCCTATTCCCCCGCCATATCGTGATTAACGa < 1:361161/147‑1 (MQ=255) cTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGc > 1:126860/1‑151 (MQ=255) tCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGcc > 1:426036/1‑151 (MQ=255) gccgcACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATAGTGATTAACGAAGTAGCCCa > 1:138879/1‑149 (MQ=255) ccgcACGGCCCCAGCCCCCCGCCTCTTTTGGCCTTCAGGGGGGGGGATGTGGGGTCGGTATCATCCAGAAAAACCCCCCCCTTTTTGGGAATAC‑CGACATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCa < 1:272279/148‑1 (MQ=255) ccgcACAGGCCCAGCTCCCCGACCCTTTTGGCCTTCAGCGGGGAGGATATGGGGTCGGTATCACCCAGAAAAATCCCCTCCTGTTTGGGGATAC‑CCGCATGAACCGGCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc < 1:266055/151‑1 (MQ=255) ggCCCAGCTCACCGCCCCTTTTGGCCTTCAGCGGGGAGGATATGGGTTCGGTATCATCCAGAAAAATCCCCTCCTGTTTTGGGATAC‑CCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCa < 1:926016/149‑1 (MQ=255) ccAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCAc > 2:930368/1‑151 (MQ=255) gCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGc > 1:853901/1‑151 (MQ=255) cTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGcc > 1:910951/1‑151 (MQ=255) tCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGcc > 1:932606/1‑150 (MQ=255) cccGACCCCTTTGCCTTTCGTGGGGGGGAATGGGGGGTCGGTATCACCCAGAAAAACCCCCCCTTTTTTTGGTATAC‑CCGCATAACCCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCa < 1:565974/148‑1 (MQ=255) agACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctc > 1:335796/2‑151 (MQ=255) cccttttGGCTTTCATGGGGGAGGTTAGGGGGTCGTTATCA‑CCAAGAAAATCCCCCCCTTTTTTGGGGTTTCCCGCATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc < 2:530301/145‑1 (MQ=255) cccttttGGCCTTCAGCGGGGAGAATAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTTTTTTGGGATCT‑CCGCATAACCCGCCACCTCTTCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcctc < 2:347985/147‑1 (MQ=255) gCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTg > 2:780119/1‑151 (MQ=255) ttCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcg > 1:923614/1‑151 (MQ=255) tCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATAccc < 1:341264/85‑1 (MQ=255) tCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATAccc > 2:341276/1‑85 (MQ=255) aTCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCgcgc > 1:850576/1‑151 (MQ=255) ggggTCGTTATCACCCAGAAAAACCCCCCCCTTTTTTGGGATAC‑CCGCAGAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgccc < 2:773542/145‑1 (MQ=255) ggggTCGGTATCATCCAGAAAAACCCCCTCCTTTTTGGGGATAC‑CCGAATGACCCGGCACCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATa < 2:53438/148‑1 (MQ=255) gggccgggTACCATCCAGAAAAAACCCCTCCTGTTGGGGAATCC‑CCGATTGACCCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc < 2:343170/145‑1 (MQ=255) gggTCGGTATCACCCAGAAAAATCCCCTCCTGTTGTGGGATAT‑CCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCa < 2:31632/151‑1 (MQ=255) gTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcgc > 2:632600/1‑151 (MQ=255) atcatcCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGcc > 1:863495/1‑150 (MQ=255) ataaaccggaaaaaCCCCCTCTTTTTTGGGAATTT‑CCCTATGAACGGCCCCCTTATACCCCGCCATATCGTGATTAACGAAGTAGCCCAATCGGTCATCCACCGCCAGCGCCTCCATCAGTGCCGCGCCCATGCCCATATTCATTGCGcct < 1:497210/141‑2 (MQ=255) ggAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTg > 1:117719/1‑151 (MQ=255) gaaaaaTCCCCTCTTGTTTTGGTATCC‑CCGAATGACCCGCCCCCTCTTCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTg < 2:476619/149‑1 (MQ=255) ccccccctTTTTTGGGATAC‑CCGCAGGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgc < 2:783988/144‑1 (MQ=255) tttgGGAATAT‑CCGCATGACCCGGCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt < 2:92732/150‑1 (MQ=255) ttGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg > 2:491468/1‑151 (MQ=255) ggaatac‑CCGCATAACCCGTCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg < 2:775821/146‑1 (MQ=255) ccgaaggAACCGCCCCCCCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTc < 1:233559/140‑1 (MQ=255) acccGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGTCCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCg < 2:656227/145‑1 (MQ=255) acccGGCACCTCCTCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATAACCTGGCTGCGCGCAGTTTTCGGATTCAGGATg < 2:743763/143‑1 (MQ=255) ccgcccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt < 2:996781/143‑1 (MQ=255) tCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacac > 2:423540/1‑150 (MQ=255) cATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacaca < 2:512614/150‑1 (MQ=255) cATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacacac > 2:734335/1‑151 (MQ=255) aCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGc > 2:3829/1‑151 (MQ=255) | CGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATAT‑CCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGC > NC_000913/298810‑299110 |
| 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: ‑ |