Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AM260479 | 3,418,719 | 0 | C | . | 27.8% | 41.8 / 13.6 | 19 | intergenic (+258/‑317) | h16_A3158/corA2 | predicted permease of the drug/metabolite transporter (DMT) superfamily/Putative Mg2+ and Co2+ transporter |
Reads supporting (aligned to +/- strand): ref base C (8/5); new base . (2/3); total (11/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.08e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGTTGCCGCAGCAGCGCTGCGACCTCTGCGAATGCCATGTCACTGCGGGCTGCGTAACGACGACCCGCCCGGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAGCAGAAC‑ACTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGCCATCCATGCGCGCCGCATCGCGTGCGGATCTTCCCTCTCCCGGCTAGCCAGCCAACCCAGGCAAACGCCATCGC > AM260479/3418610‑3418835 | aGTTGCCGCAGCAGCGCTGCGACCTCTGCGAATGCCATGTCACTGCGGGCTGCGTAACGACGACCCGCCCGGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAGCAGAAC‑ACTGCGGGCGAGTGTCttt > 2:429037/1‑134 (MQ=255) agcagcGCTGCGACCTCTGCGAATGCCATGTCACTGCGGGCTGCGTAACGACGACCCGCCCGGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAGCAGAAC‑ACTGCGGGCGAGTGTCTTTCCGTGCCAc < 1:1107800/134‑1 (MQ=255) gCGACCTCTGCGAATGCCATGTCACTGCGGGCTGCGTAACGACGACCCGCCCGGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAGCAGAAC‑ACTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGgcgc > 1:640266/1‑134 (MQ=255) cACTGCGGGCTGCGTAACGACGACCCGCCCGGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAGCAGAAC‑ACTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGCCATCCATGCGCGCCGCATcgcg > 2:908538/1‑134 (MQ=255) cTGCGGGCTGCGTAACGACGACCCGCCCGGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAGTGTTCT‑GCTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGccat > 1:208979/1‑113 (MQ=37) cTGCGGGCTGCGTAACGACGACCCGCCCGGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAG‑TGTTCTGCTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGccat > 1:517960/1‑113 (MQ=37) cTGCGGGCTGCGTAACGACGACCCGCCCGGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAG‑TGTTCTGCTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGccat < 2:208979/113‑1 (MQ=37) cTGCGGGCTGCGTAACGACGACCCGCCCGGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAG‑TGTTCTGCTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGccat < 2:517960/113‑1 (MQ=38) ggCTGCGTAACGACGACCCGCCCGGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAGCAGAAC‑ACTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGCCATCCATGCGcgccg < 1:343411/120‑1 (MQ=255) ggCTGCGTAACGACGACCCGCCCGGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAGCAGAAC‑ACTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGCCATCCATGCGcgccg > 2:343411/1‑120 (MQ=255) acgacgACCCGCCCGGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAGCAGAAC‑ACTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGCCATCCATGCGCGCCGCATCGCGTGCGGATCTTCCctct > 2:64544/1‑134 (MQ=255) gacCCGCCCGGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAGCAGAAC‑ACTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGCCATCCATGCGCGCCGCATCGCGTGCGGATCTTCCCTCTCCCgg > 2:873602/1‑134 (MQ=255) cgcccgGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAGCAGAAC‑ACTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGCCATCCATGCGCGCCGCATCGCGTGCGGATCTTCCCTCTCCCGGCTag < 2:640266/134‑1 (MQ=255) gcccgGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAGCAGAAC‑ACTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGCCATCCATGCGCGCCGCATCGCGTGCGGATCTTCCCTCTCCCGGCTagc > 2:620397/1‑134 (MQ=255) gcccgGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAG‑TGTTCTGCTGCGGGCGAGTGTCTTTCCGTGCCa > 1:1064534/1‑76 (MQ=18) gcccgGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAG‑TGTTCTGCTGCGGGCGAGTGTCTTTCCGTGCCa < 2:1064534/76‑1 (MQ=18) cTCCGCGCCTCGTCGGCCTCGCCCGCAGCAGAAC‑ACTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGCCATCCATGCGCGCCGCATCGCGTGCGGATCTTCCCTCTCCCGGCTAGCCAGCCAACCCAGGc > 1:230378/1‑133 (MQ=255) gcCTCGTCGGCCTCGCCCGCAGCAGAAC‑ACTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGCCATCCATGCGCGCCGCATCGCGTGCGGATCTTCCCTCTCCCGGCTAGCCAGCCAACCCAGGCAAACGcc < 2:1070394/134‑1 (MQ=255) gtcgGCCTCGCCCGCAGCAGAAC‑ACTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGCCATCCATGCGCGCCGCATCGCGTGCGGATCTTCCCTCTCCCGGCTAGCCAGCCAACCCAGGCAAACGCCATCGc < 2:230378/134‑1 (MQ=255) cagcagAAC‑ACTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGCCATCCATGCGCGCCGCATcg < 2:278284/66‑1 (MQ=255) cagcagAAC‑ACTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGCCATCCATGCGCGCCGCATcg > 1:278284/1‑66 (MQ=255) | AGTTGCCGCAGCAGCGCTGCGACCTCTGCGAATGCCATGTCACTGCGGGCTGCGTAACGACGACCCGCCCGGAGCACCCGGCTCCGCGCCTCGTCGGCCTCGCCCGCAGCAGAAC‑ACTGCGGGCGAGTGTCTTTCCGTGCCACCGGCGGCGCCATCCATGCGCGCCGCATCGCGTGCGGATCTTCCCTCTCCCGGCTAGCCAGCCAACCCAGGCAAACGCCATCGC > AM260479/3418610‑3418835 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |