| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,984,441 | 0 | A | C | 56.7% | 23.7 / 41.8 | 30 | G384G (GGT→GGG) | yqeF | short chain acyltransferase |
| Reads supporting (aligned to +/- strand): ref base A (13/0); new base C (0/17); total (13/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.35e-09 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-06 | |||||||||||
ATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAA‑TGCAGAGC‑TATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCC > NC_000913/2984296‑2984586 | atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAA‑TGCAGAGC‑TATTCGTCACGTTCAATGGTC‑AATGCCACACCCTg > 1:200732/1‑151 (MQ=255)atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAA‑TGCAGAGC‑TATTCGTCACGTTCAATGGTC‑AATGCCACACCCTg > 2:316453/1‑151 (MQ=255) gAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAA‑TGCAGAGC‑TATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATAc > 1:432831/1‑151 (MQ=255) aaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAA‑TGCAGAGC‑TATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACa > 1:290149/1‑150 (MQ=255) ccaaaaaaaTTTATGTCACCTTGCCAGCAATCATTATCAAAAAAATAAATGCCGGCAGATTTTTATGGGAAGGGAGAAAA‑TGCAGAGC‑TTTTCGCCACGTTAAATGGTC‑AATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCa < 1:316451/144‑1 (MQ=255) aGAAAATTTTGGTCGTCTGCCCAGCCATCATTATTCAACAGATTAATTACGGCAAATTTTTAGCGGTAAGGAGAAAA‑T‑CAAAGCTTTTTCGCCACTTTCAATGGTC‑AAGGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCt < 2:446734/151‑1 (MQ=255) aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGCTTAATGACGGCAGATATTTATCGGTAAGGAGGAAA‑TGCAGAGC‑TATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcat > 1:335851/1‑149 (MQ=255) aaatTGTATGTCGACTGCCCAGCCATCATTATTCAACAAATTAATGCCGGCAAATATTTCTCGG‑AAGGAAGAAATTGCAGTTT‑TTTTCGTCCCGTT‑AATGGGCAAATGCCACCCCCTGGCCACCGCCGATACAAAGCGTTGCCAGTCCttt < 2:199650/151‑1 (MQ=255) aaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAAATTAATGCCGGCAGATATTTATGGGTAAGGAGGAAA‑TGCAGAGC‑TTTTCGTCACGTTCAATGTTC‑AATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCttt < 2:442270/151‑1 (MQ=255) aaTGTATGTCGACTTCCCAGCCATCATTATTAAACAGATAAATGACGGAAGATATTTATCGGTAAGGAGAAAA‑TGCAGAGT‑TTTTCGTCACTTTCAATGGTC‑AATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc < 2:111907/151‑1 (MQ=255) cGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGAAAGGAGAAAA‑TGCAGAGT‑TATTCGTCACGTTCAAGGTCC‑AATCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt < 2:13547/149‑1 (MQ=255) ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAA‑TGCAGAGC‑TATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt > 1:165649/1‑151 (MQ=255) ccaTCTTTATTAAACAGATTAATGACGGCAGATATTTAGCGGAAAGGAGAAAA‑TGCAGAGC‑TATTCGTCACTTTAAATGGTC‑AATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa < 2:393326/149‑1 (MQ=255) ttattaAACAGATTAATGACGGCAGATTTTTATCGGTAAGGAGGAAA‑TGCAGAGC‑TTTTCGTCACGTTCATTGGTC‑AACCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt < 2:325211/146‑1 (MQ=255) ttattCAACAGATTAATGACGGAAGATATTTATGGGTAAGGAGAAAA‑TGCAGAGC‑TTTTCGTCACTTTCAATGGTC‑AATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt < 2:388697/149‑1 (MQ=255) aaCAGATTAATGACGGCAGATTTTTATCGGAAAGGAGAAAA‑TGCAGAGT‑TATTCGTCACGTTAAATGGTC‑AATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc < 2:432836/148‑1 (MQ=255) aaCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAA‑TGCAGAGC‑TATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc > 2:510781/1‑148 (MQ=255) aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAA‑TGCAGAGC‑TATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc > 1:524870/1‑151 (MQ=255) gcgggaaGATATTTATGGGAAGGGGGAAAA‑GGCAGAGC‑TATTCGTCACGTTCAATGGTC‑AATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg < 1:391480/143‑1 (MQ=255) cAGATATTTATCGGTAAGGAGGAAA‑TGCAGAGC‑TATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc > 2:27778/1‑151 (MQ=255) tttttAGCGGTAAGGAGAAAA‑TGCAGAGT‑TTTTCTTCTCTTTAAATGGTA‑AATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCa < 2:533150/146‑1 (MQ=255) tttATCGTAAAGGGGAAAA‑TGCAGAGC‑TTTTCGTCATGTTAAAGGGCC‑AATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc < 1:261030/149‑1 (MQ=255) tttATCGGTAAGGAGGAAA‑TGCAGAGC‑TATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg > 2:598444/1‑151 (MQ=255) tttATCGGTAAGGAGAAAA‑TGCAGAGC‑TTTTCGTCACGTTCAAGGGTC‑AATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc < 2:87408/149‑1 (MQ=255) gTAAGGAGGAAA‑TGCAGAGC‑TATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGCCATTACGTTTCACCCTTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc > 1:603962/1‑148 (MQ=255) gagC‑TTTTCGTCACGTTAAAGGGTC‑AATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt < 1:205757/151‑1 (MQ=255) ttttCGCCACGTTAAGGGGTC‑AATCCCCCCCCGTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg < 1:552526/146‑1 (MQ=255) ttCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc > 1:653157/1‑150 (MQ=255) aTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCattgacattgac > 2:495929/1‑151 (MQ=255) ccACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAccc < 2:119811/150‑1 (MQ=255) | ATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAA‑TGCAGAGC‑TATTCGTCACGTTCAATGGTC‑AATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCC > NC_000913/2984296‑2984586 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |