Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 1,592,137 | A→G | pseudogene (207/5422 nt) | BW25113_RS07915 ← | ECSE_1600 family autotransporter |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,592,137 | 0 | A | G | 89.8% | 125.9 / 8.4 | 49 | pseudogene (207/5422 nt) | BW25113_RS07915 | ECSE_1600 family autotransporter |
Reads supporting (aligned to +/- strand): ref base A (3/2); new base G (27/17); total (30/19) | |||||||||||
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.82e-01 |
CAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAATTTCGTTGTC > NZ_CP009273/1592052‑1592224 | cagccaAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCAt > 1:413072/1‑90 (MQ=255) aaaaTATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCa < 1:308321/90‑1 (MQ=255) aaaTATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCaa < 2:682547/90‑1 (MQ=255) cACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACt < 2:1248283/90‑1 (MQ=255) cTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTg > 1:1191018/1‑90 (MQ=255) cTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTg > 1:1289646/1‑90 (MQ=255) cTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTg > 2:1329021/1‑90 (MQ=255) cTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTg > 1:1417508/1‑90 (MQ=255) cTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTg > 2:1040886/1‑90 (MQ=255) gTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCa < 2:1346511/90‑1 (MQ=255) gCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGa > 1:324723/1‑90 (MQ=255) ccATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGAc < 2:849362/90‑1 (MQ=255) aTACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGc > 2:1201522/1‑90 (MQ=255) aTACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGc > 1:630803/1‑90 (MQ=255) aTACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGc > 1:817273/1‑90 (MQ=255) aTACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGc > 1:852657/1‑90 (MQ=255) aTACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGc > 2:909598/1‑90 (MQ=255) aTACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGc < 2:532113/90‑1 (MQ=255) aTACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGc > 1:453401/1‑90 (MQ=255) aCCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTg < 1:1307297/90‑1 (MQ=255) ccAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGa < 2:662626/90‑1 (MQ=255) aGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCaa > 2:1435986/1‑63 (MQ=255) aGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCaa < 1:1435986/63‑1 (MQ=255) aCCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAAcc < 1:70325/86‑1 (MQ=255) aCCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAAcc > 2:70325/1‑86 (MQ=255) gTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGc < 2:166438/90‑1 (MQ=255) tAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGcc < 2:26832/90‑1 (MQ=255) catcatAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCa > 1:396984/1‑90 (MQ=255) catAAGCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCaaaa > 1:696220/1‑90 (MQ=255) gCAACATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACAccc < 1:794306/90‑1 (MQ=255) aCATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGc > 1:838681/1‑90 (MQ=255) aCATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGc > 2:888380/1‑90 (MQ=255) aCATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGc > 2:45649/1‑90 (MQ=255) cATCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCg > 2:699348/1‑90 (MQ=255) aTCAGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGt < 1:543399/90‑1 (MQ=255) aGTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCaa < 1:1129840/90‑1 (MQ=255) gTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCa > 2:225707/1‑48 (MQ=255) gTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCa < 1:225707/48‑1 (MQ=255) gTATCAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAAt > 2:201409/1‑90 (MQ=255) cAATATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTAc < 1:485179/90‑1 (MQ=255) aaTATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACt > 2:323268/1‑90 (MQ=255) aaTATTGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACt > 1:1149860/1‑90 (MQ=255) tGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAAtt > 2:1504533/1‑90 (MQ=255) tGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAAtt > 2:1316811/1‑90 (MQ=255) tGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAAtt > 1:161369/1‑90 (MQ=255) tGGTAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAAtt > 2:959092/1‑90 (MQ=255) ggTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAAttt > 1:359127/1‑90 (MQ=255) tAATCTGGTCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAATTTCg < 1:1040886/90‑1 (MQ=255) tGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAATTTCGTTGTc < 2:1462041/90‑1 (MQ=255) | CAGCCAAAATATTAAGCACTCCAGTGCCATACCAGCCAACCAGGTAGGCATCATAAGCAACATCAGTATCAATATTGGTAATCTGATCATTATCAACTTCCAGACTTGCACCAGACGCTGAACCGCTTAGTGCCAGCAAAACACCCAGCGTCAATCTACTGAATTTCGTTGTC > NZ_CP009273/1592052‑1592224 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |