| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,953,555 | 0 | A | G | 58.8% | 10.7 / 17.4 | 19 | G26G (GGA→GGG) | ilvD | dihydroxyacid dehydratase |
| Reads supporting (aligned to +/- strand): ref base A (1/6); new base G (10/0); total (13/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.66e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-04 | |||||||||||
AAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACAC > NC_000913/3953419‑3953704 | aaaaaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGCCGTGGCGCGCCCCCCGGAAGACCCACG‑CC‑g > 2:45480/1‑150 (MQ=255) ggACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCCGTGGGGCGCCACCCGGAAGAACCACG‑CC‑gcattttgg > 1:54452/1‑143 (MQ=255) ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGGTCGTTCGCTGTGGCGGGCCCCCCGGAAGAACGCCG‑CC‑GATTTTGGTAAgg > 1:250354/1‑150 (MQ=255) ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGCTGTGGGGCGCCCCCGGGATGACCGACG‑CC‑GCTTTCGGTGAg > 2:240588/1‑150 (MQ=255) gTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGGTCGGGCGCTGTGGGGCGCCCCCCGCATGACCGGCC‑CC‑GGTTTTGGGGAGCCggatttcg > 1:264959/1‑143 (MQ=255) acacTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGCCGTTCCCTGTGGCGCGCCCCCGGCATTACCAACG‑CC‑GGCTTTGGTGAGCCCGTTATTGCGGTtgggggctccgtca > 2:217852/1‑138 (MQ=255) taaataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTg < 2:138293/149‑1 (MQ=255) aaataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTGTAc < 1:73488/151‑1 (MQ=255) ataaaGTATGCCTAAGTATCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCCGGGGCGCGCCCCCCGGAAGACCGACC‑CC‑CACTTTGGTTAGCCGAGTATCGCCGGGGTGGACTCGTTCACCCCAAtttttccg > 2:260709/1‑146 (MQ=255) ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGTGTCCTGCGCTGTTGCGCGCCCCCGGGAATACCCACC‑CC‑CGTTTTGGTAAGCCCATTATACCCGCTGTGGACTCGTTCA‑CCCAAttttttacccgtcccgcg > 1:255386/1‑135 (MQ=255) ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCGCGTTCGCCGTGGCGCGCCCCCCGGAAGAACCACG‑CC‑GCTTTTGGTAAGCCGGTTATCGCGGGTGTGGACCCGTTCT‑CCCCAATTTTACCGGGGCACgcc > 2:94109/1‑148 (MQ=255) ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCCCGTTCGCCGTGGGGCGCCCCCCGGAGGAACGACG‑CCGGGTTTTGGT‑AGCCGGGTATTGCGGGTGTGGACTCCTTCC‑CCCCATTTTTACCggggggcgtc > 2:94511/1‑143 (MQ=255) aaGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTGTACCGGGTCACGTCCATc < 1:59510/151‑1 (MQ=255) aaGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTGTACCGGGTCACGTCCATc < 1:307924/151‑1 (MQ=255) aaGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTGTACCGGGTCACGTCCATc < 1:288633/151‑1 (MQ=255) ttCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCCGGGGGGCGCCCCCGGGAAGACCCACGCCC‑GATTTGGGTAGGCGGATT‑TTGCGGTTGGGGAATTGTTCA‑CCCCATTTTTTCCGGGGCACGGCCCCCTTcgcggg > 2:7158/1‑149 (MQ=255) ccaccacTCATGGTCGTAATATGGCGGGGGGTCGTGCGCTGTGGCGCGCCCCCCGGATGACCCACG‑CC‑GCTTTTTGTAAGCCGGTTATCGCGGGTGGGGACTCGTTCC‑CCCCCTTTTTTCCCGGGCACGGCCCTCTTCGCGGtctccggga > 2:179581/1‑146 (MQ=255) gTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCgg > 2:314517/1‑150 (MQ=255) cGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAac < 2:97974/151‑1 (MQ=255) aaTGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAacac < 2:241740/150‑1 (MQ=255) | AAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACAC > NC_000913/3953419‑3953704 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |