Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 298,960 | 0 | A | C | 42.6% | 80.8 / 26.0 | 55 | Y659D (TAT→GAT) | paoC | PaoABC aldehyde oxidoreductase, Moco‑containing subunit |
Reads supporting (aligned to +/- strand): ref base A (28/3); new base C (0/23); total (28/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.39e-12 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACAC > NC_000913/298813‑299107 | tataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATAcc < 1:412212/150‑1 (MQ=255) cacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGcc > 1:69924/1‑151 (MQ=255) ccGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCat > 1:748881/1‑151 (MQ=255) cACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGccc > 1:202172/1‑151 (MQ=255) aCGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGccc > 2:1033440/1‑150 (MQ=255) aCGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGccc > 1:587551/1‑150 (MQ=255) cgccgcACGGCCCCAGCCCCCCGACCCCTTTGGCCTTCATGGGGGAGGTTAGGGGGTTGGTATCACCCAGGAAAACCCCCTCCTG‑TTTTGGGATATCCGCATGAC‑CCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCaa < 1:621770/151‑1 (MQ=255) ccgcACAGGCCCAGCCCCCCGACCCCTTTGCCCTTCATGGGGGAGGTTAGGGGGTCGGTATCACCCAGAAAAATCCCCTCCTG‑TTTGGGAATACCCGCATGAA‑CCGGCCCCTCATCCCCCGCCATATCGTGATTAACGa < 2:926561/139‑1 (MQ=255) ccccaccgcccccGCTCCCCCACCCCTTTGGCCTTCATCGGGGAGGATAGGGTGTCGGTATCATCCAGAAAAATCCCCTCCTG‑TTTTGGGATACCCGCATCAA‑CCGACACCTCATTCCCCGCCATATCGTGATTAACGAAGTAGCCCaaa < 2:124338/141‑1 (MQ=255) cgcACAGGCCCAGCTCACCGACCCCTTTGGCCTTCAGCGGGGAGGATAGGGGGTGGGTATCACCCAGAAAAACCCCCTCCTG‑TTTTGGGATACCCGCATGAA‑CCGGCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACg < 1:371704/151‑1 (MQ=255) cacaGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGt > 1:944291/1‑151 (MQ=255) ggCCCAGCTCACCGACCCCTTTGGCCTTCACGGGGGAGGATATGGGGTGGGTATCACCCAGGAAAACCCCCTCTTT‑TTTTGGGATATCCGCAGGAC‑CCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAt < 1:70652/150‑1 (MQ=255) tCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGcc > 2:362125/1‑150 (MQ=255) ccccGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGTTATCACCCAGAAAAATCCCCTCCTG‑TTTTGGGATATCCGCATGAC‑CCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCa < 1:369093/148‑1 (MQ=255) cACCGCCCCCTTTGGCTTTCATCGGGGAGGATAGGGGGTGGGTATCACCCAGGAAAATCCCCTCCTT‑TTTTGGGATACCCGCATGAC‑CCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCg < 1:75221/147‑1 (MQ=255) ccGACCCCTTTGGCCTTCATGGGGGAGGATATGGGGTCGGTATCACCCAGAAAAACCCCCCCCTG‑TTTTGGAATCACCGCAGAAC‑CCGCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGct < 2:652917/151‑1 (MQ=255) cGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctc > 2:918824/1‑151 (MQ=255) cacCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcctc > 2:563079/1‑151 (MQ=255) tttGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc > 2:631972/1‑151 (MQ=255) ttGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc > 1:58036/1‑150 (MQ=255) cctacacggggaAGGTTAGGGGGTC‑GTTTCCCCCAAGAAAACCCCCCCCTTTTTTTGGAATACCCGAATAAC‑CCGCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc < 2:150369/143‑1 (MQ=255) cTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGcc > 2:1058432/1‑151 (MQ=255) ttcaggggggggGATATGGGGTCGTTACCACCCGGAAAAACCCCCTCCTG‑TTTGGGGATACCCCCAGGAC‑CCGCCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc < 1:354197/143‑1 (MQ=255) tCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTGACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcc > 2:484921/1‑151 (MQ=255) tCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcc > 2:173404/1‑151 (MQ=255) ggggAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccc > 2:5595/1‑150 (MQ=255) gggAGGATGGGGGGTCGTTATCACCAAGAAAACCCCCCTCTTT‑TTTGGGGATCCCCACATGAC‑CCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat < 2:133382/151‑1 (MQ=255) gaggaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgc > 1:511531/1‑151 (MQ=255) gaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgccc > 2:379773/1‑150 (MQ=255) gggtgggttaccacccAGAAAAACCCCCGCCT‑‑TTTTGGGAATACCGGAAGAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc < 2:216147/137‑1 (MQ=255) tATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGcc > 1:662357/1‑151 (MQ=255) gtaccccctagAAAATCCCCCCGTT‑TTTTGAGATATCTGAAGGAG‑CGGACACCCTCTCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAACCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcg < 1:505619/139‑1 (MQ=255) tcatcCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCa > 2:534512/1‑150 (MQ=255) tcatcCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCa > 1:659281/1‑150 (MQ=255) cAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAcc > 2:463203/1‑151 (MQ=255) gaaaaaTCCCCTCCTT‑TTTGGGGATATCCGAATGAA‑CCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGa < 1:401903/145‑1 (MQ=255) aaaaTCCCCTCTTT‑TTTGGGAATACCCGAAGAAA‑CCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt < 1:429037/148‑1 (MQ=255) aaTCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTg > 2:1062577/1‑151 (MQ=255) cccctcttT‑TTTTGGGATCCCCGAAGGAC‑CCGCCCCCCCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg < 2:818136/148‑1 (MQ=255) tG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt > 1:877243/1‑151 (MQ=255) g‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt > 1:743365/1‑151 (MQ=255) ttttGGGATACCCGAATAAC‑CCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtt < 1:1058751/148‑1 (MQ=255) ttggggATACCCGAATGAC‑CCGCCCCCACACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTc < 2:511380/147‑1 (MQ=255) tggggATCCCCGAATAAA‑CCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCg < 1:233364/148‑1 (MQ=255) acccGCATGAC‑CCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCa < 1:1033761/147‑1 (MQ=255) tCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg > 1:286153/1‑150 (MQ=255) cccGCATGAC‑CCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTc < 2:105024/146‑1 (MQ=255) ccGCATAAC‑CCGCCACCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAg < 1:191389/148‑1 (MQ=255) gCAGGAC‑CCGACCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATg < 1:1170893/150‑1 (MQ=255) tGAC‑CCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATg < 2:351747/147‑1 (MQ=255) ccgcccccTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt < 2:17030/143‑1 (MQ=255) cGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc > 1:290913/1‑150 (MQ=255) aCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac > 1:242859/1‑151 (MQ=255) ccTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcaca < 1:432010/151‑1 (MQ=255) aTACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacacac > 1:879650/1‑150 (MQ=255) | TATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTG‑TTTTGGGATATCCGCATGAA‑CCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACAC > NC_000913/298813‑299107 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |