Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,440,194 | C→A | G279V (GGT→GTT) | rpoA ← | RNA polymerase subunit alpha |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,440,194 | 0 | C | A | 95.7% | 247.2 / 5.1 | 70 | G279V (GGT→GTT) | rpoA | RNA polymerase subunit alpha |
Reads supporting (aligned to +/- strand): ref base C (2/1); new base A (37/30); total (39/31) | |||||||||||
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.57e-01 |
GGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGA > NC_000913/3440147‑3440239 | ggTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAcc > 1:904501/1‑49 (MQ=255) ttAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGat < 1:4681622/50‑1 (MQ=255) ttAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGat < 1:5112676/50‑1 (MQ=255) ttAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGa > 1:4818165/1‑49 (MQ=255) tAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGata > 1:4349613/1‑50 (MQ=255) tAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata > 1:3814132/1‑50 (MQ=255) tAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata > 1:3931847/1‑50 (MQ=255) tAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata > 1:3707084/1‑50 (MQ=255) aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGata < 1:2037689/49‑1 (MQ=255) aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata < 1:615750/49‑1 (MQ=255) aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata < 1:1740767/49‑1 (MQ=255) aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata > 1:1595373/1‑49 (MQ=255) ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGatat < 1:4941528/49‑1 (MQ=255) ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGata > 1:3229384/1‑48 (MQ=255) gCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAg < 1:2934546/50‑1 (MQ=255) gCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAg < 1:2043796/50‑1 (MQ=255) gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTg > 1:1870220/1‑50 (MQ=255) gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTg < 1:4247227/50‑1 (MQ=255) ttttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTgg < 1:2282071/50‑1 (MQ=255) ttttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTgg < 1:4933949/50‑1 (MQ=255) ttttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTg > 1:4167682/1‑49 (MQ=255) tttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGa > 1:3947216/1‑50 (MQ=255) tttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGa > 1:4948198/1‑50 (MQ=255) tttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGa > 1:3845673/1‑50 (MQ=255) aaGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATa < 1:405531/49‑1 (MQ=255) aaGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAg > 1:190015/1‑50 (MQ=255) aGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGc > 1:4436282/1‑50 (MQ=255) aGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGc > 1:2126986/1‑50 (MQ=255) aGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGc < 1:3320206/50‑1 (MQ=255) ggAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCt > 1:2592147/1‑50 (MQ=255) ggAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCt > 1:1190483/1‑50 (MQ=255) aGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTc < 1:3094205/50‑1 (MQ=255) gCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCt < 1:2807718/50‑1 (MQ=255) cTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTg > 1:1774632/1‑50 (MQ=255) cTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTg > 1:3561291/1‑50 (MQ=255) cTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTg > 1:1624072/1‑50 (MQ=255) cTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTg > 1:2390179/1‑50 (MQ=255) aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt > 1:2960193/1‑50 (MQ=255) aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt > 1:793117/1‑50 (MQ=255) aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt > 1:1396223/1‑50 (MQ=255) aCCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCtt > 1:1921142/1‑49 (MQ=255) aCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCttt < 1:2088122/50‑1 (MQ=255) ccTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCttt > 1:168041/1‑49 (MQ=255) ccTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTa < 1:5011315/50‑1 (MQ=255) ccTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTa < 1:885595/50‑1 (MQ=255) cTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTaa < 1:1217759/50‑1 (MQ=255) ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg > 1:3792290/1‑49 (MQ=255) ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg < 1:2912575/49‑1 (MQ=255) ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg < 1:4690584/49‑1 (MQ=255) gTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg > 1:4753916/1‑48 (MQ=255) aCGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGt > 1:4279349/1‑50 (MQ=255) cGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGt < 1:3601995/49‑1 (MQ=255) cGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGt < 1:2851971/49‑1 (MQ=255) gCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGt > 1:4697761/1‑48 (MQ=255) cTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg > 1:2664602/1‑50 (MQ=255) cTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg < 1:2221849/50‑1 (MQ=255) cTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg < 1:173405/50‑1 (MQ=255) tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg > 1:894701/1‑47 (MQ=255) tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGc < 1:114677/48‑1 (MQ=255) tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGc < 1:4263760/48‑1 (MQ=255) tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGc < 1:4173665/48‑1 (MQ=255) aCCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCa > 1:2353051/1‑48 (MQ=255) ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag > 1:2590833/1‑50 (MQ=255) ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag > 1:255280/1‑50 (MQ=255) ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag > 1:1032281/1‑50 (MQ=255) cAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc > 1:284413/1‑50 (MQ=255) gATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCgg < 1:2156110/50‑1 (MQ=255) gATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCgg < 1:905751/50‑1 (MQ=255) aTCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa < 1:2812032/50‑1 (MQ=255) tCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa > 1:4411051/1‑49 (MQ=255) tCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa > 1:3647611/1‑49 (MQ=255) | GGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGA > NC_000913/3440147‑3440239 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |