Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,386,862 | 0 | G | A | 48.6% | 5.7 / 34.0 | 37 | intergenic (+31/‑72) | yfbL/yfbM | putative M28A family peptidase/DUF1877 family protein |
Reads supporting (aligned to +/- strand): ref base G (4/15); new base A (18/0); total (22/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.01e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAA > NC_000913/2386725‑2387004 | gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTCGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGACAAGCACatt > 1:203310/1‑151 (MQ=255) taataaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACattt < 2:312375/151‑1 (MQ=255) ataaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACattttt > 2:86860/1‑150 (MQ=255) cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTc < 2:215461/151‑1 (MQ=255) cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACatttttttttt > 2:244137/1‑150 (MQ=255) aTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTccc < 2:202599/151‑1 (MQ=255) cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGTTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGttt < 2:139079/151‑1 (MQ=255) gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGGACATTTTTTTTTTCCTCATCTGtttttaa > 2:82235/1‑148 (MQ=255) gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTTTTTTTTCTCCTCTGTTTTTAg > 2:79368/1‑151 (MQ=255) tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa < 1:183120/151‑1 (MQ=255) tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa < 1:28305/151‑1 (MQ=255) tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa < 2:256951/151‑1 (MQ=255) aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAAAGGCCCTTTTTTTTTTCCTCCTCTGTTTTTAGATGCg > 2:95678/1‑151 (MQ=255) acagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAAc < 1:320402/148‑1 (MQ=255) aTTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATGCACATTTTTTTTTTCCTCATCTGTTTTTAGATGCGAAAAAAG‑AAATTCaa > 2:26726/1‑150 (MQ=255) aTCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGCACATTTTATTTTTCCTCATCTGTTTTTAAATGCGGAACATG‑AAAATCAACggg > 2:357338/1‑150 (MQ=255) ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGACAAGCACATTTTTTTTTTTCTCATCTGTTTTTAGAAGCGAAAAATG‑AAATTCAACGGA‑GGGAAAAAAAATgggg > 1:135258/1‑150 (MQ=255) ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTATTTTTTCTCATCTGTTTTAAGATGCGAAACATG‑AAAATCAACGGA‑GGAAAAATAAATgggg > 1:146370/1‑150 (MQ=255) gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGTAAAAAAGAAATGATCATTTTTTTTTTTCTCATCTGTTTTTAGAAGCGAAACATG‑AAAATATATTGA‑GGAAAAATAAATGGGaa > 1:361825/1‑151 (MQ=255) gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGACATGCACATTTTATTTTTTCTCATCTGTTTTTAGATGCGAAACATG‑AAATTTAACGGG‑GGAAAAATAAATgggga > 1:296710/1‑149 (MQ=255) gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGAAAAAGAAATTTTTTTTTTTCTCATCTGGTTTTAGAAGAGAAAAAAGAAAATTAACCGGG‑GGAAAAAAAAATGGGa > 1:235569/1‑151 (MQ=255) gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGAAAAGCAAATTTTTTTTTTTCTCATCTGTTTTTAGATGCGAAAAAAG‑AAAATTAAAGGA‑GGAAAAAAAAATGGGGATGa > 1:400726/1‑151 (MQ=255) gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATGCACATTTTATTTTTCCTCATCTGTTTTTTGATGCGAAAAAAG‑AAAATTAACGGA‑GGAAAAAAAAATGGGAAAGa > 2:328628/1‑151 (MQ=255) tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGA‑GGAAAAATAAATGGGAATGAt < 1:426723/151‑1 (MQ=255) tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGA‑GGAAAAATAAATGGGAATGAt < 1:317702/151‑1 (MQ=255) gATGGTGTTATAACTGTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAGAAGCGACATGCACATGTTATTTGCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGA‑GTAAAAATAAATGGGAATGAtt < 1:391694/150‑1 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTTTTTTTTCTCATCTGTTTTTAAATGCGAAACATG‑AAAATCAAAGGA‑GGGAAAAAAAATGGGGATGATTGGGTa > 1:445158/1‑151 (MQ=255) gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTATTTTTTCTCCTCTGTTTTTAGATGCGAAAAAAG‑AAAATCAAAGGA‑GGGAAAAAAAATGGGGATGATTGGGTATTTTGCAGaaaat > 1:12400/1‑149 (MQ=255) ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACATGCACAATTTTTTTTTCCTCCTCTGTTTTAAGATGCGGAAAATT‑AAATTTAACGGAGGGAAAAAAAAATGGG‑ATGATTTGCTATTTTTCAGAAAATg > 2:321078/1‑151 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACATGCACATTTTATTTTTTCTCATCTGTTTTTAGATGCGAAAAAAG‑AAAATTAAAGGA‑GGGAAAAAAAATAGGGATGATTGGCGATTTTTTAGAAAATGAATCAGa > 1:240171/1‑151 (MQ=255) aTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGA‑GGAAAAATAAATGGGAATGATAGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTg < 1:74536/151‑1 (MQ=255) cGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAg < 1:102898/151‑1 (MQ=255) gCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc < 2:194043/149‑1 (MQ=255) ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa > 2:348290/1‑148 (MQ=255) ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa > 2:182896/1‑148 (MQ=255) ccAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaa < 2:346301/151‑1 (MQ=255) cAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaaa > 1:282847/1‑151 (MQ=255) | GTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAA > NC_000913/2386725‑2387004 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |