Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,096,090 | 0 | A | G | 51.0% | 54.8 / 80.0 | 51 | I83V (ATT→GTT) | hisF | imidazole glycerol phosphate synthase, catalytic subunit with HisH |
Reads supporting (aligned to +/- strand): ref base A (3/22); new base G (26/0); total (29/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.47e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-08 |
GCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGT > NC_000913/2095945‑2096240 | gCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGG‑TTaa > 1:857927/1‑150 (MQ=255) gCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGG‑TTTAGtct > 2:535761/1‑150 (MQ=255) cAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTctct < 2:375097/151‑1 (MQ=255) acgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGG‑TTAAGTCTCTggga > 2:263772/1‑144 (MQ=255) acgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGAt < 2:377988/149‑1 (MQ=255) tGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGG‑TTAAGTCTCTGGGAGATggc > 1:229918/1‑146 (MQ=255) cgcTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGG‑TTAAGGCTCTGGAAGATGCCGCGAAAATTCtt > 1:397602/1‑151 (MQ=255) cTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAGGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGG‑TTTAGTCTCTGGGAGATGCCGCGAAAATTTTTTTCTTTTgcg > 1:463583/1‑151 (MQ=255) tgtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGG‑TTAAGGCTCTGGGAGATGCCGCGAAAATTTTTTTCTTTTGCgg‑ggg > 2:723302/1‑150 (MQ=255) gtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGa < 1:771198/150‑1 (MQ=255) tCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATa < 1:86448/150‑1 (MQ=255) cGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGG‑TTAAGGCTCTGGAAGATGCCGCGAAAATTCTTTTCTTTTGCGC‑GGATAAAAtt > 1:787050/1‑151 (MQ=255) cGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAAtt < 1:832372/151‑1 (MQ=255) gATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGG‑TTAAATCTCTGGAAGAAGCCCCGAAAAATTTTTTCTTTTGCGC‑GGGTAAAAttt > 2:523399/1‑151 (MQ=255) gATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGG‑GTAAGTCTCTGGAAGATGCCGCGAAAAATTTTTTCTTTGGCGC‑GGGGAAAAatt > 1:575482/1‑151 (MQ=255) atatCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTc < 1:625673/151‑1 (MQ=255) cACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCAt < 2:229923/150‑1 (MQ=255) gCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGG‑TTAAGTCTCTGGAAGATGCCCCGAAAATTCTTTTCTTTGGGGC‑GGATAAAAATTCCAt > 1:803526/1‑146 (MQ=255) gCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAAct < 1:144035/151‑1 (MQ=255) aGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGG‑TTAAGGCTCTGGAAGATGCCGCGAAAAATTTTTTCTTTGGCGC‑GGATAAAATTTCCATCAACTCTcc > 2:302398/1‑149 (MQ=255) aGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGG‑TTTAGTCTCTGGAAGATGCCGCGAAAATTTTTTCCTTTGGCGC‑GGGTAAAATTTTCATCAACTCTcc > 1:790180/1‑149 (MQ=255) gTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCaa < 1:341353/151‑1 (MQ=255) tAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTTAAGT‑TCTGGGAGATGCCGCGAAAAATTTTTTCTTTTGGGC‑GGGTAAAATTTTCATTAACACTTCT‑CCGCTGGGAGAacc > 2:559023/1‑148 (MQ=255) aGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGG‑TTAAGTCTCTGGAAGATGCCCCGGAAATTCTTTCCTTTTGCGC‑GGATAAAATTTCCCTCAACTCTCCT‑TCGCTGGCAGACCCaa > 2:830299/1‑149 (MQ=255) aTAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGG‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTTGCGC‑GGATAAAATTTTCATCAACTCTCCT‑GCGCTGGCAGACCCAACatt > 2:837537/1‑151 (MQ=255) aTAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACatt < 2:323012/151‑1 (MQ=255) aaaaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACattaa < 1:901406/151‑1 (MQ=255) aaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGG‑TTAAGTCTCTGGAAGATGCCGCGAAAATTTTTTTCTTTTGGGC‑GGGTAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACattaa > 1:883321/1‑149 (MQ=255) gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACt < 2:750880/149‑1 (MQ=255) gTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAATCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGGTAAAATTTTCATCAACTCTCCT‑GCGCTGGCAGACCCAAAATTAATTACTCGcc > 2:895930/1‑151 (MQ=255) aTCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGG‑TTAAGTCTCTGGAAGAAGCCGCGAAAATTTTTTCCTTTGGCGC‑GGGTAAAAATTTCATCAACTCTCCT‑CCGCTGGCCGACCCCAAAATAATTACTCGcc > 1:866221/1‑149 (MQ=255) aTCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGG‑GTAAGTCTCTTGGAGATGGCCCGAAAATTCTTTCCTTTGGCGG‑GGATAAAAATTTCATCAACTCTCCTGGCGC‑GGCCGACCCAACACTAATTAATCgtccg > 1:496039/1‑147 (MQ=255) tctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGG‑TTAAGTCTCTGGGAGATGCCGCGGAAATTTTTTTCTTTGGCGC‑GGGTAAAATTTCCATCAACTCTCCT‑GCGCTGGCCGACCCCACATTAATTACTCGcc > 1:912555/1‑148 (MQ=255) cgcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCTg < 1:662858/147‑1 (MQ=255) cgcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCTg < 2:51871/147‑1 (MQ=255) cgcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGG‑TTAAGTCTCTGGGAGATGCCGCGAAAATTCTTTTCTTTGGCGC‑GGGTAAAATTTTCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCTGGcc > 2:325631/1‑150 (MQ=255) gTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATc < 1:538491/151‑1 (MQ=255) ggAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGG‑TTAAGTCTCTGGGAGATGCCGCGAAAATTCTTTTCTTTGGCGCGGGGAAAAATTTCCTTC‑ACTCTCCT‑GCGCTGGGAGACCCCACATTAATTTTTCCCCTTGCCGATCGcct > 1:168960/1‑148 (MQ=255) gAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTg < 2:579437/151‑1 (MQ=255) aGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCCGGCCGATCGCttt > 2:578809/1‑148 (MQ=255) tGATCGACATTCCGTTTTGTGTGGCGGGTGGGG‑TTAAGTCTCTGGGAGATGCCGCGAAAATTTTTTCCTTTTGCGC‑GGATAAAATTTTCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCg > 2:579148/1‑150 (MQ=255) tGATCGACATTCCGTTTTGTGTGGCGGGGGGGG‑TTAAGTCTCTGGGAGATGCCGCGAAAATTCTTTTCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCg > 2:399582/1‑150 (MQ=255) tCGACATTCCGTTTTGTGTGGCGGGGGGGG‑TTAAGTCTCTGGGAGATGCCGCGAAAATTTTTTTCTTTTGCGC‑GGGTAAAATTTTCATCAAATCTCCT‑GCGCTGGCAGACCCAACATTAATTAATCGCCTGGCCGATCGCTTTGGGGTg > 1:571547/1‑149 (MQ=255) tCGACATTCCGTTTTGTGTGGCGGGGGGGG‑GTTAGTCTCTGGAAGATGCCGCGGAAAATTTTTTCGTTGGGGC‑GGGTAAAATTTCCATTAACTCTCCT‑GGGCCGGCAGAACCAACATTAATTACTCGCCTGGCCGACCGCTTTTGCg > 1:135167/1‑147 (MQ=255) aCATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAgtg < 2:857937/151‑1 (MQ=255) aCATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAgt < 2:47201/150‑1 (MQ=255) aTTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTa < 1:762669/151‑1 (MQ=255) gTTTTGTGTGGCGGGGGGGG‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGGTAAAATTTTCATCAACTCTCCT‑GCGCTGGCAGACCCCACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATtgt > 1:604397/1‑150 (MQ=255) gCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTg < 1:830472/151‑1 (MQ=255) gggA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCt < 1:513700/151‑1 (MQ=255) gggA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCt < 2:814638/151‑1 (MQ=255) a‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGt > 2:546404/1‑151 (MQ=255) | GCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGA‑TTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCT‑GCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGT > NC_000913/2095945‑2096240 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |