Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,782,577 | C→A | 6.9% | Q514H (CAG→CAT) ‡ | glnD ← | uridylyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,782,577 | 0 | C | A | 6.9% | 178.1 / 2.8 | 58 | Q514H (CAG→CAT) ‡ | glnD | uridylyltransferase |
Reads supporting (aligned to +/- strand): ref base C (25/29); new base A (2/2); total (27/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 = 5.63e-01 |
GTCGGACAAGTCCTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTTGATGAGGTTGAGG > NC_002947/1782430‑1782724 | gTCGGACAAGTCCTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTg < 2:509979/150‑1 (MQ=255) gTCGGACAAGTCCTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTg < 1:415598/150‑1 (MQ=255) tCGGACACGTCCTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATC‑GGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCGGAACTTTTGCGCGTCCAGTGCACCGATTTCCGAGTGGTCGCCCTGGc < 1:558576/150‑1 (MQ=255) gACAAGTCCTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGa > 2:558498/1‑150 (MQ=255) gACAAGTCCTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGa > 1:522754/1‑150 (MQ=255) aCAAGTCCTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGAc > 1:569356/1‑150 (MQ=255) aGTCCTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCAGGGCGAccac < 2:529724/150‑3 (MQ=255) ggccTGGCGCTGAGCGGTGGTCGACATCACCATGGGTGGCTGCACCAGCCAGACGATCAGGCGGCTGGCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCtt > 2:250158/3‑150 (MQ=255) tCCTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTg > 1:91087/1‑150 (MQ=255) ttGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCAGGGCGAc > 1:118180/1‑142 (MQ=255) ttGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCAGGGCGAc < 2:118180/142‑1 (MQ=255) cTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTc < 1:329395/150‑1 (MQ=255) cGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGt > 2:204764/1‑150 (MQ=255) cGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGCCCCTTGCCGATGTCGTGGTACAGccc < 1:204764/150‑1 (MQ=255) cACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAg < 1:768194/150‑1 (MQ=255) aCCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAgg > 1:595484/1‑150 (MQ=255) cAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGc > 1:620810/1‑122 (MQ=255) cAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGc < 2:620810/122‑1 (MQ=255) gtggtTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGAt < 2:630671/150‑1 (MQ=255) ttCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGc < 2:569356/150‑1 (MQ=255) aCCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTg < 1:372292/119‑1 (MQ=255) aCCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTg > 2:372292/1‑119 (MQ=255) cagccagACGATCAGTCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTgg < 1:317017/97‑1 (MQ=255) cagccagACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTgg > 2:317017/1‑97 (MQ=255) gccagACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTgg > 1:364628/1‑150 (MQ=255) agACGCTCAGGCGGCTGTCACAGGCCGGCAACTGGTTCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGTTTTGGGCa > 2:480530/1‑150 (MQ=255) agACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGc > 1:764961/1‑149 (MQ=255) agACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGc < 2:764961/149‑1 (MQ=255) agACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCa > 2:400177/1‑150 (MQ=255) cAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGCCCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGccc < 2:544539/150‑1 (MQ=255) gcggcTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCTAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCAt < 2:604119/150‑1 (MQ=255) gcggcTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCAt > 2:849548/1‑150 (MQ=255) tGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGAt < 1:113706/118‑1 (MQ=255) tGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGAt > 2:113706/1‑118 (MQ=255) cccAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGc < 1:558498/150‑1 (MQ=255) ccAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGGTTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCt < 1:262827/150‑1 (MQ=255) ccAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTa > 1:20522/1‑111 (MQ=255) ccAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTa < 2:20522/111‑1 (MQ=255) cAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTCCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTg < 1:250158/150‑1 (MQ=255) gCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCg > 1:646489/1‑99 (MQ=255) gCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCg < 2:646489/99‑1 (MQ=255) gCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCtt > 1:818727/1‑143 (MQ=255) gCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCtt < 2:818727/143‑1 (MQ=255) cTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGaa < 1:264901/150‑1 (MQ=255) ttCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGa > 1:835717/1‑150 (MQ=255) cAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCg > 2:625334/1‑144 (MQ=255) cAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCg < 1:625334/144‑1 (MQ=255) ttCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCgtg > 2:140128/1‑150 (MQ=255) gTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTg < 2:254158/54‑1 (MQ=255) gTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTg > 1:254158/1‑54 (MQ=255) gTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAg > 1:130572/1‑150 (MQ=255) aCTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGTCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTg > 1:696861/1‑150 (MQ=255) aTTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGc < 1:140128/150‑1 (MQ=255) gTGGTCGCCCTGGCGACCCTTGCGGATGTCTTCGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTTGATGAg < 2:91087/150‑1 (MQ=255) gTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTTGATGAg < 2:835717/150‑1 (MQ=255) agggTCGCCAGGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTTGATGAg > 1:625545/3‑150 (MQ=255) cGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTTGATGAGGTTGa < 2:559309/150‑1 (MQ=255) gCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTTGATGAGGTTGAg < 1:400177/150‑1 (MQ=255) cccTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTTGATGAGGTTGAgg < 2:482057/150‑1 (MQ=255) | GTCGGACAAGTCCTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTTGATGAGGTTGAGG > NC_002947/1782430‑1782724 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |