Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 2,882,413 | A→C | V349G (GTT→GGT) | cysI ← | assimilatory sulfite reductase (NADPH) hemoprotein subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,882,413 | 0 | A | C | 98.1% | 170.9 / ‑4.2 | 52 | V349G (GTT→GGT) | cysI | assimilatory sulfite reductase (NADPH) hemoprotein subunit |
Reads supporting (aligned to +/- strand): ref base A (0/1); new base C (26/25); total (26/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01 |
CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGC > NZ_CP009273/2882326‑2882498 | ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc > 1:1106299/1‑90 (MQ=255) ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc < 2:1392305/90‑1 (MQ=255) ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc > 1:311248/1‑90 (MQ=255) ggTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCa < 1:446811/90‑1 (MQ=255) tCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCaa < 1:1326323/90‑1 (MQ=255) cAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAAt > 1:916720/1‑90 (MQ=255) gCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGGGCCAGTTATCATCAATGCCCTTACCCCAGCCAATAc > 1:744189/1‑90 (MQ=255) ggACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGa > 2:1473969/1‑90 (MQ=255) ggACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGa > 2:197782/1‑90 (MQ=255) cgcCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGcc < 2:540901/90‑1 (MQ=255) cgcCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGcc < 2:1443466/90‑1 (MQ=255) aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc < 2:1245599/90‑1 (MQ=255) aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc < 1:89263/90‑1 (MQ=255) aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc < 1:468380/90‑1 (MQ=255) aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc < 2:188291/90‑1 (MQ=255) aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc < 2:267520/90‑1 (MQ=255) aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc < 2:509360/90‑1 (MQ=255) aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc < 2:618901/90‑1 (MQ=255) aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc < 2:916720/90‑1 (MQ=255) aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAAc < 1:1048928/90‑1 (MQ=255) ggATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAGTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc < 2:173541/89‑1 (MQ=255) ggATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAGTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc > 1:173541/1‑89 (MQ=255) ggATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACt > 1:1158862/1‑90 (MQ=255) aTGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCa > 1:922861/1‑90 (MQ=255) aTGCGACCATTTTCGATAAACAGCGTCAGGGGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCa > 1:1445727/1‑90 (MQ=255) gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc > 1:1062458/1‑55 (MQ=255) gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc < 2:1062458/55‑1 (MQ=255) gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCata < 2:1333493/90‑1 (MQ=255) gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCata > 1:906744/1‑90 (MQ=255) gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCata > 1:634879/1‑90 (MQ=255) cGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCatat > 1:284171/1‑90 (MQ=255) gACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATg > 2:162190/1‑90 (MQ=255) tCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTAACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATc < 2:1410061/90‑1 (MQ=255) aaaCAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTc < 2:532449/90‑1 (MQ=255) aaaCAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTc < 2:284171/90‑1 (MQ=255) gCGTCAGGTGCCAGTTATCATCAATGCCCTTAcccc < 2:210131/36‑1 (MQ=37) gCGTCAGGTGCCAGTTATCATCAATGCCCTTAcccc > 1:210131/1‑36 (MQ=37) gCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAg > 2:1420539/1‑38 (MQ=38) gCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAg < 1:1420539/38‑1 (MQ=38) gCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAAtt > 2:871746/1‑90 (MQ=255) gCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAAtt > 2:262277/1‑90 (MQ=255) cGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAAttt > 2:672618/1‑90 (MQ=255) cAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGAt > 2:231294/1‑90 (MQ=255) aGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATc > 2:946484/1‑90 (MQ=255) ggTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATcc > 2:911396/1‑90 (MQ=255) cAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATc > 2:475029/1‑69 (MQ=255) cAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATc < 1:475029/69‑1 (MQ=255) aGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCcgcgc < 2:341467/90‑1 (MQ=255) aGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCcgcgc < 1:911396/90‑1 (MQ=255) aGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCcgcgc < 1:112213/90‑1 (MQ=255) catcaATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCt < 1:271295/90‑1 (MQ=255) catcaATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCt < 1:1210931/90‑1 (MQ=255) aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt > 2:241418/1‑90 (MQ=255) aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt > 2:416263/1‑90 (MQ=255) aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt > 2:503816/1‑90 (MQ=255) aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt > 1:372300/1‑90 (MQ=255) aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt > 2:1070749/1‑90 (MQ=255) cccTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTcc < 2:1445727/90‑1 (MQ=255) cTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGc > 1:1443269/1‑90 (MQ=255) | CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGC > NZ_CP009273/2882326‑2882498 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |