Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,158,080 | G→T | 12.1% | R19L (CGG→CTG) | opuA → | glycine betaine ABC transporter ATP‑binding protein/permease |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,158,080 | 0 | G | T | 12.1% | 87.2 / 4.1 | 33 | R19L (CGG→CTG) | opuA | glycine betaine ABC transporter ATP‑binding protein/permease |
Reads supporting (aligned to +/- strand): ref base G (15/14); new base T (2/2); total (17/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01 |
AAAACCCCATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGC‑CTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGTGGCGA > NC_002947/3157931‑3158221 | aaaaCCCCATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGC‑CTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGcgcg > 1:608123/1‑150 (MQ=255) aCCCCATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGC‑CTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGAt < 1:1137043/150‑1 (MQ=255) cccATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGC‑CTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCTGATcc > 2:685987/1‑149 (MQ=255) ccATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGCCGCTCTTC‑ACACTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCGGTTAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAg > 1:969483/1‑150 (MQ=255) cATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGC‑CTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAgc < 2:1094554/150‑1 (MQ=255) gcgcTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGC‑CTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCTGATCCGCGCAgtg < 2:440381/150‑3 (MQ=255) gcgcTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGC‑CTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCTGATCCGCGCAgtg < 1:685987/150‑3 (MQ=255) ttGATTACGCGGGTTTTGGCCGATTGGCCGGGAGC‑CTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGatcgga > 1:37162/1‑145 (MQ=255) ggCCGATTGGCATGGAGC‑CTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGAtt > 1:1407050/1‑150 (MQ=255) aTGGAGC‑CTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCGTGTCCCTTTATCGATTCACTATATTgc > 2:1374592/1‑150 (MQ=255) gAGC‑CTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCt < 1:1374592/150‑1 (MQ=255) tCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCg < 1:620066/150‑1 (MQ=255) cTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTTGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGc < 1:1139354/150‑1 (MQ=255) cTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCa > 1:437714/1‑150 (MQ=255) ttGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTTCGGTCACCAACCGGTa > 2:169041/1‑150 (MQ=255) cGTGATCGCTGATAAATACCCACAGATAACCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGa > 1:398554/1‑150 (MQ=255) tCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTcca < 2:608123/150‑1 (MQ=255) gATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAAc > 2:484278/1‑150 (MQ=255) cccACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAACCGCCCTTGGCTCTGGATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGa < 1:169041/150‑1 (MQ=255) ccACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAg < 2:1407050/150‑1 (MQ=255) cacaGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTg > 1:1317353/1‑57 (MQ=255) cacaGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTg < 2:1317353/57‑1 (MQ=255) ccAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAAcc > 2:1049270/1‑150 (MQ=255) cAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCgg > 1:857966/1‑115 (MQ=255) cAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCgg < 2:857966/115‑1 (MQ=255) cAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCt > 1:1239976/1‑150 (MQ=255) aaaCCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGTg < 1:1114533/150‑1 (MQ=255) gtcTGCGCTGATCCGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGTGGCGa > 1:996822/3‑150 (MQ=255) tGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTcgcg < 2:429452/133‑1 (MQ=255) tGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTcgcg < 2:429774/133‑1 (MQ=255) tGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTcgcg > 1:429774/1‑133 (MQ=255) tGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTcgcg > 1:429452/1‑133 (MQ=255) tGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGt > 2:1037135/1‑135 (MQ=255) tGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGt < 1:1037135/135‑1 (MQ=255) | AAAACCCCATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGC‑CTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGTGGCGA > NC_002947/3157931‑3158221 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |