| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,506,020 | 0 | C | A | 41.9% | 57.6 / 16.9 | 43 | G21C (GGC→TGC) | ydcO | BenE family inner membrane putative transporter |
| Reads supporting (aligned to +/- strand): ref base C (18/7); new base A (0/18); total (18/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.66e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01 | |||||||||||
TGAGAACAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGG > NC_000913/1505870‑1506162 | tGAGAACAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGcc > 1:511715/1‑151 (MQ=255) gATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTAATACCGCC‑Ag < 2:416978/151‑1 (MQ=255) tGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGcc > 2:71996/1‑151 (MQ=255) gTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCa < 2:670681/151‑1 (MQ=255) gTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCa < 1:126370/151‑1 (MQ=255) agTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCAGTAGc > 2:154609/1‑151 (MQ=255) agTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGACGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGTCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTATTACCGCC‑AGAAA‑CACCGCCAGTAGc > 2:637760/1‑151 (MQ=255) tCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCAGTAGCgt > 1:471128/1‑151 (MQ=255) tACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCAGTAGCgtgggtgg > 2:239976/1‑151 (MQ=255) cgcCGTCACCCAGCCAGAGTTTTGTGCGGTGGTGGCCCGGGCGACAATC‑GCTGCTTCCAAGATTATTGCCGCTGACC‑GGGCGTAGCAAATAAATCCCGCA‑AAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTc < 1:623822/151‑1 (MQ=255) gcCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTc > 2:674596/1‑150 (MQ=255) ccGTCACCCAGCCAGAGTTTTGTGCAGTGGTGGCTCCGGCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑GGGCGTAGCAAATTAATCCCGCC‑AGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCa < 1:268770/150‑1 (MQ=255) ccGTCACCCAGCAAGAGATTTGTGAGGGGGGGGCCCCGGCGAAAATC‑GCTGCTTGCCAGTTTTTTCCCGCTGAAC‑TGGCGTACCAAATTAATCCCGCC‑AAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCa < 1:55476/150‑1 (MQ=255) ccGTCAACCGGCCAGAGGTTTGTGCAGTGGTGGCTCCGGCGTCAATC‑GCTGCTTGCCAGCTTATTGCCGCTGAAC‑TGTCGTAGCCAATTAATACCGCC‑AGGAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCAc > 2:685715/1‑151 (MQ=255) catccgccaaggtaTTTGTGCAGGTGTGCCTCCGACGAAATCC‑GCTTCTTGCTAGTTTCTGGCCGGTGAGG‑GGGCGTAAAAAAATAACACCCCA‑AGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTc < 2:268596/137‑1 (MQ=255) ccagccagAGTTTTGTGCGGTGGGGGCCCCGGCGACAACC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑GGGCGTAGCAAATAAATCCCGCC‑AAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCg < 1:54932/148‑1 (MQ=255) ccagcaaaAGATTGGGGCGGGGGGGGCTCCGGCACAAACC‑GCTGCTGGCCAGATTATTGCCGGGAAAG‑GGGGTTAGCAAATAAACACCGCC‑AAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTc < 2:707/147‑1 (MQ=255) ccaccaagatttttGTGCGGTGGGGGCCCCGGCAAAACTC‑GCTTCTGCCAAGTTTTTTGCCGCTGACG‑TGGCGTAGCAAATAAACCCCCCC‑AAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGt < 2:686648/138‑1 (MQ=255) gagaTTGGTGGGGTGGTGCCCCCGGCGACAACC‑GCTCCTTCCCAGATTTTCGCCGCTAACC‑GGGGGTAGCAAATAAACACCGCC‑AAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtata < 2:211552/150‑1 (MQ=255) ttttggcgggggtggC‑CCGGCAACAACCGGCTGCTTGCCAGATTGTCCCCG‑TAAACGGGGCGGGGCCAA‑TAATCCCCCCAAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtata < 1:241128/136‑1 (MQ=255) tttGTGCGGTGGGGGCCCGGGCGAAAATC‑GCTGCTTGCCAGATTATTGCCGCGAACC‑GGGCGAACCAAATAAATCCCGCC‑AAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat < 1:255851/145‑1 (MQ=255) gggCAGTGGTGCCCCGGGAAACACTC‑GCGGCTTCCCAGATTATTCCCGCTACCC‑TGGCGTA‑CCAAATAAACCCCCC‑AAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTAc < 1:570000/147‑1 (MQ=255) ggggtggtgCCCCGGGCAAAAATC‑GCTGCTGGCCAGTTTTTTGCCGTTGACG‑GGGCGTAGCAAATAAACACCGCA‑AAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGc < 2:714810/148‑1 (MQ=255) gCAGTGGTGGCTCCGGCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGc > 2:582557/1‑151 (MQ=255) ggggtggCCCCGGCAAAAACC‑GCTGTTGGCAATATTTTTGCCGCTGAGC‑GGGGGTAGCAAATAAAACCCGCC‑AAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCa < 2:634388/147‑1 (MQ=255) gggggCCCCGGCAACACTC‑GCTGCTTCCCAGTTTATTCCCGCTAAAC‑GGGCGTACACAATAAATCCCGCC‑AGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGt < 1:143339/148‑1 (MQ=255) tggCTCCGGCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGt > 1:318502/1‑149 (MQ=255) cGGCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAttt > 2:136311/1‑150 (MQ=255) cGGCAACAATC‑GCTGCTGGCCAGATTTTGGCCGTTGAAC‑GGGGGTAGCAAATTAACCCCGCC‑AAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAttt < 2:440920/150‑1 (MQ=255) gCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTc > 1:557434/1‑151 (MQ=255) gACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCg > 1:721102/1‑151 (MQ=255) aCAATC‑GCTGCTTGCCAGATTTTTCCCGCGAAAC‑GGGGGTAGCAAATTAATCCCCCC‑AAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTcc < 1:295799/149‑1 (MQ=255) c‑GCTGCTTGCCAGATTATTCCCGCTGAAC‑GGGCGAAGCAATTAAACCCCGCA‑AAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTtata < 2:214524/150‑1 (MQ=255) ttGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcac > 1:494607/1‑150 (MQ=255) ttGCCAATTTATTCCCGCTAAAC‑GGGCGTAGCAAATAAACCCCCCA‑AAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGca < 2:738367/149‑1 (MQ=255) caagtttattGCCGCTGAAC‑GGGCGAAGCAAATAAATCCCGCA‑AGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcaca < 2:248843/143‑1 (MQ=255) ttGCCCCTGAAC‑GGGCATAGCAAATAAATCCCGCC‑AAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATaaa < 1:635136/149‑1 (MQ=255) aaC‑TGTCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCCGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCCTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACATGTGATAAACTTCCg > 2:681444/1‑146 (MQ=255) aaC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTg > 1:522599/1‑151 (MQ=255) aaC‑TGGCGTACCCAATTAAACCCCCC‑AGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtt < 1:43873/148‑1 (MQ=255) aaC‑GGGCGTACCAAATTAACCCCGCA‑AAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtttt < 1:304990/150‑1 (MQ=255) aC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGc > 1:177559/1‑151 (MQ=255) tGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGc > 2:705695/1‑149 (MQ=255) ggggtagcAATTAAACCCCCC‑AAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGc < 2:678324/141‑1 (MQ=255) cGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCgg > 2:142471/1‑149 (MQ=255) | TGAGAACAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATC‑GCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATTAATACCGCC‑AGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGG > NC_000913/1505870‑1506162 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |