| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,578,143 | 0 | A | C | 48.2% | 70.0 / 24.1 | 56 | L79V (TTA→GTA) | maeB | malic enzyme: putative oxidoreductase/phosphotransacetylase |
| Reads supporting (aligned to +/- strand): ref base A (24/5); new base C (0/27); total (24/32) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.79e-11 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCC > NC_000913/2577999‑2578291 | cGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAg < 1:160529/151‑1 (MQ=255) gAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGc > 1:650722/1‑151 (MQ=255) gcgcGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCAc > 2:161532/1‑151 (MQ=255) gcgcGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCAc > 2:152401/1‑151 (MQ=255) aataatTTTGCCCGGGTCGATTTCGTCAACTCCAATGTCAAATCCACAACCCCCGGCGAATTTCTAAAACAGACCGCCCTTG‑CCT‑TCCCCCCCGGGTTGGCCGGCCGGCGCGCCATTTTGGCCAAC‑CCCCAGCACCGCCGTACCGTTagag < 1:320852/151‑1 (MQ=255) aaTTTGCCCGGGTCGATTTCGTCACCTTAAATGTAAAATACATCAACCCCGGCGATTTTCTTAAACAAAACGCCCTTC‑CCT‑CCCATCCCGGTTTTGCCTCCCAGCGCGCCAATGTTGCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCa < 1:300805/151‑1 (MQ=255) gggTGGATTTCGTAAATTTAAATGAAAAATACATCAACCCGGGCAATTTTTTTAAAAAGACCCCCTTTC‑CCT‑TCCATCCCCGTTTGCCCCGCCGGGGCGCCAATG‑TGCCTAACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGccacc < 2:930046/150‑1 (MQ=255) gTCGAGTTCTACACCTCAAATGCAAAATAC‑TCAAACCCCGGCAATTTTTTAAAACAAACCCCCTTGCCCT‑CCCACCCCCGGTTTGCCGCCCGGCGCCCAAATTTTCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAgg < 2:465848/151‑1 (MQ=255) tCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGt > 1:502147/1‑151 (MQ=255) gAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGtt > 2:701629/1‑150 (MQ=255) ttCGTAACCTTCAATGTCAAATACATCAACCCGGGCAATTTTCTAAAACAGACCGCCTTTC‑CCT‑CCCACCCCCGTTTGGCCTGCCAGCGCGCCAATGTGGCCAAC‑CCCCGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTa < 1:830321/148‑1 (MQ=255) gtcaattcaaaTGTCAAATACACAACACCCGGAGAATTTCTGAAAAGGAGCGCCTTGG‑CCT‑GCAAGCCCCGGTGTGTCAGCAACGGCGCTATTGGTGCGAAA‑CCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGtt < 1:957215/135‑1 (MQ=255) tCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat > 2:186577/1‑151 (MQ=255) attGTCAAATACATCAATCCCGGGAAATTTCTTAACCAGAACGCCCTTG‑CCT‑CCCACCACGGGTTGGCCTGCCAGCGCGCAAATGTTGCCAAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGt < 2:703082/145‑1 (MQ=255) aattaaaaaTACACCACCCCCGGCGATTTTCTAAAACAGACCCCCTTTG‑CCT‑CCCATCCCTGTTTTGCCTCCCAGCGCCCCAATGTTGCCTCC‑CCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtata < 1:724092/144‑1 (MQ=255) aaTGTAAAATACATCATCCCCGGCGAATTTCTAAACCAGACCGCCCTTC‑CTT‑CCCATCCCTGTTTTGCCTGCCGGCGCGCCATTGTTGCCAAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTatat < 2:262171/151‑1 (MQ=255) aaaaacaacacccccGGCAATTTTCTAAAAAAGACCCCCCTTG‑CCT‑CCCACCCCCGGTTTGCCTGCCAGCGCCCAATTTTTCCCTAC‑CCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg < 2:518664/136‑1 (MQ=255) aaaTACATCAACCCGGGCAATTTTCTTAAACAGAACGCCTTTG‑CCT‑CCCATCACCGGTTTCCCTGCCAGCGCGCCAATGTTCCCAAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGt < 2:185525/149‑1 (MQ=255) aaTACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAgg > 1:712050/1‑151 (MQ=255) acaacaACCCGGGCGAATTTTAAAAACAG‑ACCCCCTTC‑CCTCTCCATCCCCGTTTTCCCCGCCACGGCACCAATTTTGCCCAA‑CCCCAGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCt < 1:241332/150‑1 (MQ=255) aCATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtt > 1:8169/1‑151 (MQ=255) cccGGCGATTTTTTAAAACAGACCGCCCTTC‑CCT‑CCCACCCCCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACg < 1:310911/149‑1 (MQ=255) cccGGCGAATTTCTTAAACAGACCCCCCTTC‑CCT‑CCCACCCCCGGTTTCCCTGCCGGCGCCCCATTGTTCCCAAC‑CCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAAc < 2:579065/148‑1 (MQ=255) cGGCGAATTTCTTAAACAGAACGCCCTTG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGt > 1:38157/1‑150 (MQ=255) ggCGATTTTCTAAAACAGACCGCCTTTC‑CCT‑CCCACCCCGGTTTTCCCTGCCAGCGCCCCATTGTTGCCAAC‑CCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCt < 1:432515/151‑1 (MQ=255) aaTTTCTAAAACAGACCCCCCTTC‑CCT‑CCCACCCCCGTTTTCCCGGCCAGCGCCCCATTTTTGCCAAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttt < 1:152964/149‑1 (MQ=255) tttCTTAAACAGAACGCCCTTG‑CCT‑CCCATCCCCGTTTTGCCTGCCAGCGCGCAAATGTCGCTCAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTc < 2:760174/150‑1 (MQ=255) tttCTAAAACAGACCGCCCTTG‑CTT‑CCCATCCCGGTTTTGCCTCCCAGCGCGCCATTGTTGCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttt < 1:465428/147‑1 (MQ=255) aaacaaacCGCCTTTG‑CTT‑ACCTTCACCGCTTTCCCCGCCAGCGCCCTTATTTTCCCCAC‑CCCCTGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAttt < 1:488524/149‑1 (MQ=255) aaCAGACCGCCCTTG‑CCT‑CCCACCCCCGTTTTCCCTCCCAGCGCGCCATTGTTCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCaa < 1:12775/151‑1 (MQ=255) aCAGACCCCCCTTG‑CCT‑CCCATCCCCGTTTTGCCTGCCAGCGCGCCAATGTTGCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAg < 2:328073/151‑1 (MQ=255) gAACGCCCTTG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGAc > 2:416218/1‑150 (MQ=255) gAACGCCCTTG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGAc > 2:796915/1‑150 (MQ=255) cGCCCTTG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAgg > 1:75443/1‑151 (MQ=255) cGCCCTTG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAgg > 2:465122/1‑151 (MQ=255) cctttG‑CCT‑CCCATCCCCGTTTTGCCGCCCGGCGCCCAATTGTTCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTg < 2:949875/148‑1 (MQ=255) cccttc‑cct‑tCCACCCCCGTTTTGCCGGCCGGCGCCCAATTGTTCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa < 1:185679/146‑1 (MQ=255) cctgcccT‑CCCACCCCCGTTTTGCCTGCAGGCCCGCAATTGTTCCCTAC‑CCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa < 1:797237/140‑1 (MQ=255) ttG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCgg > 2:146803/1‑151 (MQ=255) ccT‑CCCACCCCCGGTTTCCCTGCCAGCGCGCAATTGTTCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCgg < 1:783096/148‑1 (MQ=255) cT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAAc > 1:587423/1‑151 (MQ=255) t‑tCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAAc > 2:509290/1‑150 (MQ=255) ccacccccGTTTTGCCTGCCAGCGCCCCATTGTCCCCTAC‑CCCCCGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCaa < 1:784311/141‑1 (MQ=255) ccacccccGTTTTCCCTGCCGGCGCCCAAATTTTCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGcc < 1:476106/145‑1 (MQ=255) cATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCt > 2:36486/1‑151 (MQ=255) cccccGTTTTCCCTGCCACCCCCCCATTGTTGCTTCC‑CCCCCCCTCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTg < 1:790293/147‑1 (MQ=255) ccACGGGTGTGCCTGAAAGGCCGCGATTGTTTCTTAC‑CCCCGCCGCCGCCGTGCCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTg < 2:583892/149‑1 (MQ=255) gTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt > 1:229668/1‑150 (MQ=255) ttGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAgg > 2:240285/1‑151 (MQ=255) cctgccAGCGCGCAATTGTTCCCAAC‑CCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGcc < 1:672662/150‑1 (MQ=255) ctgccAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGcc > 2:942004/1‑149 (MQ=255) ccAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAgcgc > 1:301065/1‑151 (MQ=255) cgcgcCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc > 1:72471/1‑148 (MQ=255) cgcCATTTTTCCCTCC‑CCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATc < 1:187246/151‑1 (MQ=255) tGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCtg > 2:192873/1‑151 (MQ=255) ccTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTg > 1:507322/1‑151 (MQ=255) a‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTgtc > 2:305840/1‑147 (MQ=255) | CGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTG‑CCT‑TCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCC > NC_000913/2577999‑2578291 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |