| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,740,746 | 0 | C | T | 43.9% | 40.7 / 14.5 | 42 | V135V (GTG→GTA) | aroF | 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tyrosine‑repressible |
| Reads supporting (aligned to +/- strand): ref base C (4/19); new base T (18/0); total (23/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.49e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01 | |||||||||||
ATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCAGT‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATT‑AACCCTTTCCAGCCGACAGTGGTACGGGGTTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATC > NC_000913/2740606‑2740895 | aTGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTCACCAGCTCaa < 2:830816/151‑1 (MQ=255)aTGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCACATTTACCCGcc > 1:640191/1‑147 (MQ=255) aGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTCACCAGCTCAAgcagca < 1:15069/151‑1 (MQ=255) aGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTCACCAGCTCAAgcagca < 1:78130/151‑1 (MQ=255) gAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTCACCAGCTCAAGCAGCAATTTAc < 2:4141/151‑1 (MQ=255) cACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCCTAATTACCCACCCCAGCAGCAAATTACGCGCGATc > 1:13189/1‑149 (MQ=255) aGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTTACCAGCTCAAACAGCAATTTACGCGCGAGCTGCAG‑ccc > 2:290328/1‑150 (MQ=255) cGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCCTATTTACCAGCCCAAAAAGCAAATTTCGCGCGGTCTGCCG‑CCCGGc > 1:846125/1‑147 (MQ=255) ccGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTTACCAGCCCAAGCAACAATTTTCGCGCGGTCTGCAC‑CCCGGCtt > 2:237597/1‑144 (MQ=255) gTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCCCTTTTAACCAACCCAA‑CCGCAATTTTCTCGCGGTCTGCGG‑CCCCGCctct > 2:96814/1‑144 (MQ=255) gTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCCTATTTACCAACCCAAACAACAATTTACGCGCGATCTTCAG‑CCCCGCTTCTACATc > 1:886360/1‑149 (MQ=255) tGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCaaa < 1:479008/150‑1 (MQ=255) gTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCaaaa < 1:654664/150‑1 (MQ=255) gTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAag < 2:808504/151‑1 (MQ=255) cGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGc < 1:394535/151‑1 (MQ=255) ccAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGccatcc < 2:666519/150‑1 (MQ=255) tGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGCCATCCATATg < 1:607301/150‑1 (MQ=255) cTGACTAGCTCAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCTGT‑CCCATATTCGCCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGCCATCCATATgg < 1:235775/149‑1 (MQ=255) cTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTTCCCCGCTCCAGCAGCAACTTTCCCCCGATCTGCCG‑CCCCGCCTCTTCCTCACAAAAGCCCTCCATCTtggg > 2:179453/1‑151 (MQ=255) gACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGa < 1:625709/150‑1 (MQ=255) aCCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGAt < 2:109400/150‑1 (MQ=255) ccAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCCTATTCACCAGCACAAGCAGCAATTTACGCGCGGTCTGCAC‑CCCGGCCTCTACATCAAAAAAACCCTCCATATgggggt > 1:664751/1‑147 (MQ=255) cGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATT‑AACCCttt < 2:698860/151‑1 (MQ=255) cGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATT‑AACCCttt < 2:595382/151‑1 (MQ=255) gCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑TCCATATTCAACAGCTCAAACAGCAATTTTCCCGCGGTCTCCAC‑CCCGGCGACTAACTCATAACAACCACACATAATGGGGTC‑GGTTATT‑AACCCttt > 2:180807/1‑150 (MQ=255) cccAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCCTATTTACCAGCTCAAACAGCAATTTAAGCGCGATCTGCAG‑CCCCGCCTCTACATCAAAAAAGGCATCCCTATGGGGGTC‑GTTAATT‑AAACCctt > 2:652395/1‑149 (MQ=255) cccAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATT‑AACCCttt < 1:652444/149‑1 (MQ=255) cAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTTCCCAGCTCAAGCAGCAATTTACGCGCGATCTTCAG‑CCCCGCTTCTACATCAAAAAAGCCATCCATATGGGGGTC‑GTTAATT‑AACCCTTTcc > 1:592212/1‑149 (MQ=255) ggTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCCTATTTACCCGCTCAAGCAGCAATTTACCCGCGATCTGCAG‑CCCCGCCTCTTCATCAAAAAAGCCATCCATATTGGGGGC‑GGTAATT‑AACCCTTTTCAGcc > 1:389241/1‑151 (MQ=255) ttGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCCTATTTACCAGCTCAAGCAGCAATTTTCGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAAAGACATCCATATTGGGGTC‑GTTAATTAAACCCTTTccgcccac > 1:577035/1‑144 (MQ=255) gCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCCTATTTACCCGCCCAAGCAGGAATTTAAACGCGGTCCTCCG‑CCCGGCCTTTACATCAAAAAAGCCATCCAAATGGGGATC‑GTTAATT‑AACCCTTTTCAGCCGACAGTgg > 1:593443/1‑149 (MQ=255) cTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑CCCATATTTACCAACCCCAGCAGCAATTTAAGCGCGGTCTGCAA‑CCCGGCTTCTACATCAAAAAAGCCATCCATATGGGGGTC‑GGTAATT‑TAACCTTTTCAGCCGAAAGTgg > 1:56242/1‑148 (MQ=255) tATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGT‑TCCATATTTATCAGCTCAAACAACAATTTTTGCGCGCGCTGCGG‑CCCGGCCTCTTCATCAAAAAAGCCATCCCTATGGGGGTCGGTTAATT‑AAACCTTTCC‑GCCCACAGTGGGACg > 2:324268/1‑151 (MQ=255) aaCGCTTCCGTCGCCAGTGGCGGT‑CCCCTATTTCCCCGCTCAAGCAGCAATTTTAGCGCGGTCTTCAG‑ACCGGCTTCTACATCAAAAAAGCCCTCCCTATTGGGGTC‑GTTAAAT‑TAACCTTTCCAGCCCACCGTGGGACGGGGgtttt > 1:678662/1‑148 (MQ=255) aCGCTTCCGTCGCCAGTGGCGGT‑CCCCTATTTACCAGCTCAAACAGCAATTTTCGCGCGCGCTGCCG‑CCCCGCCTCTTAATCCAAAAAGCCCTCCCCATGGGGGGC‑GGTTATT‑AAACCCTTCCCGCCGACAGTGGGAAGGGGgttttt > 2:672279/1‑147 (MQ=255) aCGCTTCCGTCGCCAGTGGCGGT‑CCCCTATTTACCAACCCCAACAGCAATTTATGCGCGGTCTGCAGCCCCGGCTT‑TTCATCAAAAAAACCATCCATATGGGGGTC‑GTTTATT‑TAAACTTTTCACCCGACAGGGGGACggggggttttaa > 2:166906/1‑142 (MQ=255) gCTTCCGTCGCCAGTGGCGGT‑CCCATATTCACCAGCTCAAGCAGCAAATTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAAAGCCATCCATATGGGGGTC‑GTTTATT‑AAACCTTTCCAGCCGACAGGGGGACGGGGGTTTTTaaaa > 1:265887/1‑150 (MQ=255) cTTCCGTCGCCAGTGGCGGT‑CCCCTTTTTCCCAGCCGAAACCGGCATTTATGCGCGGTCGGCCG‑CCCGGCTTTTTCATCAAAAAAGACCTCCCTTTTGGGGTT‑GTTAATT‑AACCCTTTCCAGCCGCCAGAGGGACGGGGgttttttaaa > 2:345781/1‑144 (MQ=255) cGCCAGTGGCGGT‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATT‑AACCCTTTCCAGCCGACAGTGGTACGGGGTTTTTCAAAATAGAcgcg < 1:288243/150‑1 (MQ=255) aGTCGCGGT‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGCCATCCACATGGGGATC‑GTTAATT‑AACCCTTTCCAGCCGACAGTGGTACGGAGTTTTTCAAAATAGAcg < 2:832764/144‑1 (MQ=255) cATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATT‑AACCCTTTCCAGCCGACAGTGGTACGGGGTTTTTCAAAATAGACGCGCATTACCAGATAGAgg < 2:71165/151‑1 (MQ=255) ttCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATT‑AACCCTTTCCAGCCGACAGTGGTACGGGGTTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTAt > 2:104285/1‑151 (MQ=255) cACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGTTC‑GTTAATT‑AACCCTTTCCAGCCGACAGTGGTACGGGGTTTTTCAAAATAGACTCGCATTACCAGATAGAGGCTATc < 2:592224/150‑1 (MQ=255) cACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATT‑AACCCTTTCCAGCCGACAGTGGTACGGGGTTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATc < 2:453651/150‑1 (MQ=255) | ATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCAGT‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAG‑CCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATT‑AACCCTTTCCAGCCGACAGTGGTACGGGGTTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATC > NC_000913/2740606‑2740895 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |