Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 1,568,198 | A→G | F417S (TTC→TCC) | pqqL ← | insulinase family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,568,198 | 0 | A | G | 97.9% | 149.8 / ‑3.8 | 48 | F417S (TTC→TCC) | pqqL | insulinase family protein |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (17/30); total (18/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.75e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01 |
AACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAAC > NZ_CP009273/1568111‑1568274 | aaCTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGAc < 1:747037/90‑1 (MQ=255) aCTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACg < 2:627936/90‑1 (MQ=255) ctgctgCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg < 1:107845/90‑1 (MQ=255) ctgctgCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg < 2:531752/90‑1 (MQ=255) ctgctgCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg < 2:1301776/90‑1 (MQ=255) ctgctgCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg < 1:1364015/90‑1 (MQ=255) ctgctgCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg < 1:475266/90‑1 (MQ=255) gctgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGGCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCg > 1:1077087/1‑90 (MQ=255) tgcCATTTTCCCGCCAGCCTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGta < 1:779406/90‑1 (MQ=255) cATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGtatta < 1:301946/90‑1 (MQ=255) cAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCt < 2:1157367/90‑1 (MQ=255) gACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTgg < 2:1181878/90‑1 (MQ=255) cTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGAt < 2:497420/90‑1 (MQ=255) ttGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg > 1:881977/1‑66 (MQ=255) ttGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg < 2:881977/66‑1 (MQ=255) cACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAg < 2:85490/90‑1 (MQ=255) tAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGAc < 2:363476/90‑1 (MQ=255) tAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGAc < 1:1066006/90‑1 (MQ=255) aTTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTg > 1:1431059/1‑90 (MQ=255) ttGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGt > 1:1331905/1‑90 (MQ=255) tgctgcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTc < 1:557051/90‑1 (MQ=255) tgctgcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTc < 2:622565/90‑1 (MQ=255) ctgcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGAc < 1:1145160/83‑1 (MQ=255) ctgcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGAc > 2:1145160/1‑83 (MQ=255) gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg > 1:189022/1‑49 (MQ=255) gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg < 2:189022/49‑1 (MQ=255) gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTaa < 2:919414/61‑1 (MQ=255) gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTaa > 1:919414/1‑61 (MQ=255) gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCa > 2:769334/1‑90 (MQ=255) ccACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCAt < 1:769334/90‑1 (MQ=255) caGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCa > 1:438080/1‑86 (MQ=255) caGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCa < 2:438080/86‑1 (MQ=255) tttCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATc < 1:1405584/90‑1 (MQ=255) tttCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATc > 2:1155111/1‑90 (MQ=255) aGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCa < 2:12570/90‑1 (MQ=255) gTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAg < 1:591924/90‑1 (MQ=255) tctcTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTa < 2:422727/69‑1 (MQ=255) tctcTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTa > 1:422727/1‑69 (MQ=255) ctctTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGAc > 2:1319985/1‑53 (MQ=255) ctctTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGAc < 1:1319985/53‑1 (MQ=255) ctctTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATc > 1:1454881/1‑90 (MQ=255) ctctTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATc > 2:825532/1‑90 (MQ=255) ctctTCCGGCGACAAGGACGGCGTACTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATc > 1:1118483/1‑90 (MQ=255) ctctTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATc > 1:570127/1‑90 (MQ=255) ctTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCaa > 1:732172/1‑90 (MQ=255) ttCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAAc > 2:50645/1‑90 (MQ=255) ttCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAAc < 1:1229957/90‑1 (MQ=255) ttCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAAc < 2:732172/90‑1 (MQ=255) | AACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAAC > NZ_CP009273/1568111‑1568274 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |