Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 501,564 | 0 | A | C | 30.8% | 30.8 / 15.6 | 26 | *559E (TAA→GAA) | ybaL | putative transporter YbaL |
Reads supporting (aligned to +/- strand): ref base A (16/2); new base C (0/8); total (16/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.88e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCG > NC_000913/501426‑501699 | ttAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAc > 2:493564/1‑140 (MQ=255) ttAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAc > 1:159567/1‑140 (MQ=255) tAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAcc > 1:82206/1‑140 (MQ=255) tATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGTGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTc > 2:454351/1‑140 (MQ=255) aTCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTca < 1:493564/140‑1 (MQ=255) aTCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTca < 2:733320/140‑1 (MQ=255) cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcacc > 1:108748/1‑140 (MQ=255) ttCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCa > 1:1157913/1‑140 (MQ=255) aGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTg > 1:20988/1‑140 (MQ=255) gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTc > 2:399491/1‑140 (MQ=255) tccccgccgcATCTGGTTTTGGGCTCTCGATGCCTGATGCGCCGCGGGCGCGTTTTATCATCCCTACTTTTTTTTCATTTTTTCCACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCcag < 2:1157913/136‑1 (MQ=255) cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca > 2:644368/1‑140 (MQ=255) cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca > 1:780153/1‑140 (MQ=255) aagcATCTGGCTTTTGGCTCTCAGGGCCTGATGCGCCGCGGGCGCGTTTTATCATCCCTACATATTTTTCATATTTTCCCCCCGGAACCCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTc < 2:732446/138‑1 (MQ=255) cttttGGCTTTTGCCCCTCGTTCCTTGATGGGGCGCGGGCGCGTTTTATCATCCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAg < 2:219864/136‑1 (MQ=255) aTCTGGCATTTGGCTCTCGATGCCTGATGCCACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGc > 2:15919/1‑140 (MQ=255) tGGCTCTCGGTGCCTGATGCGACGCGGGCGCGTTTTATCATCCCTACATATTTTTCATTTTTTACATCCGGCACCCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg < 2:780153/140‑1 (MQ=255) tGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTAcctg > 1:732446/1‑137 (MQ=255) tGGCTCTCGATGCCTGATGCGACGCGGGCGCGTCTTATCATGCCTACATTTTTTTCATTTTTTACACCCGGCACCCACCGTTTCCCCCGTCCCCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg < 2:1106784/140‑1 (MQ=38) tGATGCGCCGCTGGCGCTTCTTATCATGCCTACATTTTTTTCATTTTTTCCACCCGGCACCCACCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCg < 2:849407/140‑1 (MQ=255) gATGCGCCGCTGGCGCGTTTTATCATGCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc < 2:297620/140‑1 (MQ=255) gCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCa > 1:224314/1‑140 (MQ=255) ttCTTATCATCCCTACATTTTTTTCATTTTTTCCACCCGGCACCCACCGTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt < 2:82206/139‑1 (MQ=255) cTACATTTTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTc < 1:454351/140‑1 (MQ=255) tACATTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg < 1:399491/140‑1 (MQ=255) aCATTTTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg < 2:159567/140‑1 (MQ=255) tttttCATATTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGata < 2:20988/140‑1 (MQ=255) ttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc > 2:474405/1‑140 (MQ=255) ttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc > 1:317671/1‑140 (MQ=255) tACATCCGGCACCCCCCGTTTCCCCCCCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa < 2:224314/140‑1 (MQ=255) ctcccGGCAACCACTTTTCCCCCCGTCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATc < 1:644368/137‑1 (MQ=255) cccGGCAACCACCGTTTCCCCCGTCCCCACCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGt < 1:561746/139‑1 (MQ=255) ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt > 2:872060/1‑140 (MQ=255) ttttCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATTGGcgcg < 2:676471/139‑1 (MQ=255) | TTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCG > NC_000913/501426‑501699 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |