Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,222,233 | 0 | G | A | 46.4% | 25.9 / 12.4 | 29 | pseudogene (60/75 nt) | ycgH | pseudogene; putative ATP‑binding component of a transport system |
Reads supporting (aligned to +/- strand): ref base G (13/2); new base A (0/13); total (13/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.83e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGG‑CCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG > NC_000913/1222088‑1222382 | aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGG‑CCTGAAat > 2:738199/1‑150 (MQ=255) aGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGG‑CCTGAAATATAAGTTCTGAt < 2:431116/151‑1 (MQ=255) aGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGG‑CCTGAAATATAAGTTCTGAt > 1:76864/1‑151 (MQ=255) gAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGG‑CCTGAAATATAAGTTCTGATCCCGCCGtt > 1:344123/1‑151 (MQ=255) aaTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGG‑CCTGAAATATAAGTTCTGATCCCGCCGTTAGCt > 1:790521/1‑151 (MQ=255) ttgggggaAAATAAAATTTGGCGTAAATGGCAACCTGAATCCAGCGCCCGGCGTCTGGGAAATGTGGGCGTGCAGCTGGGTGATAGTGGCTACAATGACCCCGAAGTGGTGGTGGG‑CATGAAATATAATTTCTGATCCCGCCGTTAGCTa < 1:149571/147‑1 (MQ=255) aaCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGG‑CCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcgt > 1:379684/1‑150 (MQ=255) aCTTGGAGAAAAGGGCACCCGAAACCCAGGGCCCAGCGTCGGGGGAAGGGGGGCGGGCAGCGGGGTAATAATGTCTACATTGACCCGGCAGGGTTGGTGGG‑CCTAAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAt < 2:295221/151‑1 (MQ=255) ggggaaaaGGGCACCCGAAACCCGGCGCCCGGCTTGTGGGGAATGGGGGGGTGAAGCGGGGTAATAAGGGCAAAAAGGCCCCCGAGGTGATGGGGGG‑CTTAAAATATAAGTTCTAACCCCGCCGTTAGCTAAAAAACCGCGTCGTa < 1:463910/141‑1 (MQ=255) aactagacccccgCGGCATGCGTGTGGGGAGGGTGGGGGGGAGGGGGGGTAATAATGGAAAAAATGACCCCGCTGTGGTGGTGGG‑CCAAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgcg < 1:433640/134‑1 (MQ=255) aaTCCGGCGGCCGGCGTCGGGGGAATGGGGGCGTGCGGCTGGGTGATAAGGGCTAAAATGCCCCCGGGGTGAGGGGGGG‑CCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa < 2:153296/151‑1 (MQ=255) ccAGCGGTCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGG‑CCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat > 2:223803/1‑150 (MQ=255) ggttttGGGGAGTGGGGGCGTGCTGTAGGATGAAAGTGGCAACAAGAACTCCGGGGT‑ATGGGGGGACCGTAAAAATAACTTCTTTGCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTg < 1:127075/146‑1 (MQ=255) cGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGG‑CCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc > 2:447937/1‑151 (MQ=255) ggggatggggggcgggcGGTTGGGTAATATGGTTTACAATGCCCCCGCGGTGGTGGTGGG‑CATAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa < 1:693778/142‑1 (MQ=255) ggggATTGGGGGCGTGAAGCTGGGTAATAATGGCTACAATGCCCCCGCAGGGTTGGGGGG‑CCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa < 2:374896/151‑1 (MQ=255) ggggAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGG‑CCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAt > 2:345604/1‑150 (MQ=255) gggAGTGGGGGCGTGCGGTGGGGTGATAAGGGCTAAAATGACCCCGCGGGGATGGTGGC‑CCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAg < 2:76863/151‑1 (MQ=255) gtgGGGGTGCTGTGGGTTAATAATGGCAAAAAAGCCCCGGCGGTGTTGGGGGG‑CTTAAAATATATTTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGc < 1:280814/146‑1 (MQ=255) cTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGG‑CCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcat > 1:558405/1‑148 (MQ=255) gggTAATATTGGCTACAATGCCCCCGAAGTGAGGGTGGG‑CTTAAAATAAAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAt < 1:187671/150‑1 (MQ=255) aaTGGCTACAATGCCCCCGCGGTGGTGGGGGG‑CCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg < 2:233061/150‑1 (MQ=255) tGACACCGCAGTGATGGTGGG‑CCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAtt > 1:190141/1‑151 (MQ=255) cacCGCAGTGATGGTGGG‑CCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGc > 2:139511/1‑151 (MQ=255) ccGCAGTGATGGGGGG‑CCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCg < 2:101817/150‑1 (MQ=255) gCAGTGATGGTGGG‑CCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt > 1:493821/1‑151 (MQ=255) gTGAGGGGGGG‑CCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt < 2:757010/150‑1 (MQ=255) gtggT‑CCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt < 1:368655/145‑1 (MQ=255) cTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCt > 2:763482/1‑151 (MQ=255) gAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTg > 2:3751/1‑150 (MQ=255) | AAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGG‑CCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG > NC_000913/1222088‑1222382 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |