Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,386,862 | 0 | G | A | 63.6% | ‑3.0 / 15.3 | 22 | intergenic (+31/‑72) | yfbL/yfbM | putative M28A family peptidase/DUF1877 family protein |
Reads supporting (aligned to +/- strand): ref base G (1/7); new base A (14/0); total (15/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.69e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-03 |
TTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAA > NC_000913/2386717‑2386972 | tttttATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACa < 1:139619/151‑1 (MQ=255) aaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACattttattt < 1:290263/151‑1 (MQ=255) tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCA‑TCTGTTTTAAGa < 1:130958/151‑1 (MQ=255) tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCA‑TCTGTTTTAAGa < 2:256062/151‑1 (MQ=255) tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCA‑TCTGTTTTAAGa < 1:360121/151‑1 (MQ=255) aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGAAATGCACATTTTTTTTTTCCTCA‑TCtttttttaaaaggg > 1:205061/1‑143 (MQ=255) aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTTTTTTTCctcc‑tc > 2:80650/1‑137 (MQ=255) aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCA‑Tc < 1:80648/137‑1 (MQ=255) cagacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGGAAAAAGGAAATGCTCAATTTATTTTTTCTCATTCTGtttttaaagcggaa > 1:227579/1‑140 (MQ=255) ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGAAATGAACATTTTATTTTTTCTCA‑TCTGTTTTTAAAAGC‑GAAACATG‑AAAATAAAAGGAGGGAAAAAAAATAGGa > 1:447668/1‑151 (MQ=255) ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAGCGAAAATCAAATTTTTTTTTTTCCCA‑TTTGTTTTTAAAAGA‑GAAAAATG‑AAATTTAACGGAGGAAAAAAAAATgggg > 1:484926/1‑150 (MQ=255) ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGACAAGGAAAATTTTTTTTTCCTCC‑TCTGTTTTTAGATGT‑GAAAAATG‑AAAATTAACCGGGGAAAAATAAATgggg > 1:157446/1‑150 (MQ=255) ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATGGACATTTTTTTTTTCCTCA‑TCTGGTTTAAGATGC‑GGAACATG‑AAAATCAAAGGGAGAAAAAAAAATgggg > 2:272017/1‑150 (MQ=255) ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGATTCATAAAAAAAGACATACCAATTTTTTTTTTTCTCA‑TCTGTTTTTAGATTG‑GAAACATG‑AAAATCTACGGAGGGGAAATAAATaggg > 1:127686/1‑150 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAACGACAAGAACATTTTATTTTTCCTCA‑TCTGTTTTTAGATGC‑GAAAAATG‑AAATTTAAAAGAGGAAAAAAAAAATGGAATGATTGGCTa > 1:199273/1‑151 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGAAATAAACAATTTATTTTTCCTCC‑TCTGTTTTTAAAAGA‑GGAAAATG‑AAAATTAAATGAGGGAAAATAAATGGGGATTATTgggtt > 1:324893/1‑148 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCA‑TCTGTTTTTAGATGG‑GAAAAATG‑AAAATTAAAGGGGGGAAAAAAAATGGGGATGATTGGGTa > 2:447352/1‑151 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGGAATGCACAATTTATTTTTCCTCA‑TCTGTTTTTAAAAGCGGAAAAATA‑AAATTAACGGGGGGAAAAAAAAATGGG‑ATGATTggttt > 2:203900/1‑148 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAAAAAAAAAGAAATGCACATTTTATTTTTTCTCC‑TATGTTTTAAGAAGA‑GAAAAATGAAAATTTAAAGGGGGAAAAA‑AAATTGGAAAGAATGGGTa > 1:550217/1‑151 (MQ=255) gtgtataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCA‑TCTTTTTTTAGATGC‑GAAAAATG‑AAAATCAACGGAGGAAAAATAAATGGGAATGATTGGCTAt > 2:366477/4‑151 (MQ=255) ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCCCATTTTATTTTCCCTCA‑TCTGTGTTAAGATGC‑GAAACATT‑GAATTCAACGGAGGAAAAATATATGGTAATGATTGGCTATTTTGCAGAAATTg < 1:214278/151‑1 (MQ=255) caGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACAAGCACATTTTTTTTTTTCTCC‑TCTGTTTTTAGATGC‑GAAAAATG‑AAAATCTAAGGAGGGAAAAAAAATTGGGATGATTTGGTATTTTTCAGAAAATTACTCAGaa > 1:249617/1‑151 (MQ=255) | TTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAA > NC_000913/2386717‑2386972 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |