Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 993,917 | 0 | A | C | 63.0% | 14.3 / 28.5 | 27 | V43G (GTG→GGG) | ssuB | aliphatic sulfonate ABC transporter ATPase |
Reads supporting (aligned to +/- strand): ref base A (9/1); new base C (0/17); total (9/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.13e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-05 |
GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG > NC_000913/993768‑994064 | gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCcca > 1:476802/1‑150 (MQ=255) cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCcaccac < 2:80263/150‑1 (MQ=255) ttgtCTTCCTGAATTCCAGCCGGCGGTGTGGTCCCCGCTAACCCACCGCCTGCGGTTGGCTTTCCCAGCCCTGCCAGCAGGCGCAGCAGGGTTCTTTTCCCACCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTa < 2:378490/150‑1 (MQ=255) ctcccTGAATTCCGCCCGGGGGTGGGGTCCCCCCTAACACACCCCCTGGGGTGGGCTTTCCCACCCCTGCCAGCGGGCGCAGCAGGGTTCTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATc < 1:10672/148‑1 (MQ=255) cTTCCTAAATTTCACCCGGGGGTGTGGTGCCCGCTAACACATCGCCTGCGTTGGGCGTTCCCAGCCCTGCCAGCAGGCGCAGCAGGGTACTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATc < 2:37290/151‑1 (MQ=255) ttCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGAGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGACACCACCGCTGCGGCCCACCACCGACACAAACTGACCTTCAGGAATCTGTAAATc > 1:17228/1‑150 (MQ=255) gcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat > 2:16819/1‑151 (MQ=255) ggtgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTc > 1:216664/1‑150 (MQ=255) ggcccccGCAAAACCACCGCCGGCGGTTGGCTTTCCCAGCCCTGCCGGCGGGCGAAGCGGGGTCTTTTCCCCCCCCCCGCTGGGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTc < 1:588873/146‑1 (MQ=255) aaaCACCCCGCGTGCGGTTGGCGTTCCCAGCCCTCCCAGAGGCCGCAGCAGGGTCCTTTCCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt < 1:399391/150‑1 (MQ=255) cggCGTTTGGCTTTCCCAGCCCTGCCAGCGGGCGCAGCGGGTTACTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc < 2:728586/149‑1 (MQ=255) cGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAtt > 2:638637/1‑151 (MQ=255) tttCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa > 1:597328/1‑151 (MQ=255) tttCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa > 2:392339/1‑151 (MQ=255) ttcccaccccTGCCAGCAGGCGCAGCAGGGTATTTTGCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa < 2:388405/143‑1 (MQ=255) ttcccAGCCCTGCCAGCAGGCGCAGCAGGGTCCTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa < 2:718361/148‑1 (MQ=255) ttcccAGCCCTCCCAGCAGGCGCAGGGGGGTTCTTTTCCCCCCCCCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa < 2:140795/148‑1 (MQ=255) cccTGCCAGCAGGCGCAGCAGGTTATTTTTGCCCCCACCGCTGCGCCCCCCCACCGCCGCAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg < 2:288953/148‑1 (MQ=255) ccTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg > 2:422583/1‑150 (MQ=255) gCCAGCGGGCGCAGCGGGTTACTTTGCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGccc < 2:309485/150‑1 (MQ=255) agcagGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt > 1:691566/1‑151 (MQ=255) cagcgggttctttttGCCCCCCCCGCTGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg < 1:298299/139‑1 (MQ=255) cagcagGTTACTTTCCCCACCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg < 1:583382/150‑1 (MQ=255) gTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt > 1:588555/1‑151 (MQ=255) tgcccccccccGCTGCGGCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc < 1:118755/142‑1 (MQ=255) tccccccccccGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc < 1:148029/142‑1 (MQ=255) tccccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc < 1:8597/142‑1 (MQ=255) cccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg < 2:566640/144‑1 (MQ=255) cccaccaCCGGTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg < 1:584723/150‑1 (MQ=255) ccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATg < 1:70123/151‑1 (MQ=255) | GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG > NC_000913/993768‑994064 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |