Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 1,844,118 | A→G | E342G (GAA→GGA) | sppA → | signal peptide peptidase SppA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,844,118 | 0 | A | G | 94.7% | 51.4 / ‑3.4 | 19 | E342G (GAA→GGA) | sppA | signal peptide peptidase SppA |
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (4/14); total (4/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 = 4.93e-01 |
ATTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCT > NZ_CP009273/1844044‑1844202 | aTTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGt < 1:134468/90‑1 (MQ=255) aaaCGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCg < 2:87375/90‑1 (MQ=255) aaaCGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCg < 2:201858/90‑1 (MQ=255) aaaCGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCg < 1:245475/90‑1 (MQ=255) gCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGa > 2:110303/1‑90 (MQ=255) gCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGa > 2:309698/1‑90 (MQ=255) gCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGATAcc < 1:158858/90‑1 (MQ=255) gCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATAcc < 1:93051/90‑1 (MQ=255) tACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGc < 2:309814/90‑1 (MQ=255) gTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGcac < 2:416837/90‑1 (MQ=255) aTTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGaaa < 1:521776/90‑1 (MQ=255) tggatggCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGa < 2:414951/90‑1 (MQ=255) ggatggCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGAt < 1:551262/90‑1 (MQ=255) gCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCtt < 2:105839/65‑1 (MQ=255) gCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCtt > 1:105839/1‑65 (MQ=255) gCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTcc > 2:473922/1‑90 (MQ=255) gCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTcc < 2:86464/90‑1 (MQ=255) cGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCt < 1:90953/90‑1 (MQ=255) cGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCt < 2:73689/90‑1 (MQ=255) | ATTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCT > NZ_CP009273/1844044‑1844202 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |