Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,922,736 | C→A | 10.5% | R29L (CGC→CTC) | PP_4331 ← | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,922,736 | 0 | C | A | 10.5% | 110.8 / 3.6 | 38 | R29L (CGC→CTC) | PP_4331 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (16/18); new base A (2/2); total (18/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01 |
AGGGTGGCCACCGCCGCACGCAGCTCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTCCTCGGTCATCATCTTGTCCAGTTCGACC > NC_002947/4922592‑4922880 | aGGGTGGCCACCGCCGCACGCAGCTCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTgg < 1:244843/150‑1 (MQ=255) ccACCGCCGCACGCAGCTCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGc > 1:777546/1‑150 (MQ=255) cgccgcACGCAGCTCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTa < 1:517561/150‑1 (MQ=255) gCTCATGGATGGCCTCACCCGTGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGaaaa < 1:815949/150‑1 (MQ=255) tCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTTGCGCTGTCTCTTGCTGTAGTTTAGGAAAAACCACAGGCTCAgc < 1:322751/150‑1 (MQ=255) tCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCATTGAGGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAgc < 2:542702/150‑1 (MQ=255) aTGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAg > 2:303811/1‑150 (MQ=255) gTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGaaa < 2:39053/109‑1 (MQ=255) gTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGaaa > 1:39053/1‑109 (MQ=255) cAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCTCCCAGGCCAgc < 2:634891/127‑1 (MQ=255) cAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAgc > 1:634891/1‑127 (MQ=255) cAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACctc > 2:292962/1‑150 (MQ=255) aCGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAg < 1:303811/150‑1 (MQ=255) aCGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAg < 1:292962/150‑1 (MQ=255) ccAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAg < 2:224458/110‑1 (MQ=255) ccAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAg > 1:224458/1‑110 (MQ=255) ccAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACt > 1:82955/1‑150 (MQ=255) tCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGAt < 2:82955/150‑1 (MQ=255) cTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAg > 1:803231/1‑114 (MQ=255) cTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAg < 2:803231/114‑1 (MQ=255) cTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGt > 2:735010/1‑150 (MQ=255) cTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGt > 2:648270/1‑150 (MQ=255) gcgcTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTGTCCa < 2:638298/150‑1 (MQ=255) tagTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTc > 2:807090/3‑150 (MQ=255) cACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGAc < 1:735010/150‑1 (MQ=255) cccTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCAc < 1:74373/139‑1 (MQ=255) cccTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCAc > 2:74373/1‑139 (MQ=255) ccTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCa > 1:312212/1‑127 (MQ=255) ccTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCa < 2:312212/127‑1 (MQ=255) ccTCGGCCTGCTGCGCAGCCAGCTCGAGCTGGCGCTTGCTGTAGTTGAGGAAAAAc > 1:361178/1‑56 (MQ=255) ccTCGGCCTGCTGCGCAGCCAGCTCGAGCTGGCGCTTGCTGTAGTTGAGGAAAAAc < 2:361178/56‑1 (MQ=255) tCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGa > 1:335847/1‑46 (MQ=255) tCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGa < 2:335847/46‑1 (MQ=255) tgctgcGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTCCTCGGTCATCATCtt > 2:461521/1‑150 (MQ=255) ctgcGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTCCTCGGTCATCATCTTGt < 1:648270/150‑1 (MQ=255) aagcgccagcGCGAGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCGGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTCCTCGGTCATCATCTTGTCCa > 2:491720/5‑150 (MQ=255) agcgccagcGCGAGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCGGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTCCTCGGTCATCATCTTGTCCAg < 1:491720/147‑1 (MQ=255) gcTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTCCTCGGTCATCATCTTGTCCAGTTCGAcc > 1:170840/1‑150 (MQ=255) | AGGGTGGCCACCGCCGCACGCAGCTCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTCCTCGGTCATCATCTTGTCCAGTTCGACC > NC_002947/4922592‑4922880 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |