Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,237,723:1 | +G | 10.5% | coding (904/1551 nt) | PP_1975 ← | EmrB/QacA family drug resistance transporter |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,237,723 | 1 | . | G | 10.5% | 133.5 / 11.3 | 38 | coding (904/1551 nt) | PP_1975 | EmrB/QacA family drug resistance transporter |
Reads supporting (aligned to +/- strand): ref base . (18/16); new base G (2/2); total (20/18) | |||||||||||
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 |
GCAAAACCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCTTGGAAAACTGCAGGATGACGAA > NC_002947/2237583‑2237867 | gCAAAACCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGGCG‑CAGGTACAg < 2:212715/150‑1 (MQ=255) cccAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGc > 1:449463/1‑150 (MQ=255) ccAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTG‑GCGCAGGTACAg > 2:751597/1‑142 (MQ=255) ccAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTG‑GCGCAGGTACAg < 1:751597/142‑1 (MQ=255) aGCGCGCATAGTACCTTGGGCGATATCACCTTCACCAGTTTGGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCaga < 1:225201/150‑1 (MQ=255) gcgcgcATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCagat < 1:502909/150‑1 (MQ=255) cATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATc < 1:675385/150‑1 (MQ=255) cTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTGACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTa > 2:744985/1‑150 (MQ=255) cTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTa < 1:727854/150‑1 (MQ=255) aCCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTa > 2:81519/1‑118 (MQ=255) aCCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTa < 1:81519/118‑1 (MQ=255) ccTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCa > 1:848520/1‑150 (MQ=255) gCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTgg < 1:236255/150‑1 (MQ=255) aGCGGGATCAGGACGAGCTACGGAACACGCATCCACGTGATCACCTCGCCAATCAGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTc < 2:398542/150‑1 (MQ=255) ggATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTGACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACa < 1:744985/150‑1 (MQ=255) aTCAGGATTAGCCGCGGAATACCCATCCACAGGAGCGCGTCGCTAATGTGCAAGGCGTCGTAACCCCGCAGCTACGC‑CAGGTACAGCTGTAGCAGATAGTTGGAGCTGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACatg < 2:299982/150‑3 (MQ=255) cGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGAGGtt > 2:308203/1‑150 (MQ=255) tACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGa < 1:560401/150‑1 (MQ=255) ccATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCca > 2:205472/1‑96 (MQ=255) ccATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCca < 1:205472/96‑1 (MQ=255) ccATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGc < 2:631629/105‑1 (MQ=255) ccATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGc > 1:631629/1‑105 (MQ=255) ccATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCaaa > 2:4769/1‑133 (MQ=255) ccATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCaaa < 1:4769/133‑1 (MQ=255) aTCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATa > 2:86478/1‑150 (MQ=255) cTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAg > 1:587230/1‑150 (MQ=255) aaTCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTg > 1:655692/1‑150 (MQ=255) aTCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGa < 1:86478/150‑1 (MQ=255) aTCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGa < 2:655692/150‑1 (MQ=255) tCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAg > 2:172535/1‑150 (MQ=255) ggCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCTTGGa > 1:395992/1‑150 (MQ=255) ttGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCTTGGAAAAc > 1:148207/1‑150 (MQ=255) tGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCTTGGAAAACt > 2:223603/1‑150 (MQ=255) ccTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTAGAAGCCCACACCAATGCCCAGGCTGGCGATGCTCGACAGGAAAAAGGTGCGGTTGCCGAGGATACGCAGGTTGATCAGCTGGTCAGCCTTGGAAAACTGCAggg > 2:155660/1‑149 (MQ=255) cTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCTTGGAAAACTGCAGGAt > 2:364691/1‑150 (MQ=255) cTGTACCTG‑GCGCAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACa > 2:812361/1‑83 (MQ=255) cTGTACCTG‑GCGCAGGTACAGCGGTAGCAGATAGATCGAGCCCTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACa < 1:812361/83‑1 (MQ=255) cTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCTTGGAAAACTGCAGGATGACGaa > 2:373737/1‑150 (MQ=255) | GCAAAACCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAGGAACAGCTGCGGAATACCCATCCACATGATCACCTCGCCAATCTGCAAGGCGTTGTAACCCTGTACCTGCGC‑CAGGTACAGCGGTAGCAGATAGATCGAGCCGTACAGGCCCACACCCATGCCCAGGCTGGCGATGCTCGACAGGCCAAAGTTGCGGTTGCCGAGGATACGCAGGTTGATCAGCGGGTGAGGCTTGGAAAACTGCAGGATGACGAA > NC_002947/2237583‑2237867 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |