| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | ACB143 | 4,772,546 | 0 | T | G | 29.6% | 22.4 / 15.3 | 56 | intergenic (‑140/+21) | grsB/lgrB_1 | Gramicidin S synthase 2/Linear gramicidin synthase subunit B |
| Reads supporting (aligned to +/- strand): ref base T (4/34); new base G (14/2); total (18/38) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.28e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTG > ACB143/4772422‑4772662 | ccgggCACCGGGCTTTGCCTGCGGAACGTGCCAGCCCGCCACAGGCTCATGAGGCTGCCGAACAAATGACAGAGCCGGACGGCGGCCGGCTCGGGGTGTACTGAAGTGGACGAG‑CC‑GGG‑CCCGGTCATACC‑GGGc < 7:286867/133‑1 (MQ=255) aGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑GCCGGGGCTAAC‑GGGcc > 2:302650/1‑135 (MQ=255) gCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGGCCCCGGGAAAACC‑GGGccgc > 8:102528/1‑133 (MQ=255) aCCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGCG‑CCGGGG‑CCCGCGCACACC‑CGGcccccc > 3:362209/1‑132 (MQ=255) aCCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGTCATACC‑GTGc < 1:118207/129‑1 (MQ=255) aCCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGTCATACC‑GTGc > 2:118207/1‑129 (MQ=255) aCCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGGGAGTCC‑GGGCCCGcgg > 2:45303/1‑133 (MQ=255) aCCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCCGTCCTTCC‑GGGCCCGccg > 5:95938/1‑133 (MQ=255) aCCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCCGGCCTACC‑GGGCCCGccc > 8:224229/1‑133 (MQ=255) aCCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAAC‑AATTCAAGCGCCCGCCCGCAGCCGGCGAGGGGGGGACAGAAGGTGCCGCGCCC‑GGG‑CCCGGGGAAACC‑GGGccccccc > 4:120346/1‑131 (MQ=255) gacTTTGCCTGCGGACCGGGCCCGCCCGCGACTGGCTAATGATCCTGCCTAACAAATTAACGCGCCGGACGGCGGCCGGCTCTGGGTGTACTTAAGTTGCCGAG‑CC‑GGC‑CCCGGTCATACC‑GGGCCCGCAGAgcg < 2:235210/133‑1 (MQ=255) tttGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCCGGCAGGCC‑GGGCCCCCCgcgcgcga > 7:91047/1‑133 (MQ=255) gCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑GCCGGGCATAAC‑CGGCCCCCACCgcgcgaagg > 7:68158/1‑130 (MQ=255) cTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGGGCTGAC‑GGGCCCCCAGAGCGCGACTggg > 2:110206/1‑135 (MQ=255) cgcaaCGTGCCCGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGGGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCa < 7:343485/132‑1 (MQ=255) ggAACGTGCCAGCCCGCTACAGGCTCATGATGCTGCCGAACAAATGAAAGCGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCa < 3:297464/134‑1 (MQ=255) aaCGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑GCCGGGCAATCC‑GGGCCCCCCCAGCCCGACTggggggaag > 8:78928/1‑129 (MQ=255) aaCGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCCGGCATACC‑GGGCCCCCCGCGCGCGACAGGGGGAACt > 8:194667/1‑135 (MQ=255) cGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAAAAAATGAACGAGCCGGAGGGCGGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACtcg < 5:283894/135‑1 (MQ=255) tcccAGCCCGCGACAGGCTCATGAGGCTGCCCAAAAAAGGAAACCGCGGGACGGCGGCCGGCTCTGGGTGTACTGATGTTGACGAG‑CC‑GGG‑CCCGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACtcgtc < 2:4282/133‑1 (MQ=255) tGCCCGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAACGCGCCGGAGGGCGGCCGGCTCTGGGTGTACTGAAGTTCACGCG‑CC‑GGG‑CCCGGTCATCCC‑GGGCCCGCAGAGCGCTACTGGGTCAACtcgtc < 1:229148/135‑1 (MQ=255) tGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACtcgtc < 7:227552/135‑1 (MQ=255) tGCCACCCGGCTAAAGGCTCTGAATGCTGCCGAACAATGCAATGCGCGGGAGGGCTGCCGGCTCGGGGGGTACTTAAGTGCCCGCG‑CC‑GGG‑CCCGGTAATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACtcgtc < 4:232993/135‑1 (MQ=255) ggCCAGCCCGCGACAGGCTCATGATGCGGCCGAACAAATGAACGCGCCGGACGGCGGCCGGCTCTGGGTGTACTGAAGTGGACGAG‑CC‑GGC‑CCCGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACtcgtc < 1:41439/134‑1 (MQ=255) aGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGGGCTACC‑GGGCCCGCCGAGCGCGACGGGGTGAACtcgggggtc > 2:35589/1‑129 (MQ=255) cccGCGACAGGCTCATGATCGTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTTTACTGAAGTTGACGAG‑CG‑GGG‑CCGGGTAATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGAtctc < 7:261491/135‑1 (MQ=255) cccGCGACAGGCGCATGATGCTGCCGAACAAATGAACGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGAtctc < 1:50133/135‑1 (MQ=255) ggCTAATCATGGCGCCCAACAAATCAACCCCCCGGCCGGCGGCCGGCTCTGGGTGAACGGATGTTGACGAG‑CC‑GGC‑CCCGGTAAAACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGt < 5:96594/135‑1 (MQ=255) cATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGGCCTACC‑CGGCCCGCCGCGCGCGACGGGGGGCACACGGCGCTCACCGCGGCAGTAAGa > 1:308430/1‑135 (MQ=255) gctgcCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTc < 2:207236/135‑1 (MQ=255) gAACAAATGAACGCGCCGGACGGCGGCCGGCTCTGGGTGAACTGAAGTTGCCGAG‑CC‑GGG‑CCCGGTAATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggcg < 1:303582/135‑1 (MQ=255) aacatttGAATGAGCCGGAGGGCTGCGGGCTCTGGGTGTACTGATGTTGCCGAG‑CG‑GGG‑CCGGGTCATACG‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggcgg < 2:290900/129‑1 (MQ=255) aaCAAGTGAACGCGCCGGCGGGCGGCGGGCTCGGGGTGTACTGAGGTGGCCGCG‑CG‑GGG‑CCGGGCCATCCC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggcgg < 7:230334/135‑1 (MQ=255) aaCAAATTAAAGCGCCGTACGGCTGCCGGCTCTGGGTGTACTGAAGTGCACGAG‑CC‑GGG‑CCCGGTAATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggcgg < 1:203882/135‑1 (MQ=255) aaCAAATGAAAGAGCTGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGTGATACC‑GGGCCCGCAGCGCGCGACTGGGGGAACGCGTCGATCTCCGCTTGAGTAAGACGTTTGAAggcgg > 2:250582/1‑135 (MQ=255) aTGAACCAGCCGGACGGGTGCCGGCTGTGGGTGTACTTAAGTGGACGCG‑CC‑GGG‑CCGGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAg < 1:277527/135‑1 (MQ=255) aTGAAAGCGCCGGACGGCTGCGGGGTCTGGTTGAACTTATTTTCCCGAG‑CC‑GGG‑CCCGGTAATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAg < 2:90478/135‑1 (MQ=255) aaaGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGGCAAACC‑GGGCCCCCAGCGCGCGACTGGGGGACCTCGGCGATCTCCGCTGGAGTAACGCGTTTGAAAGCGGCCCATGAt > 7:190283/1‑135 (MQ=255) gacccGGAGGGCTGCCGGCTCGGGGTGTACTGACGTTCACGAG‑CG‑GGG‑CCGGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTa < 1:22171/132‑1 (MQ=255) ccggccggCGGCCGGCTGGGGGGGAAGGCACGTTGCCGCG‑CC‑GGG‑CCCGTTCATCCC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGcc < 8:205984/135‑1 (MQ=255) gggCTGCGGGCTGTGGGTGAAGTCATGTTCCCGCC‑CC‑GGG‑‑CCGGTTATACCGGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTc < 7:55329/134‑1 (MQ=255) ctggtgtGGGTGAACTCATGTTTCCGAG‑CG‑GGG‑CCCGTTCATACG‑GGCCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGAAGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCt < 2:249330/127‑1 (MQ=255) ggtgggggTGTAATGAGGTGCCCGAG‑CG‑GGG‑CCGGGTCATCCG‑GGCCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAAt < 4:99536/130‑1 (MQ=255) cgtttGGGGGTACTGAAGTTCACGCG‑CC‑GGG‑CCCGGTAATACG‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAAt < 6:319703/131‑1 (MQ=255) tctGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGTCAAACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCg < 5:162209/135‑1 (MQ=255) ggTGTACTGACGTGGACGAG‑CC‑GGG‑CCCGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGAcgc < 5:124845/135‑1 (MQ=255) gtgtACTGAAGTTGCCGAG‑CC‑GGG‑CCCGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGAcgc < 3:146200/134‑1 (MQ=255) tttACTGAAGTGCACCAG‑CC‑GGG‑CCCGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACgcgc < 3:303574/133‑1 (MQ=255) tgtACTTAAGTGCCCGAG‑CC‑GGG‑CCCGGTCAACCC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACgcgc < 5:139425/135‑1 (MQ=255) tgtACTGATGTGCACGAG‑CC‑GGG‑CCCGGGCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACgcgc < 2:226343/135‑1 (MQ=255) tACTTACGTTCACGAG‑CC‑GGG‑GCCTGTAATACG‑GGGCCCGCAGAGCGCTACTGGGTGAACTCGTAGATCTCCTCTGCAGTAAGACGTTTCATGGCGGCCTAGGATTTAGCCTACTCGTCCTGAATCGACGCgccc < 8:280951/135‑1 (MQ=255) aCTGACGTTGACGCG‑CC‑GGG‑CCCGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCTGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccc < 8:99217/134‑1 (MQ=255) gAAGTTGACGAG‑CC‑GGG‑ACCGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTACTCTGCAGGAAGAAGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgcccagcc > 5:65036/1‑135 (MQ=255) aGTTGACGAG‑CC‑GGG‑CCCGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccctgccct < 5:143467/135‑1 (MQ=255) gttcccGAC‑CC‑GGG‑CCCGGGCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccctgccct < 4:100889/129‑1 (MQ=255) tgccGAG‑CG‑GGG‑CCCGGCCAACCG‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTg < 2:32473/131‑1 (MQ=255) | CAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAG‑CC‑GGG‑CCCGGTCATACC‑GGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTG > ACB143/4772422‑4772662 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |