| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,027,233 | 0 | T | A | 27.8% | 62.2 / 10.1 | 36 | intergenic (+5/+376) | fre/fadA | NAD(P)H‑flavin reductase/3‑ketoacyl‑CoA thiolase |
| Reads supporting (aligned to +/- strand): ref base T (20/6); new base A (0/10); total (20/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.15e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTT‑AACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTAT‑ATTGC‑CGGACGTTTTGAGATGGCGAAAATTGCCCGC‑GATCTGTTTTGCAG‑TGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACC > NC_000913/4027014‑4027469 | aaGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTT‑AACGGCGGGCTTGCCGGATCACGGGACGCCGGCAGAGGCTGATATCTATTATT‑C‑CGGGCGGTTTGTGATGGCGGAAAATGCCCCC‑GACCTGTTTTTCAG‑TGGGCGCAATGTGCGGGGAGAAAGCCTGTTTGGGGAGGGGGTTGCATTTATCTTGGATATAAAAAAA‑ccccccccctccagggggggaga > 2:91289/1‑235 (MQ=255) gggtggAACAACGGGAAGGGGGTGGGGGTGGGCGTACTGGCACCGTGTT‑AAGGGCGGTATGGCTGGATCACGGTACGCTGGCAGAGAATGATATTTAT‑ATTGC‑GGGACGTTTTGAGAGGGCGAAAATTGCCCGC‑GATCTGTTTTGCAG‑GGGGCGTAATGCGCGGGAATATCGCTTGTTTGGCGATGCGTTTGAATTTATCTGAGATAAAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAAc < 1:23320/246‑1 (MQ=255) gggggggggCGTACTGTCCCCGTGTT‑AAGGGCGTTAGTGCAGGATCACGGTACCCTGGCAGAGCATGATCTCTAT‑ATTGC‑CGGACTTTTTGAGTTGGGGAAAATTGCCCAC‑GATCTGTTTTGAAG‑TGAGGGTAATGCGCGGGAAGATCGTGTTTTTGGCGTTGTGTTTGGATTTATCGGAGAAATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAGCTAAACTGTTATTCAGTGGCATTTAGAt < 1:454855/244‑1 (MQ=255) ggCGTGGGCGTACTGGCACCGTGTT‑AACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTAT‑ATTGC‑CGGACGTTTTGAGATGGCGAAAATTGCCCGC‑GATCTGTTTTGCAG‑TGAGCGTAATGCGCGGGAAGATCGCCTGGTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGc < 1:22911/196‑1 (MQ=255) gggggCGTACGGGCAC‑GTGTTTAACGGCGGTATTGCAGGATCACGGTACGCGGGCAGAGCATGATATCTAT‑ATTGC‑CGGACGTTTTGAGATGGCGAAAATTGCCCGC‑GATCTGTTTTGCGG‑TGAGCGTAATGCGCGGGAAGATCGCTTGTTGGGCGATGGGTTTGCATTTATCTGAGATATAAAAAAC‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATg < 1:141936/249‑1 (MQ=255) cacggtgtA‑AACGGCGGATTTGCAGGATCACGTAACTGTGACAAAGTATTATATCTAT‑ATTGC‑GGGCGGTTTTGAG‑GGGGGAAAAATGCCCCCCGATGTGTTTTGCAG‑GGAGCGAAGTGCGGGGGAAGACCGCCGTTTTGGCGATGCGTTTGCATTTATC‑GAAAAAAAAAAAAAGCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATGCAGTGGCATTTAGATCTATGACGTATCTGGCaa < 2:246298/247‑1 (MQ=255) ggtgtT‑AACGGCGGTATGGGAGGATCACGGTACGCTGGCAGAGAATGATTTCTAT‑ATTGC‑GGGACGTTTGGAGAGGGCAAAAATTGCCCGC‑GATCTTTTTGGCGG‑GGAGGGTAATGGGGGGGAAGAGCGCCGGTGGGGGGATGCGTTGGCATTTATCGGAGATAAAAAAAAC‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAAcc < 2:207258/250‑1 (MQ=255) aaaCGGCGGTATTGCAGGATCACGGTACGCGGGCAGAGTATGATATCTAT‑ATGGC‑CGGCCGTTTTGAGATGGCGAAAATTGCCCGC‑GATCTGTTTTGCAG‑TGAGCGTAATGCGCGGGAAGATCGCTTGTTGGGCGATGCGTTTGCATTTATCTGAGATAAAAAAAAC‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGc < 1:186565/249‑1 (MQ=255) aGGATCACGGTACGCTGGCAGAGTATGATATTTAT‑ATTGC‑CGGACGTTTTGAGAGGGCGAAAATTGCCCGC‑GATCTGTTTTGCAG‑TGAGCGTAAGGCGCGGGAAGATCGCCTGTTTGGCGATGGGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTc < 1:12512/251‑1 (MQ=255) aTCACGGTACGCTGGCAGAGCATGATATCTAT‑ATTGC‑CGGACGTTTTGAGATGGCGAAAATTGCCCGC‑GATCTGTTTTGCAG‑TGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGGAGAACGGCAACTAAACTGTTCTTCAGGGGCATTTAGGTCTATGGCGTATCTTGCAAACCATGCCCGCTGCggccccgtcgcg > 1:223708/1‑240 (MQ=255) cccggtccggtgacAGAGCTTGATATCT‑T‑TTTGCGCGGACGTTTTGAGGTGGAGAAAGTTGCCCGC‑GTTTTGTTTGGCAG‑TGGGCGTAATGCGCGGGAAGCCCGCTTGTTGGGCGGGGCGTTTGCATTTATGTGAGATAAAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTc < 2:439713/238‑1 (MQ=255) cACGGTACGCTGGCAGCGAACGATATCTAT‑ATTGC‑CGGACGTTTTGAGATGGGGAAAATTGCCCGC‑GATCTGTTTTGCCG‑TGAGAGTAATGCGCGGGAAGATCGCCTGATTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAAACGTTATTTAGGGGCATTTATGTCTATGACGTATCTGGCAAAACATCCCCGATGCGACGATGTCGTGTc > 2:341764/1‑251 (MQ=255) agagCATGATATCTAT‑ATTGC‑CGGACGTTTTGAGATGGCGAAAATTGCCCGC‑GATCTGTTTTGCAG‑TGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCCTTTAGATCTATGACGTAACTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCTTGCCTAc > 2:192513/1‑250 (MQ=255) gagCATGATATCTAT‑ATTGC‑CGGACGTTTTGAGATGGCGAAACTAGCCCGC‑GATCTGTTTTGCAG‑TGAGCGTAATGCGCGGGAAGCTCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGCACGGCCACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCCAACCATGCCCCATGCGACGCTGTCGCGTCTTATCGTGCCTaaa > 1:252034/1‑248 (MQ=255) atgttGC‑TGGACTTTTTGAGATGGAGAAAATGGCCTGG‑GAT‑TTTTTGGGAGGTGAGTGAAATGCGCGGAAAGACCGCCTGTTTGGCGAGGCGTTTGCATTTATCGGAAAGATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTAGCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGaa < 2:179028/245‑1 (MQ=255) gttctgagacgcggAAAATTGCCCGC‑GATCTGTTTTCCAG‑TACGTTGAACGGGCGAACAGATCGCCTTTTTGGCCATGTTTTTGCATTCATCTGCTTAATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAATCTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATaa < 1:103330/238‑1 (MQ=21) ttggagagggagaaaatgccccGC‑GATCTTTTTTGCAG‑TGAGCGTAATGCGCGGAAAGATCGCTTGTTGGGCGATGCTTTTGCATTTATCTGAGATAAAAAAAAC‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAg < 1:132370/231‑1 (MQ=38) ttGAGATGGCGAAAATTGCCCGC‑GATCTGTTTTGCAG‑TGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGGTAAGGc > 2:140751/1‑251 (MQ=255) gagggagAAAATTGCCCGC‑GATCGTTTTTGCAG‑GGGGGGTAAGGCGCGGGAAGATCGCCTTTTTGGGGATGCGTTTGCTTTTATCGAAAATAAAAAAAAC‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGttt < 1:261588/245‑1 (MQ=37) gaTGGCGAAAATTGCCCGC‑GATCTGTTTTGCAG‑TGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGttt > 1:176871/1‑251 (MQ=255) aaaaTTGCCCGC‑GATCTGTTTTGAAG‑TGAGGTTAGGGCGGGGAAAGATCGCGTTTTGGGCGGTGCGTGTGAATTGATCGGAAATAAAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGGGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTAcgccgc < 1:274154/251‑1 (MQ=25) ccGC‑GATCTGTTTTGCAG‑TGAGCGTAATGCGCGGGACGATAGCCTGTTTGGCGATGCGTTTGCCTTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGAGGGGAGCACGGCAACAAAACTGTTATACAGTGGCATTTAGATCTATGACGCATCTGGCAAACCATGCCCGATGCGACGCCGTCGCGGCTTATCGTGCCTACAAATAGTCCGCACCGTAGGCCGGATAAGGGGTTTTCGCCGCATCCGGc > 1:528035/1‑250 (MQ=255) gATCTGTTTTGCAG‑TGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGCCGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGCAGGCCGGCTAAGGCGTTTACGCCGCATCCGGCaactg > 1:445169/1‑248 (MQ=255) gATCTGTTTTGCAG‑TGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAAATGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCCATGCGACGCTGTCGGGTCTTATCGTGCCTACAAAAAGTCCGAACCGCAAGCCTGAGAAGGGGGTTACGCCCCATCCGGccatcg > 2:192759/1‑246 (MQ=255) cTGTTTTGCAG‑TGAGAGTACTGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAAATGTTATTCAGGTGCCTTTAGATCTATGACGTATCTGGCAAACCATTCCCGATGCGACCCTGTCGCGTCTTATCGTGCCTTCAAAAAGTCCGAACCGTAGGCCGGATAAGGCGTTTTcgcaccctcccgcactttgtg > 2:285082/1‑233 (MQ=255) ttGCAG‑TGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGACGAACGGCAAATAAACTGTTATTTAGTGGCATTTAGATCTATGACGTATCTGGCAAACCAAGCCCGATGCGACGCTGTTGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTggggcctgg > 2:180421/1‑244 (MQ=38) gAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTCTCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACAGCAACTAAACTGTTATGCAGAGGCCTTTCGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGGAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGAATGATGCCTgt > 1:523759/1‑250 (MQ=37) aTGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCCTTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCAATTGGTGCATGATGCCTGATGCGACg > 1:233321/1‑251 (MQ=255) aTGCGCGGGAAGATCGCCTGTTTGGCGATGAGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTTAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCAGCTGCGACGCTGTAGCGTCTTTTCGTGCCTCCAAATAATCCGAAACGCACGACGGATAAGGCGGTTTCGCCGCACCCCGCAATTGGTGcctggtaccgtgtacgacc > 2:489625/1‑233 (MQ=255) ggcgGGAAGAGCGGCTGTTGGGCGGTGCTTTGGAATTTATCTGAGAAAAAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTg < 2:176871/250‑1 (MQ=25) cgGGAAGATCGTCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAACA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGAGGCATTTTGATCTAAGACGTATCTGGCAAACCATGCCCCATGCGCCGCTGTCGCGTATTATCGTGCCTACAAAACGTCCGAACCGGAGGCCGGGAAAGGCGTTTACGCCGCATCCGGCAATTTGTGCATGATGCCTGAAGCGACGCTGtt > 2:325398/1‑250 (MQ=255) gggAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGGATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCCATTGGTGCATGATGCCTGATGCGACGCTGTcg > 2:449219/1‑251 (MQ=255) ggAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTcgc > 1:119100/1‑251 (MQ=255) tttGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGGAACTAAACTGGTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTTTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGGGCATGATGCCTGATGCGACGCCGTCGGGGCTTATAGTGcc > 2:252008/1‑251 (MQ=25) aTTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGGTAAGGCGTTTTCGCCGCATCCGGCCATTGGGGCATGGTGCCCGGTGCGACGCTGTCGCGTCTTATCGGGCCTACaaaatgtccgggat > 2:225992/1‑240 (MQ=25) aTTTATCGGAGATAAAAAAAAC‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCg < 1:225992/247‑1 (MQ=255) | AAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTT‑AACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTAT‑ATTGC‑CGGACGTTTTGAGATGGCGAAAATTGCCCGC‑GATCTGTTTTGCAG‑TGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAA‑CCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACC > NC_000913/4027014‑4027469 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |