Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,679,688 | 0 | T | G | 28.2% | 64.1 / 10.1 | 39 | G44G (GGT→GGG) | ydgI | putative arginine:ornithine antiporter |
Reads supporting (aligned to +/- strand): ref base T (9/19); new base G (11/0); total (20/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.45e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATGG > NC_000913/1679552‑1679823 | tatCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTAGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGGAtt > 2:504956/1‑140 (MQ=255) atCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTa > 1:721042/1‑140 (MQ=255) gAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGgcgc > 2:801461/1‑140 (MQ=255) gAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGgcgc < 1:908234/140‑1 (MQ=255) aCTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCAttt > 1:70917/1‑140 (MQ=255) aCTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCAttt < 1:801461/140‑1 (MQ=255) gggACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGGATTACTGGCGCTGGCATTTTAt > 1:573889/1‑140 (MQ=255) gACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTg < 2:721042/140‑1 (MQ=255) cTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGGTGGGGGTTTACTTGCCCTGGCAATTTATTGCt > 2:847668/1‑140 (MQ=255) aGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGc < 1:436816/140‑1 (MQ=255) aGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGc > 1:798563/1‑140 (MQ=255) cattCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCCTCGGCTGGGGGTTTATTGGCGCTGGCATTTTAATTCTGGGCttt > 2:1006537/4‑140 (MQ=255) aCTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTg < 1:1006537/140‑1 (MQ=255) ccGCGCTGGGATGAAGCTCAATGCTGGTCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGTCGTTACCAGCCCTCCAGCACTGCTCATCCGCTGGGGAACTACTGGCGCTGCCGCTTTTTTGGTTGGCTTTGCCATg > 2:32269/1‑140 (MQ=255) cgcTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTg < 1:686258/140‑1 (MQ=255) cgcTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGGATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTg > 2:664991/1‑140 (MQ=255) cAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGGATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCAtt > 2:910311/1‑140 (MQ=255) aaTGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTc > 1:460625/1‑140 (MQ=255) aaTGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGGATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTc > 1:707378/1‑140 (MQ=255) aaTGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGGATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTc > 1:276286/1‑140 (MQ=255) tGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGt > 1:715947/1‑140 (MQ=255) ggCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGGATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGa > 2:236321/1‑140 (MQ=255) gggTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTTTTGCTGGCCTTTGCCCTGCTGATCCTCACGCGCATTCGTCCCGAACttt > 2:1035316/1‑139 (MQ=255) ggTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGGATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGa > 2:876641/1‑140 (MQ=255) tttCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAcggcg < 1:1140369/140‑1 (MQ=255) tttCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAcggcg < 2:798563/140‑1 (MQ=255) tttCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAcggcg < 2:707378/140‑1 (MQ=255) tCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGt < 2:573889/140‑1 (MQ=255) gTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGGATTACTGGCGCTGGCATTTTATTGCTGGGCTTTTCCCTGCTGATCCTCACGCGCCTTCGTCCCGAACTTGACGGCGGGATc > 2:499420/1‑140 (MQ=255) gCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTa > 1:2136/1‑140 (MQ=255) aaaaTATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGGATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGGATCTTTACCTAt > 1:328900/1‑140 (MQ=255) ggcggcAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCcgcg < 1:876641/140‑1 (MQ=255) gcggcAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGa < 2:682542/140‑1 (MQ=255) gcggcAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGa < 1:910311/140‑1 (MQ=255) cggcAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGaa < 1:217191/140‑1 (MQ=255) ggcAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGGATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGGTCCTCACGCGCATTCGTCCCGAACTTGACGGCGGGATCTTTTCCTATGGCCGCGAAg > 2:427685/1‑140 (MQ=255) gCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTgg < 1:1035316/140‑1 (MQ=255) ccAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGc < 2:1068346/140‑1 (MQ=255) cagcaCTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATc < 1:414115/140‑1 (MQ=255) cTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTg < 1:847668/140‑1 (MQ=255) cTCATCGGCTGGGGGATTACTGGCGCCGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTg > 2:651943/1‑140 (MQ=255) tGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATg < 1:499420/140‑1 (MQ=255) ggggTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATgg > 2:127240/1‑140 (MQ=255) | TATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATGG > NC_000913/1679552‑1679823 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |