Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,630,340 | 0 | A | C | 44.4% | 57.1 / 13.4 | 46 | G61G (GGT→GGG) | ettA | energy‑dependent translational throttle A |
Reads supporting (aligned to +/- strand): ref base A (15/10); new base C (0/20); total (15/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.23e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTT > NC_000913/4630196‑4630485 | cgcATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑g > 2:751607/1‑151 (MQ=255) aTACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATg < 2:96792/151‑1 (MQ=255) caaGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAAt < 1:835631/149‑1 (MQ=255) cAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATg < 2:184061/151‑1 (MQ=255) gCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGccc > 1:319833/1‑151 (MQ=255) ttCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGcc < 1:329416/150‑1 (MQ=255) ttCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGcc < 1:463842/150‑1 (MQ=255) ttCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGcc > 2:162483/1‑150 (MQ=255) aGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGCCATa > 2:579518/1‑150 (MQ=255) gggCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGCCATAAt > 2:611819/1‑151 (MQ=255) gggCGTAACCCACTTCTGAACCCGTTTCTTCAAGGGACCCCCGCACGGTGTTTCCCGGGTTCAGCTGCGGTTCCTGCGGCAAAAACCCAT‑TCTTGATGTCTGCCGGCGGACGCGCTTCCCCTTG‑GATGTTTTTATCAAGGCCCGCCATa < 1:421474/149‑1 (MQ=255) acccactctggAAACCGTTTTTTAAAGGGACTCACGCACGGTGTGTTCGGGTTTCAGCTGCGGTTCCTGCGGCAGATACCCAA‑TCTTGATGTCGGGCTGCGGCCGCGCTTCCCCTTC‑GATGTCTTTATCAATGCCCGCCATAATGCgcag < 2:687906/140‑1 (MQ=255) ccACTTCTAAAACCGTTTCTTCAATGGACCCACGCACGGTGTGTCCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCCCCTTC‑GATGTTTTTACCAATCCCCGCCATAATGCGcagcag < 2:847532/151‑1 (MQ=255) cctccTGAACCCGCTTCTTCAAGGGCCCCACGCCGGGTGTGTTCGGGGTTCAGCGGCGTTCCCGGCGGAAAATAACAAT‑TTTTGTTGTCTGGCGGCGGACGCGCTTC‑CCTTCGGATGTTTTTATCATTCCCCGCCATAATGCGCAGCAgg < 1:337915/146‑1 (MQ=255) aCTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGCCATAATGCGCAGCAggg > 2:196997/1‑151 (MQ=255) ctaaacccGCTTCTTCAATGGCCTCCCGCCCGGTTTGTTCCGGTTTCAGCTGGGGTCCCTGCGGCAGATACCCAA‑TTTTGTTGTCGGGCTGCGGACGCGTTTCCCTTTC‑GATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTg < 2:626228/143‑1 (MQ=255) acccGTTTCTTAAGTGGACTCGCGCGCGGTGTTTTCCGGTTCCGGCGGCGTTTCCTGGGGCA‑AAAACCAATTTTTTATTTCTGGTTGGGGACGCCGTTCACCTTG‑GTTGTTTTTAAAAATGCCCGCCATAATGCGCAGCAGGGTGGAc < 2:669625/146‑1 (MQ=255) ccGCTTCTTCAATGGCCCCTCCCGCGGTGTGTTTCGGTTTCAGCTGCGGTCCCTGCGTCATATAGCCAT‑TCTTGATGCCTGACTCCGCAGCCGCTTCCCCCTC‑GATGTCTTTAGCAATGCCCGCCATAATGCGCAGCAGGGTGGACtt < 2:675661/148‑1 (MQ=255) ccGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTAcc > 2:664194/1‑151 (MQ=255) aTGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTc > 1:570600/1‑150 (MQ=255) gACTCACCCACGGTGTTTCCCGGGTTCAGCTGGGGTTCCTGCGCCAGATAACCAT‑TCTTGATGTCTGGCTGCGGGCGCGCTTCCCCTTC‑GTTGTTTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGAc < 2:194244/151‑1 (MQ=255) cccgcacgGTGTGTCCCGGTTTCAGCTGCGTTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATTTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCa < 2:396456/148‑1 (MQ=255) cgccgggTGTGTCCGGGGTCCTGCTGCGGTTCCTGCGGCAGAAAACCAA‑TCTTGATGTCTGGCGGCGGCCGCGCTCCACCTTG‑GATTTCTTTACCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAg < 1:419606/144‑1 (MQ=255) cgcgGTGTTTTCCGGGTTCTGCTGCTGTTCCGGGGGCAGATACCCAT‑TCTTGATGTCGGGCTGGGGAGGCGCTTCCCCTTC‑GATGTCTTTATAAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGaca < 1:777049/149‑1 (MQ=255) cacgGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAT‑TCTTGATGTCTGGCTGCGGACGCGCTTCCCCTTC‑GATTTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGaca < 2:385359/151‑1 (MQ=255) gggTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGccc > 2:758/1‑151 (MQ=255) tCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAgg > 2:507790/1‑151 (MQ=255) gCGGTCCCTGCGAAAAAAAACCAA‑TTTTGATGTTTGGCTGGGGGCGCGCTTCCCCTTC‑GTTGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGga < 1:197934/147‑1 (MQ=255) gtcccTGCGGCAGAAAACCAT‑TTTTGTTTTCGGGCTGGGGCCGCGTTTCCCTTTG‑GTTTTCTTCACAAATCCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAggg < 1:512798/140‑1 (MQ=255) gTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACt > 2:739527/1‑151 (MQ=255) ccTGCGGCAGATAACCAA‑TTTTGATGTCTGGCTGCGGCCGCGCTTCCCCTTC‑GATGTTTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCa < 2:876943/150‑1 (MQ=255) ccTGCGGCAGATAACCAA‑TCTTGATGTCTGGCGGCGGACGCGCTTCCCCTTC‑GATTTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCa < 2:405693/150‑1 (MQ=255) gCGGCAGATAACCAA‑TTTTGATGTCTGGCTGCGGCCGCGCTTCCCCTTC‑GATGTTTTTATCAATCCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCagaga < 2:505777/151‑1 (MQ=255) ggCAGATAACCAA‑TCTTGATGTCTGGCTGCGGCCGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAgagaga < 2:180979/151‑1 (MQ=255) ggCAGATAACCAA‑TCTTGATGTCTGGCGGCGGCCGCGCTTCCCCTTC‑GATGTTTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAgagaga < 2:507821/151‑1 (MQ=255) aGAAAACCAA‑TCTTGATTTCTGGCTGCGGCCGCCTTCCCCCTTG‑GATGTCTTTATAATTGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATg < 2:398856/150‑1 (MQ=255) aaCCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGttttt > 2:254691/1‑151 (MQ=255) tGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATg > 2:747210/1‑151 (MQ=255) tttCTGGCTGGGGACGCCCTCCCCCTTC‑GATTTTTTTATCAATCCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGa < 1:25156/147‑1 (MQ=255) tGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACg > 1:534304/1‑151 (MQ=255) gTCGGGCTGCGGACGCGCTTCCCCTTC‑GATGTTTTTAACAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACg < 1:751017/150‑1 (MQ=255) ggcggcGGACGCGTTTCACCTTC‑GATTTCTTTATAAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTcg < 1:873460/151‑1 (MQ=255) ggCTGCGGACGCGCTTCCCCTTC‑GATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAatat < 2:41371/142‑1 (MQ=255) tGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTcggcg > 1:27620/1‑151 (MQ=255) ggccgctCTTCCCCTGG‑GTTGTCTTTACCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGtt < 1:583883/139‑1 (MQ=255) ttCCCCTTC‑GATGTTTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACtt < 2:546390/149‑1 (MQ=255) | CGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAA‑TCTTGATGTCTGGCTGCGGACGCGCTTCACCTTC‑GATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTT > NC_000913/4630196‑4630485 |
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: ‑ |