| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,590,331 | 0 | C | G | 22.7% | 39.0 / 11.1 | 22 | D317H (GAC→CAC) | yjiY | putative transporter |
| Reads supporting (aligned to +/- strand): ref base C (8/9); new base G (2/3); total (10/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AACGCGCGTCGGTTTCATTAGCCAGCAGTTTTGGCGTCGTACCGGAAGAGATCAGCGCGTGGAAGCCAGATACCGCACCACAGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTCGCGCGGTGCGAGGATCAGCCACACTGGCAGCAGTGC > NC_000913/4590190‑4590472 | gcactgctgccagtgtggctgatcctcgcaccgcgcgactatctggcaaccttcctgaaaatcggcgttatcgtcggcctggcgctgggtatcgtggtgctgaacccggaactgaaaatgcctgccatgacCCAGTACCGTCAATGTAc < 8:21264‑M1/18‑1 (MQ=255) ctgctgccagtgtggctgatcctcgcaccgcgcgactatctggcaaccttcctgaaaatcggcgttatcgtcggcctggcgctgggtatcgtggtgctgaacccggaactgaaaatgcctgccatgacCCAGTACCGTCAATGTACTgg > 7:21264‑M1/129‑149 (MQ=255) cagcagTTTTGGCGTCGTACCGGAAGAGATCAGCGCGTGGAAGCCAGATACCGCACCACAGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCa < 2:6560/149‑1 (MQ=255) ggAAGAGATCAGCGCGTGGAAGCCAGATACCGCACCACAGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGa > 6:111243/1‑149 (MQ=255) aGCGCGTGGAAGCCAGATACCGCACCACAGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTc < 4:85492/129‑1 (MQ=255) gCCAGATACCGCACCACAGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTcc < 3:112908/112‑1 (MQ=255) gCCAGATACCGCACCACAGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTcc > 4:112908/1‑112 (MQ=255) ccacAGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGAtttt > 1:123822/1‑149 (MQ=255) aGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGc < 3:103835/138‑1 (MQ=255) aGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGc > 4:103835/1‑138 (MQ=255) cGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAg < 5:111243/149‑1 (MQ=255) cGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTcgcgcg < 2:123822/149‑1 (MQ=255) ccggaactgaaaatgcctgccatgacCCAGTACATTGACGGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTc < 2:29829‑M1/116‑1 (MQ=255) ccggaactgaaaatgcctgccatgacCCAGTACATTGACGGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTc > 1:29829‑M1/27‑142 (MQ=255) atgcctgccatgacCCAGTACATTGACGGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAgc > 2:170770‑M1/15‑83 (MQ=255) atgcctgccatgacCCAGTACATTGACGGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAgc < 1:170770‑M1/69‑1 (MQ=255) catgacCCAGTACATTGACGGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTCGCGCGGTGCGAGGATCAGCCACACTgg < 2:125752‑M1/143‑1 (MQ=255) aGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTCGCGCGGTGCGAGGATCAGCCACACTgg > 3:47252/1‑149 (MQ=255) CAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTCGCGCGGTGCGAGGATCAGCCACACTGGCAGCAGt < 5:127710‑M1/149‑1 (MQ=255) GTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTCGCGCGGTGCGAGGATCAGCCACACTGGCAGCAGTGc > 6:127710‑M1/1‑149 (MQ=255) cGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTCGCGCGGTGCGAGGATCAGCCACACTGGCAGCAGt < 2:90565/143‑1 (MQ=255) cGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTCGCGCGGTGCGAGGATCAGCCACACTGGCAGCAGt > 1:90565/1‑143 (MQ=255) | AACGCGCGTCGGTTTCATTAGCCAGCAGTTTTGGCGTCGTACCGGAAGAGATCAGCGCGTGGAAGCCAGATACCGCACCACAGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTCGCGCGGTGCGAGGATCAGCCACACTGGCAGCAGTGC > NC_000913/4590190‑4590472 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |