| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | AP009048 | 4,017,701 | T→C | 12.9% | N1175D (AAC→GAC) | rhsB ← | rhsB element core protein RshB |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | AP009048 | 4,017,701 | 0 | T | C | 12.9% | 158.9 / 14.2 | 62 | N1175D (AAC→GAC) | rhsB | rhsB element core protein RshB |
| Reads supporting (aligned to +/- strand): ref base T (32/22); new base C (3/5); total (35/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.79e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CATCATACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGC > AP009048/4017632‑4017758 | catcatACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTg > 1:4319726/1‑71 (MQ=34) atcatACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGtt < 1:3480034/69‑1 (MQ=35) atcatACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGc > 1:3627337/1‑71 (MQ=34) atcatACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGc < 1:1752779/71‑1 (MQ=34) atcatACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGc < 1:4405417/71‑1 (MQ=34) atACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGccc > 1:4517920/1‑70 (MQ=35) tACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCa > 1:2578004/1‑71 (MQ=34) ctgcCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGccc > 1:3718923/1‑64 (MQ=35) ctgcCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtcat > 1:3985109/1‑71 (MQ=34) ctgcCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtca < 1:4087639/70‑1 (MQ=35) ctgcCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCACtc < 1:2328411/69‑1 (MQ=25) tgcCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtcat > 1:3631792/1‑70 (MQ=35) gcCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtcatca > 1:678614/1‑71 (MQ=35) ccccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtcatca > 1:870433/1‑70 (MQ=37) ccccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtca < 1:3976256/67‑1 (MQ=35) cccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCata < 1:4454545/71‑1 (MQ=39) ccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCata > 1:1092352/1‑70 (MQ=255) gCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGcccc < 1:2886092/57‑1 (MQ=255) ccGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTgcg > 1:3466368/1‑70 (MQ=37) tAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcaccac > 1:1792733/1‑71 (MQ=35) tAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcaccac > 1:4471175/1‑71 (MQ=35) tAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcaccac > 1:7008/1‑71 (MQ=35) aaGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACg > 1:1515996/1‑71 (MQ=38) gctgctGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCt > 1:1029605/1‑71 (MQ=40) gctgctGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCt < 1:3989112/71‑1 (MQ=40) gctgctGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCt < 1:1847099/71‑1 (MQ=255) ctgcttCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTg < 1:2567677/71‑1 (MQ=255) tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcacca > 1:786884/1‑65 (MQ=35) gctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGtt > 1:195973/1‑71 (MQ=255) ctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCTTCATATTCTGCGCACCACGCTGt < 1:42191/69‑1 (MQ=255) tgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGc > 1:3718220/1‑71 (MQ=255) gcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCt > 1:1176722/1‑53 (MQ=35) gcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcaccac < 1:455822/62‑1 (MQ=37) cTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCt < 1:4036418/61‑1 (MQ=255) cTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCtt > 1:3893048/1‑71 (MQ=40) tGATGCGGGTTCTCTTCATTCAGCAGGTTG‑‑CCATTCATCATATTc > 1:2079367/1‑45 (MQ=11) gATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtcatca > 1:4422842/1‑41 (MQ=32) gATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTcc > 1:3352164/1‑71 (MQ=39) aTGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCg < 1:1778066/71‑1 (MQ=39) tGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTcc < 1:704746/69‑1 (MQ=40) tGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGt < 1:1970642/71‑1 (MQ=39) gCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGccc < 1:1983706/63‑1 (MQ=40) cGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGc > 1:1585624/1‑71 (MQ=39) gggTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCt > 1:1924373/1‑42 (MQ=38) ggTTCTCTTCATTCAGCAGGTCGCCCCATTCATCata < 1:3801000/37‑1 (MQ=255) ggTTCTCTTCATTCAGCAGGTCGCCCCATTCATCata < 1:1766175/37‑1 (MQ=255) gTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATAtt < 1:609393/38‑1 (MQ=40) gTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGt > 1:210751/1‑54 (MQ=255) gTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGt > 1:205019/1‑54 (MQ=255) ttCTCTTCATTCAGCAGGTTGCCCCATTCATCata < 1:3240154/35‑1 (MQ=40) ttCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcacca < 1:3888224/47‑1 (MQ=38) tctcTTCATTCAGCAGGTTGCCCCATTCATCATATTCTg > 1:3102495/1‑39 (MQ=40) tctcTTCATTCAGCAGGTTGCCCCATTCATCATATTCTg > 1:2718997/1‑39 (MQ=40) ctctTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCtgatg < 1:3704978/70‑1 (MQ=40) ctctTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCtgat > 1:252781/1‑69 (MQ=40) ttcattcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCtgatga < 1:563/68‑1 (MQ=37) ttcattcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAgcg > 1:1505312/1‑71 (MQ=35) tcattcaGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGc > 1:1334903/1‑44 (MQ=255) cattcaGCAGGTTGCCCCATTCCTCATATTCTGCGca > 1:1014866/1‑37 (MQ=255) cattcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAgcg < 1:3246588/69‑1 (MQ=35) cattcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAgc < 1:2265494/68‑1 (MQ=35) cattcaGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAgcgc > 1:157131/1‑70 (MQ=25) cattcaGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAgcgc < 1:1142593/70‑1 (MQ=25) | CATCATACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGC > AP009048/4017632‑4017758 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |