Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 869,107 | 0 | A | G | 43.9% | 55.0 / 10.9 | 41 | I530V (ATT→GTT) | gsiA | glutathione ABC transporter ATPase |
Reads supporting (aligned to +/- strand): ref base A (3/20); new base G (18/0); total (21/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.52e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAG > NC_000913/868973‑869253 | tCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTCATCGTGTgg > 1:10946/1‑150 (MQ=255) ccGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGGTTAGTCATCGTGTGGC‑gggg > 1:277489/1‑147 (MQ=255) cgacgaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGGCATCGTGTGGC‑gggg > 1:1010975/1‑146 (MQ=255) gaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGGGCGGGTTATTCATCTTGGGGC‑GGGGgtgt > 2:384614/1‑148 (MQ=255) gaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGAGGGAGGGGGGGGTAGTTTTCGTGGGGC‑gggggtttt > 2:250867/1‑142 (MQ=255) gTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGc < 2:799117/151‑1 (MQ=255) tCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCaa < 2:832585/150‑1 (MQ=255) ggCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAAtt < 2:36984/151‑1 (MQ=255) aTGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTTGTCATCGTGTGGG‑GGGGGTGT‑T‑T‑T‑TCC‑GGGGAAATTTTTTaa > 1:1053015/1‑150 (MQ=255) tttCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGGGGTAGAGTGGGGTTGGCATCCTGTGGG‑GGGGGAGT‑A‑TGT‑T‑C‑GGGGAAATTTTTTGAAtt > 1:982797/1‑150 (MQ=255) cgcgGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGgc < 2:789895/151‑1 (MQ=255) gcgGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGgcg < 1:1195/151‑1 (MQ=255) cgGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGGGGGTTAGTCATCGTGTGGG‑GGGGGGGT‑A‑T‑G‑TCC‑GGGCAAAATTTGTAAATTGTCCCAccc > 1:362417/1‑147 (MQ=255) gggCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGgcgcg < 2:416767/151‑1 (MQ=255) ggCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGgcgcg‑c < 2:32248/151‑1 (MQ=255) tCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCATCGTGTGGG‑GGGGGGGG‑A‑T‑T‑TCC‑GGGGAAATTTTTGAAATT‑GTGCCCCGCGCG‑C‑gg > 1:263637/1‑144 (MQ=255) ggtgCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGAGCG‑C‑GGTC‑GTTG‑AAAAccc < 2:919428/148‑1 (MQ=255) ttgCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCCGGTTAGTCATCGTGTGGCGGGGGATGT‑T‑T‑T‑C‑C‑GGGGCAATTTTTTAAATTGGGCCCCGCGGCG‑C‑GGCC‑TCCG‑AAAAcc > 2:1019531/1‑149 (MQ=255) tgCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGTGGGGGGGTTAGTAATCGTGTGGC‑GGGGGGGTGT‑T‑T‑T‑C‑GGGCAAATTTTTTAAATT‑GTGCCCCGCGCG‑C‑GGGG‑TTCG‑GAAAacc > 1:789053/1‑148 (MQ=255) gCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCgca < 2:697458/150‑1 (MQ=255) cGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCgcagca < 2:282145/150‑1 (MQ=255) tcCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCATCGTGTGGC‑GGGGGTGT‑A‑T‑T‑T‑CGGGGGAAATTTTTAAAATTGTTCC‑CGGCCCG‑C‑GGGCTTTCG‑AAAACCC‑CAGCAAcc > 1:88883/1‑148 (MQ=255) tcCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGGGGTAGAGGGGGTTGTTCAACGTGTGGC‑GGGGGGGT‑T‑T‑CTT‑C‑GGGGAAATGGTTTAAATTGTTCCCCGGGGCG‑G‑GGGC‑TTGG‑AAAACCC‑CCGCATcc > 2:10718/1‑148 (MQ=255) ccAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCCTCGTGTGGC‑GGGGGTGT‑T‑T‑A‑TCC‑GGGGAAATGTTTTAAATTGG‑GCACCGCGCG‑C‑CGGC‑TTTG‑gaaaaaccccgcacccct > 1:403054/1‑135 (MQ=255) aGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGGTTAGTCATCGGGTGGC‑GGGGGTGT‑A‑T‑T‑T‑CGGGGGAAATTTTTTAAATT‑TTGCCCCGCGCG‑C‑GGGC‑TTTGCAAAAcccccccaac > 1:601515/1‑137 (MQ=255) cGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCGCAGCATCCTTATAcgc < 1:802523/151‑1 (MQ=255) cGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCGCAGCATCCTTATAcgc < 2:523152/151‑1 (MQ=255) gATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTTGTCATTTTGTGGCGGGGGTTTT‑A‑T‑T‑T‑C‑GGGCAAATTGTGTGAATTGGT‑CACGGGGGGCC‑GGCC‑TTGG‑AAAAACC‑CAGCAACCTtttaatc > 2:679096/1‑145 (MQ=255) tGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCGCAGCATCCTTATACGCGTAAATTACTggc < 2:1055051/150‑1 (MQ=255) gCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCGCAGCATCCTTATACGCGTAAATTACTggcg < 1:650644/150‑1 (MQ=255) cGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTTGTCATCGTGTGGC‑GGGGGTTT‑T‑T‑T‑T‑CGGGGGAAATTTTTTAATTTGGGTC‑CCGCGCG‑C‑GGGCTTTTG‑AAAACCCCCGGAACCCTTTAACCGGTAAAttcttggg > 2:233840/1‑143 (MQ=255) gTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTcc < 1:805718/150‑1 (MQ=255) gTTTATCTCCCACGATATGGCGGTGGGAGAGCCGGTTCGTCATCGTGTTGC‑GGGCGTGG‑ATT‑C‑C‑C‑GGGCAAATTTTTGAAATTGGT‑CACCGCGCG‑CGGGCC‑TTCG‑AAAACCC‑CCCCAACCCTATACCCGTAAAATaactggggccggttc > 2:443930/1‑136 (MQ=255) tctcCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTc < 2:677724/149‑1 (MQ=255) cccACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCg < 1:278286/147‑1 (MQ=255) atGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTccc < 2:344627/151‑1 (MQ=255) tGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTccc > 2:489071/1‑150 (MQ=255) tAGAGCGGATTAGTCATCGTGTGGC‑GGCGATGG‑A‑T‑C‑C‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCGCAGCATCCTTTTACGCGTAAATTACTGGCGGCCGTTTCGGTCGCTGAACCGTCCCGACAaccacc > 2:629031/1‑151 (MQ=255) agagCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGAcc < 2:217756/150‑1 (MQ=255) gCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAg < 1:486925/151‑1 (MQ=255) ggATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGGCGCTGAACCGTCCCGACAACGACCGCAg > 1:331845/1‑149 (MQ=255) | TCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGT‑A‑T‑C‑T‑C‑GGGCAAATTGTTGAAATTGGTCCACGGCGCG‑C‑GGTC‑TTCG‑AAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAG > NC_000913/868973‑869253 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |