| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,222,233 | 0 | G | A | 45.7% | 30.9 / 16.4 | 35 | pseudogene (60/75 nt) | ycgH | pseudogene; putative ATP‑binding component of a transport system |
| Reads supporting (aligned to +/- strand): ref base G (12/7); new base A (0/16); total (12/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.26e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAAT‑G‑TGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAC > NC_000913/1222085‑1222379 | gtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAAT‑G‑TGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑a > 1:235825/1‑150 (MQ=255) cGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAAT‑G‑TGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGt < 2:780824/151‑1 (MQ=255) gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAAT‑G‑TGGGCGTGCCGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGtt > 1:449342/1‑151 (MQ=255) gATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAAT‑G‑TGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATc < 2:230989/151‑1 (MQ=255) aTGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAAT‑G‑TGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATcc < 2:736334/151‑1 (MQ=255) aGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAAT‑G‑TGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGtt > 1:653754/1‑150 (MQ=255) gaaaTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAAT‑G‑TGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGtt > 1:421515/3‑149 (MQ=255) aaTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAAT‑G‑TGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCt > 2:449028/1‑151 (MQ=255) aTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAAT‑G‑TGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCt > 2:512956/1‑150 (MQ=255) ggaaaaTTAAACTTGGCGTAAATGGCACCCTAAATCCAGGGGCAGGCGTGGGGGGAAT‑G‑TGGGCGTGCAGCTGGGTGATAATGGCTACAATGACCCCGCAGTGTGGGGGGGCTTG‑AAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAAc < 2:43269/148‑1 (MQ=255) ggaaaaTTAAACTGGGCGAAAATGCAACCCTGAACCCGGCGGCCAGTTTCTGGGGAAT‑G‑GGGGGGGGCAGTTGGGTGAAAATGGATACAATGACCCCGAAGGAATGGTGGGCCTA‑AAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAAc < 1:377961/148‑1 (MQ=255) aaaTTACATTGGGCGTAACTGGCAACCTAACCCCCGCGCCCGGCGTCGGGGAAAT‑G‑GGGGCGTGCGGCTGGTTAAGAAGGACTACAATACCCCCGCAGTGATGGGGGCTCGG‑AAATAAAACTTTTGCCCCCGCCGTTAGCTAAAAAAcc < 2:280651/149‑1 (MQ=255) aaaCTTGGCGAAAATGCCACCCTGAACCCAGCGGCCAGCGTCTGGGGAAT‑G‑GGGGCGTGCAGCGGGGTGATAATGCCTACAATGACCCCGCGGTGAGGGGGGGCCTA‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg < 1:643708/150‑1 (MQ=255) aaaCTGGGAGAAAAGGGCACCCTGAATCCAGCGGCCGGGGTGGGGGGATG‑G‑GGGGGGGGCGGTTGGGTAATAATGGCAACAAGGACACCGCAGTGATGGGGGGCTTA‑AAATAAAATTTCTGACCCCGCCGTTAGCTAAAAAACCGcgtcg < 1:87964/150‑1 (MQ=255) ggCGAAAATGGCAACCTGAATCCAGCGGCCGGCGTCGGGGGAAT‑G‑TGGGCGTGCAGCGGGGTGATAATGGCTACAATGACCCCGCAGTGATGGTGGGCTTA‑AAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCa < 2:65129/150‑1 (MQ=255) aagggCAACTTGAAACCAGCGCCCGGCGTGGGGGGATT‑G‑GGGGGTGGCAGTGGGGTAAAAGGGGCAAAAATGCCACCG‑AGGTAAGGTGGGGCAGAAAATATAGGTTCGGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc < 1:520678/148‑1 (MQ=255) tGAATCCAGCGGCCAGCGTCTGGGGAAT‑G‑TGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg > 1:523960/1‑151 (MQ=255) aacccggCGGCCAGCGTCGGGGGAAT‑G‑GGGGCGTGCAGCTGGGTAATAATGGCTACAATGACCCCGCAGTGTTGGTGGGCCTA‑AAATATAGGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa < 2:534228/145‑1 (MQ=255) cggCGGCCGGCGTCT‑GGGAAT‑GTTGGGTGTCCAGCTGGGTGATAATGTCTAAAATAACCCCTCTTTGATGGGGGGCCTA‑AAATATACTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAatat < 2:372283/149‑1 (MQ=255) gCGGCCAGCGTCTGGGGAAT‑G‑TGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt > 2:50794/1‑150 (MQ=255) ggcaggCGTGTGGGGAAT‑G‑TGGGCGTGAAGCGGGGTGATAATGGCTAAAATAACACGGCAGTGAGGGTGGGCTTA‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAatat < 2:879185/142‑1 (MQ=255) ggCCGGCTTCGGGGAAGT‑G‑GGGGCGGGAATCGGGATGATAATGGCTAAAATCCCCCCGAAGTGGGTGTGGGCCAA‑AAATTTAATTCCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata < 1:497069/146‑1 (MQ=255) ccAGCGTCGGGGGATG‑G‑TGGGGGTGCAGCTGGGTGATAGGGGCTAAAATGACACCGCAGTGAGGGTGGGCTTA‑AAAAATAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt < 2:558039/146‑1 (MQ=255) cggcgtgggggAATGG‑TTGCCGTGCAGCGGGTTAATATTGCTTGCAACGAATCCACGGTGGTGGTGTGATCG‑AAATATAAGTTGTTTCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatta < 1:199052/139‑1 (MQ=255) ggggATT‑G‑TGGGCTTGCAGCGGGTTAATAATGGCTAAAATGCCCCGGCATTGATGGTGGCCAAA‑AAATTTAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa < 1:433507/151‑1 (MQ=255) ggggatgggggggGTGCAGCTGGGAGAAAATGACTAAATAAACCCCAAGGGGATGGGGGTCTTA‑AAATATAATTTTTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa < 2:649709/142‑1 (MQ=255) aaggggggCGTGCAGCTGGTTAATAATGTCAAAAATGCCACCGGGGTAGTGGGGGGCTTA‑AAATATATGTTCTTATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCa < 1:434381/145‑1 (MQ=255) tGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATc > 1:866035/1‑151 (MQ=255) aCAATGACCCGGCAGTGATGGTGGGCAAA‑AAATTTAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGcc < 1:834425/140‑1 (MQ=255) aaagACACCGCAGGGAGGGTGGCCCTA‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc < 2:314026/147‑1 (MQ=255) gCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt > 2:522989/1‑151 (MQ=255) gCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt > 2:122410/1‑151 (MQ=255) ggggggCTTA‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTa < 2:217022/145‑1 (MQ=255) ggggggCGTA‑AAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcg < 2:743600/148‑1 (MQ=255) gCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAc > 2:501015/1‑151 (MQ=255) ctag‑AAATAAAGGTTCGAAGCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCaaa < 1:456900/146‑1 (MQ=255) | GTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAAT‑G‑TGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAC > NC_000913/1222085‑1222379 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |