Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,547,989 | G→A | 12.9% | A371A (GCC→GCT) | rnr ← | exoribonuclease R |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,547,989 | 0 | G | A | 12.9% | 71.6 / 4.0 | 31 | A371A (GCC→GCT) | rnr | exoribonuclease R |
Reads supporting (aligned to +/- strand): ref base G (14/13); new base A (2/2); total (16/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01 |
TCCGGCACGACCTCCTTGTACTCTTCGCGCAGCGCCTTGCCCTCACGGGTACGCGAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCGGGGAAGTACACCGAGTTACCGCGCTGCTGGGCTTCGACATC > NC_002947/5547843‑5548138 | tCCGGCACGACCTCCTTGTACTCTTCGCGCAGCGCCTTGCCCTCACGGGTACGCGAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCt > 1:425110/1‑150 (MQ=255) ccGGCACGACCTCCTTGTACTCTTCGCTCAGCGCCTTGCCCTCACGGGTACGCGAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCtt < 2:65033/150‑1 (MQ=255) aCGACCTCCTTGTACTCATCGCGCAGCGCCTTGCCCTCACGGGTACGCGAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTc > 1:423053/1‑150 (MQ=255) ctTCGCGCAGCGCCTTGCCCTCACGGGTACGCGAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGa > 1:587704/1‑150 (MQ=255) cgcAGCGCCTTGCCCTCACGGGTACGCGAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATg < 2:719918/150‑1 (MQ=255) cTCACGGGTACGCGAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGAc < 1:683028/150‑1 (MQ=255) gcgAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATg > 2:186360/1‑107 (MQ=255) gcgAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGGGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATg < 1:186360/107‑1 (MQ=255) gAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGc < 1:408112/150‑1 (MQ=255) cTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGca < 2:423053/150‑1 (MQ=255) cTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATgg > 2:708181/1‑101 (MQ=255) cTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATgg < 1:708181/101‑1 (MQ=255) cTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCAGCTTTGTTCATGGTCATTTcaca > 1:791974/1‑93 (MQ=255) cTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCAGCTTTGTTCATGGTCATTTcaca < 2:792407/93‑1 (MQ=255) cTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCAGCTTTGTTCATGGTCATTTcaca < 2:791974/93‑1 (MQ=255) cTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCAGCTTTGTTCATGGTCATTTcaca > 1:792407/1‑93 (MQ=255) ttgttgTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGctct > 2:629663/1‑150 (MQ=255) gTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGc > 1:491714/1‑150 (MQ=255) ggCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGcac < 2:425110/150‑1 (MQ=255) ggAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCAcgcg > 2:792367/1‑150 (MQ=255) tGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCg < 2:587704/150‑1 (MQ=255) tCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAg < 2:302135/103‑1 (MQ=255) tCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAg > 1:302135/1‑103 (MQ=255) cGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCGGGGAAGt > 1:305812/1‑150 (MQ=255) gCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCGGGGAAGTa > 1:186630/1‑150 (MQ=255) gCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGATCGCTAAGTTGTTTCAGCGAGCACAGGCCGTTGGACACCTCTTCGGGGCGGATCGGCACCACGCGCTCGGGGAAGTa > 1:574698/1‑150 (MQ=255) ttCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGAc > 2:461612/1‑57 (MQ=255) ttCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGAc < 1:461612/57‑1 (MQ=255) gTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGTCATGCGGGCTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCGGGGAAGTACACCGa > 1:592281/1‑150 (MQ=255) cTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCGGGGAAGTACACCGAGTTACCGCGCTGCTGGGCTTc < 1:583790/150‑1 (MQ=255) cTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCGGGGAAGTACACCGAGTTACCGCGCTGCTGGGCTTc < 2:305812/150‑1 (MQ=255) ggCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCGGGGAAGTACACCGAGTTACCGCGCTGCTGGGCTTCGACATc > 1:737218/1‑150 (MQ=255) | TCCGGCACGACCTCCTTGTACTCTTCGCGCAGCGCCTTGCCCTCACGGGTACGCGAGTGCTCGAGCATGCTGCTGACCTTGTTGTAGGTCAGACGGGCATGGGAGTGGATCACGCCTTCGTAGAACTGATAGTCGACCATCTGGCCGGCTTTGTTCATGGTCATTTCACAGACCATGGCCAACCGGTCGACATGCGGGTTCAGCGAGCACAGGCCGTTGGACAGCTCTTCGGGGAGCATCGGCACCACGCGCTCGGGGAAGTACACCGAGTTACCGCGCTGCTGGGCTTCGACATC > NC_002947/5547843‑5548138 |
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: ‑ |