Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,903,567 | G→C | 17.9% | intergenic (+50/‑221) | PP_1703 → / → nirB | assimilatory nitrate reductase/sulfite reductase/nitrite reductase large subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,903,567 | 0 | G | C | 17.9% | 95.3 / 12.7 | 41 | intergenic (+50/‑221) | PP_1703/nirB | assimilatory nitrate reductase/sulfite reductase/nitrite reductase large subunit |
Reads supporting (aligned to +/- strand): ref base G (16/16); new base C (3/4); total (19/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTA > NC_002947/1903420‑1903715 | ggATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTgggcggg < 2:954619/149‑1 (MQ=255) gAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgccaatg < 1:1259337/150‑6 (MQ=255) aTGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgc > 1:875883/1‑141 (MQ=255) aTGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgc < 2:875883/141‑1 (MQ=255) gtcGGTTGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGa > 2:330313/3‑150 (MQ=255) atggatgCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCggg < 1:213103/150‑1 (MQ=255) tgCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgcc < 2:635808/129‑2 (MQ=255) tgCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgcc > 1:635808/1‑128 (MQ=255) gcggcggcgGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTAcc < 2:395511/150‑1 (MQ=255) ggctgcggctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgccaa > 1:1246186/1‑119 (MQ=255) ggctgcggctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgccaa < 2:1246186/122‑4 (MQ=255) gctgcggctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCgg < 1:763584/134‑1 (MQ=255) gctgcggctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCgg > 2:763584/1‑134 (MQ=255) gctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCggg < 1:1087602/129‑1 (MQ=255) gctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCggg > 2:1087602/1‑129 (MQ=255) gAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggcaac < 1:1282901/150‑1 (MQ=255) ttGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggcaacccaggc > 1:370860/1‑150 (MQ=255) gCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggcaacccaggcaacc > 2:1096183/1‑150 (MQ=255) cAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCgg > 2:344575/1‑108 (MQ=255) cAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCgg < 1:344575/108‑1 (MQ=255) aaGGCCAGGCGTTATCGGCGCGATGTGTACTTAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAggcaacccaggcaaccca < 2:265979/150‑1 (MQ=255) ccAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTTTCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCa > 1:1168669/1‑150 (MQ=255) cAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCaa > 1:408252/1‑150 (MQ=255) gTTAGCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGa > 2:548362/1‑150 (MQ=255) tATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGCCTGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGAtt < 1:737290/150‑1 (MQ=255) gcgcGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggcaacccaggcaa > 2:97850/1‑129 (MQ=255) gcgcGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggcaacccaggcaa < 1:97850/129‑1 (MQ=255) gtgtACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACTCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCt > 2:404814/1‑150 (MQ=255) gtACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACTCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGc < 1:404814/150‑1 (MQ=255) tACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCAc < 2:1335210/140‑1 (MQ=255) tACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCAc > 1:1335210/1‑140 (MQ=255) gCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCgg < 1:1283769/150‑1 (MQ=255) tCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGGGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcag > 1:904946/1‑95 (MQ=255) tCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGGGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcag < 2:904946/95‑1 (MQ=255) ccAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCcgc < 1:71645/78‑1 (MQ=255) ccAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCcgc > 2:71645/1‑78 (MQ=255) ggggTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGTAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACaaa < 1:330313/150‑1 (MQ=255) ggTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcag < 1:269664/85‑1 (MQ=255) ggTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcag > 2:269664/1‑85 (MQ=255) gttgtGCACCTTGAATGGGCCTTGGCTGGCCATGCCCAACGCAAGTGTAGGAGCGGGTTTACGCAATTGCGAGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCa < 1:1096183/150‑1 (MQ=255) gCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTCCCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGAt < 1:57931/150‑1 (MQ=255) tGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAGCCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGccc > 2:1109207/1‑150 (MQ=255) gcgggcCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggc < 1:544581/55‑1 (MQ=255) gcgggcCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggc > 2:544581/1‑55 (MQ=255) ggggcCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAGCCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTa < 1:1109207/149‑1 (MQ=255) cgggcCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCt > 1:282232/1‑93 (MQ=255) cgggcCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCt < 2:282232/93‑1 (MQ=255) | TGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTA > NC_002947/1903420‑1903715 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |