| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | Exported | 834,575 | G→C | 26.1% | P1011P (CCG→CCC) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | Exported | 834,575 | 0 | G | C | 26.1% | 103.6 / 40.5 | 65 | P1011P (CCG→CCC) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base G (30/18); new base C (10/7); total (40/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.81e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01 | |||||||||||
GAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCT > Exported/834506‑834640 | gAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGa > 1:2512401/1‑71 (MQ=255) cacaaaGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCAc > 1:2355412/4‑71 (MQ=255) aaTAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCg < 1:1233782/69‑1 (MQ=255) aaTAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGa < 1:4425568/70‑1 (MQ=255) acaaaGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCAcc < 1:5475291/69‑1 (MQ=255) agaggaCGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctggt > 1:1276809/1‑71 (MQ=255) gaggaCGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctggtc < 1:2587385/68‑1 (MQ=255) gaagaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctg > 1:1353679/1‑68 (MQ=255) gaagaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACct > 1:5152866/1‑67 (MQ=255) aagaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtct > 1:3823411/1‑71 (MQ=255) aCGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctggtccggt < 1:4856072/71‑1 (MQ=255) aCGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctggtccggt < 1:1970290/71‑1 (MQ=255) aCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtctggt > 1:3877035/1‑71 (MQ=255) cGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCtggtctggtc < 1:3241924/71‑1 (MQ=255) cGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCtggtctggtc < 1:237431/71‑1 (MQ=255) gagaAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtctgg < 1:2699864/68‑1 (MQ=255) agaTCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGtctggtctgg < 1:1168819/69‑1 (MQ=255) gaTTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCAc > 1:2737495/1‑54 (MQ=255) gaTTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACCTGGTCCGGTCTGGaaa < 1:3602679/71‑2 (MQ=255) aTCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGtctcgtctgg < 1:728429/67‑1 (MQ=255) tCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGAt < 1:2317952/71‑1 (MQ=255) tcGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACCTGGTCCGGTCTGGaaag > 1:445746/2‑67 (MQ=255) tcGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACCTGGTCCGGTCTGGaaag > 1:4281304/2‑67 (MQ=255) tcGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACCTGGTCCGGTCTGGaaag > 1:4075744/2‑67 (MQ=255) cgcgATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACg < 1:1001419/70‑1 (MQ=255) cgcgATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGt > 1:4414918/1‑71 (MQ=255) cgcgATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGt > 1:4343050/1‑71 (MQ=255) gATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTc < 1:1515092/71‑1 (MQ=255) atatTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtc < 1:3785735/46‑1 (MQ=255) atatTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCt < 1:4720274/71‑1 (MQ=255) ttCAGGCACAGCCGCGCAAAATTATCTCCAGCCCg < 1:4831188/35‑1 (MQ=255) ttCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTa > 1:3980592/1‑70 (MQ=255) cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGAc > 1:1180945/1‑35 (MQ=255) cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCtggtctggtc > 1:562999/1‑46 (MQ=255) cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACaa > 1:913106/1‑71 (MQ=255) cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACa < 1:3276936/70‑1 (MQ=255) cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACa < 1:19916/70‑1 (MQ=255) aGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTggtctggtct > 1:4521034/1‑46 (MQ=255) aGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGgtctggtctg > 1:4440270/1‑47 (MQ=255) gCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGAt > 1:3828071/1‑51 (MQ=255) cacaGCCGCGCAAAATTATCTCCAGCCCGACctggtc > 1:548338/1‑37 (MQ=255) cacaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGcc > 1:5485810/1‑71 (MQ=255) gcaGCCGCGTAAAATTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGcgg < 1:1426495/70‑3 (MQ=255) acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCg > 1:775879/1‑71 (MQ=255) acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCg > 1:2810182/1‑71 (MQ=255) acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCg > 1:3431265/1‑71 (MQ=255) aGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGt < 1:3199465/71‑1 (MQ=255) aGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGt < 1:4438733/71‑1 (MQ=255) ccGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGaa > 1:1951803/1‑48 (MQ=255) ccGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGaa > 1:314518/1‑48 (MQ=255) gcAAAATTATCTCCAGCCCGCCCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTacac > 1:1316838/1‑70 (MQ=255) gcAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTacac > 1:4136311/1‑70 (MQ=255) taaaaTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACAcc > 1:2529333/2‑70 (MQ=255) cAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACa > 1:1576026/1‑56 (MQ=255) aaaaTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCggg > 1:2979045/1‑62 (MQ=255) aaaTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAAc > 1:692219/1‑71 (MQ=255) aaTTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACg > 1:5147618/1‑71 (MQ=255) ttATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGtt < 1:4566101/71‑1 (MQ=255) tctcCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCAc > 1:3109015/1‑71 (MQ=255) tctcCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCAc > 1:422689/1‑71 (MQ=255) cAGCCCGACCTGGTCTGGTCTGGAAGATGAACATGTTTCTTACAACGCCGGTTACACCAACGTTCAc < 1:3874263/67‑1 (MQ=255) cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAga < 1:3325402/71‑2 (MQ=255) cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGc > 1:2279472/1‑71 (MQ=255) cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGc > 1:3771642/1‑71 (MQ=255) cAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAgt > 1:3245122/1‑70 (MQ=255) aGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAgtt < 1:3154329/71‑3 (MQ=255) cGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGa > 1:5564691/1‑64 (MQ=255) | GAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCT > Exported/834506‑834640 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |