Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,984,718 | 0 | G | C | 36.5% | 72.5 / 24.6 | 56 | P292R (CCG→CGG) | yqeF | short chain acyltransferase |
Reads supporting (aligned to +/- strand): ref base G (7/26); new base C (19/0); total (26/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.72e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01 |
CCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCT > NC_000913/2984571‑2984867 | ccaTTGACATTGACCCGACGCTCATCCCACTCAAGCATCGTGCCAACCGAAAGCGCCTGGGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGgc < 2:180947/150‑1 (MQ=255) cattgacattgaCCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGgcg‑c < 1:1157220/151‑1 (MQ=255) cattgacattgaCCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGgcg‑c < 1:897247/151‑1 (MQ=255) attgacattgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCgg > 1:1107353/1‑147 (MQ=255) gacattgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAAtt < 2:384739/151‑1 (MQ=255) acattgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACCCCGCCG‑Caa > 1:291386/1‑148 (MQ=255) tgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCa < 1:645994/150‑1 (MQ=255) gacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCAt < 1:1212929/150‑1 (MQ=255) cGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGGCTAAACCCGCG‑CCATTTCCCTCAATGcc > 1:248048/1‑150 (MQ=255) tCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCAACCACCCCCG‑CACCTCCCACCAATACCgg > 2:493226/1‑146 (MQ=255) aaGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACAACCGCG‑CAAATTCCCCCAAAGCCGGATCTACACCa > 1:1187912/1‑148 (MQ=255) aaGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCGTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCGACGGCTGGTCGCATTCTCCGGCG‑CAATTCGCAGCAATGCCGGATCTACACCaa < 2:1019415/149‑1 (MQ=255) aGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGc < 1:643847/151‑1 (MQ=255) tCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCgg > 1:872638/1‑112 (MQ=255) cTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑Caa < 2:872086/116‑1 (MQ=255) gCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCaaa < 1:829899/150‑1 (MQ=255) gCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATAAACCGGCG‑CAATTCCCACCAATTCCGGATCTACACCAACGCTGGCcaa > 1:628892/1‑150 (MQ=255) cAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCG‑CCAATCCCCTCAACACCCGGGCTACACCCACGCCgggcaatac > 2:161779/1‑144 (MQ=255) cAACCGAAAGCCCCTGTGCAGTAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACTCTCCAGGCAACGGCGGGTCGCATACACCGTCA‑CAATTCCCATCAATGCCAAATCTACACCAACGCTGGCa < 2:109757/146‑1 (MQ=255) aaCCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATgcg < 1:403390/151‑1 (MQ=255) acggAAAGCGGCTGTGCAGCAAACGATTCATTAACATCGATAACATCGACTTGCGCCACCCTGCC‑GCCTACACGCGCCACGCAACGGCTGGTCGCATACACCTTGA‑CAATTCCGCTCAATCCGGGATCTACACCAACGCTGGCAAATgcgc < 2:1218723/148‑1 (MQ=255) gAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATAAACCCGCG‑CAATTCCCCTCAATGCCGGATCTACACCAAAGCTGGCAAAAgcg > 1:883028/1‑147 (MQ=255) aaaGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CCATTCCCATCAATGCCGGGTCTACACCAACGCCGGCAAAAGCGCGaa > 1:841959/1‑150 (MQ=255) aaGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCAAACACCCGCG‑CAAATCCCCTCAAAGCCGGGTCTACACCAAACCCGGCAAAAgc > 1:504365/1‑144 (MQ=255) gcCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCC‑CAATTCCCCTCAATGCCGGGTCTTAACCAAACCTGGCCAAATCGCGAATGCggg > 1:1077751/1‑151 (MQ=255) ccTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATAAACCCGCG‑CAATTCCCATCAATGCCGGGTCTACACCAACCCCGGCAAAAAAGCGAATGCggg > 1:48913/1‑150 (MQ=255) gtgCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCa < 1:920271/150‑1 (MQ=255) gcaAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CCATTCCCATCAAAGCCGGGTCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑c > 1:137862/1‑147 (MQ=255) aaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCCCATACAACCGCC‑CCATTCCCCTCACAGCCGGATCTACACCAACGCTCGCACAAGCGCGGAAGCGGGGCAG‑ccacgg > 2:484827/1‑144 (MQ=255) gCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCACACACCCGCG‑CCATTCCCCTCAATTCCGGAGCTACACCAACGCTGGCCAAAGCCCGAACGCGGGCCGC‑CACGGGGaaa > 1:753671/1‑150 (MQ=255) ttCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CAATTCCCATCCATGCCGGGTCTACACCAACCCTGGCAAATGCGCGAAAGCGGGCCAG‑CACCGGGaaa > 1:387810/1‑148 (MQ=255) cATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CAATTCCCATCAATGCCGGATCTTCACCAACGCCGGCAAAAGCGCGAAAGCGGGCCAG‑CACGGGGAAATTCa > 1:1118611/1‑150 (MQ=255) cATTAGCCTCGATAAGAGCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTCACCGGCA‑TAATACGCAGCACTGCGAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCaa < 2:1111702/151‑1 (MQ=255) tAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACg < 2:1041450/150‑1 (MQ=255) ccTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCg < 1:754842/151‑1 (MQ=255) cGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCGCCAGGCAACGGCGGGTCGCATACTTCTGAG‑CGATACCCCGCACAGCCAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCGTg < 2:360743/150‑1 (MQ=255) cGATAAGATCGACTTCAGCAA‑ACTGCCAGCCTGCACGCTCCAGGCAACGGCGGGTCGCATTCTCCGGCG‑TAATACCCAGAAAAGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCt < 2:1260448/146‑1 (MQ=255) gATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CAATTCCCCTCAATGCCGGATCTACACCAACGCCGGCAAATGCGCGCATGAGGGCCCG‑CACCGGGAAATTCAACGCTCGTGc > 1:91499/1‑150 (MQ=255) aTAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGGATAAACCCGCC‑CCATTCCCCTCAAAGCCCGGACTACACCCACCCTGGGAAAAGCGCGAATGCGGGGCAG‑CACAGGGGAAATCAAAGCTCGTGCtt > 1:106471/1‑151 (MQ=255) aGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCGTGCTTTgg < 1:45018/151‑1 (MQ=255) aTCGACTGCAGCCA‑ACTGCCAGCCTTCACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCGTGCTTTGGCt < 2:750925/151‑1 (MQ=255) tCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CAAATTCCATCAATGCCGGGTCTACACCCACGCCGGCAAAAGCGCGAATGCGGGCCAG‑CACCGGGGAATTTAAAGCCCGTGCTTTTGCtt > 1:868336/1‑151 (MQ=255) cAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CCATTCCCCTCCATGCCGGATCTACACCAACGCTGGCAAAAGCGCGGAAGCGGGCCAG‑CACGGGGAAAATTAAACCTCGTGCTTTTGCtt > 1:28935/1‑144 (MQ=255) gCCA‑CCTGCCAGACGACACGCCCCAGGCAACGGCGGGTCGCATTCACCGGCG‑CGATTCCCATCCATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtc < 2:293048/147‑1 (MQ=255) cA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGCGGGCTACACCCGCC‑CCACTCCCCTCCATGACCGGTCTCCACCCAAGCTGGCCAATGCGCGGATGCGGGCCAG‑CCCGGGGGAATTTCACGCCCGTGGttttgggtctctc > 2:1259405/1‑135 (MQ=255) cA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGCCGCCTACACCGGCG‑CCATTCCCCTCAATGCCGGACCCCCACCAACCCTGGGCAATCCGCGGATTCGGGCCCG‑GCCGGGGCAATTTAACGCACCGGCTTTGGCTTCACtcctcat > 1:1153613/1‑150 (MQ=255) a‑aCTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCTGCG‑CGATTCCCATACGTGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcatcatca < 2:1280234/151‑1 (MQ=255) aGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGc < 1:901718/151‑1 (MQ=255) gCCTACACGCTCCAGGCAACGGCGGGGTGCATCAACCGGCG‑CAATTCCCACCCATGGCGGATCTACACCAACGACGGCCAAATCCCGGATTCGGGGCCG‑CCCGGGCAAAATCAACCCTCGTGCCTTGGCGTGGCTCATCAACATTAcccc > 2:1026671/1‑148 (MQ=255) acGCTCCAGGCAACGGCGGGGCGCATACACCGGCGCCAATTCCC‑TCAACGCCCGATCTACACCAACGCTGGGAAAAGCGGGAATTGGGGCCAGCCACGGG‑AAAATCAAAACCCGTGCTTTTGCTTCTCTCCTTATCATTACCCCAGCTg > 2:221909/1‑149 (MQ=255) tCCAGGCAACCGCGGGTCGCATTCTCCTTAG‑CAATACCCCGCCCACCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTgcg < 2:313491/147‑1 (MQ=255) gCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCa > 2:1236766/1‑150 (MQ=255) gCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCAt < 1:946367/151‑1 (MQ=255) acacCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTAcc > 1:487064/1‑151 (MQ=255) cacCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTAccc < 1:423978/151‑1 (MQ=255) cacCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTAccc < 1:1261552/151‑1 (MQ=255) ccGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCg < 1:1030329/150‑1 (MQ=255) cGGCG‑CAATTCCCATCAATCCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCTTTCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCt < 2:706433/151‑1 (MQ=255) | CCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCA‑ACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAG‑CACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCT > NC_000913/2984571‑2984867 |
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: ‑ |