| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 16,749 | T→A | intergenic (+192/+2) | insL1 → / ← hokC | IS186/IS421 family transposase/protein HokC |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 16,749 | 0 | T | A | 92.6% | 70.2 / ‑3.0 | 27 | intergenic (+192/+2) | insL1/hokC | IS186/IS421 family transposase/protein HokC |
| Reads supporting (aligned to +/- strand): ref base T (1/1); new base A (23/2); total (24/3) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.14e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAAC‑GCGAGGCATC‑GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCA‑GATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACA‑GATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAG‑CATATAGGAGGCCTCGGGTTGATG‑GTAAAATATCACTC > NC_000913/16506‑16995 | gCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAAC‑GCGAGGCATC‑GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCA‑GATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTc > 1:88519/1‑251 (MQ=38) cATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAAC‑GCGAGGCATC‑GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCA‑GATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCgg > 1:350047/1‑251 (MQ=39) ccAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAAC‑GCGAGGCATC‑GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCA‑GATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGaaa > 1:76393/1‑250 (MQ=255) gAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTAGAGGTGGCCATTCGTGGGACAGTATCCCTGAAAGCCTACAAAACGCAATTGAAGAAC‑GCGAGGCATC‑GTCTTAACGAGTCACCTAGGCGTCGCATTCTTCA‑GATAGTTCAACCCTTAAGATAGCGCTTATGGGATCACTCCCCGCCAGTGCACTTACTCGGATTCGAAAACCGTGAAAAc > 1:268843/1‑249 (MQ=14) aacagaaaagaaaaATTACCTCGTTGGGGAAAAAAAAAAAAATGGTCATCTGGATTTTACGGGGGGCCATCGGGGGGACAGTCCCCCTAACCCCCTAAAAAACGAAATTGAAGACC‑GGGGGGCATC‑GTCTTACCGAGGCCCGGGGGTGCCGCTTTCTTAA‑GAGGGTTCACCCTTTAATTGAGCGTTTAGGGGACCCCCCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGt < 2:186276/247‑1 (MQ=255) accgaaaaaaagaattacccggTTGTGGAGAAAAAAAAAAATGGTCATGGGGAGATAAGAGGTGGCCATGCTGGGGAAATTACCCCTGACACCCTAAAAAACGCAAT‑GAAAAACCGGGAGG‑ATCCGTTTAAACGAGCCACCGAGGCGTCGC‑TTCCTCAAGTTGGTTCACCCGTTAATTTAGCGCTTAGGGAATCCCGCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTc < 2:204826/230‑1 (MQ=255) gaaCTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAAC‑GCGAGGCATC‑GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCA‑GATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCAAGTTCg > 1:315133/1‑251 (MQ=255) tCTGGAGCTTACAGGTGGCCATTCGTGGGGCAGTTTCCCTGACAGCCTACTAAACGCAATTGTAGAAC‑GCGAGGCATC‑GTCTTACCGCGGCAACGAGGAGTCGCATTCTCCA‑GATGGTTCAACCATTAAGTTAGCGCTTATGGAATCACTCCCCGCCGATGCACTTACTCGGATTCGTAAGCCGTGAAAAAAGCAACCTCCGTCTGGCCAGTTCGGTTGTGAAACTCACAGAGGTCTTTTCTCGTTACCac > 2:385463/1‑250 (MQ=21) aCAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAAC‑GCGAGGCATC‑GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCA‑GATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCtttt > 2:335252/1‑230 (MQ=255) cAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAAC‑GCGAGGCATC‑GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCA‑GATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTCCTCGGATTCGTAAGCCGTGTAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAGCCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACt > 1:342386/1‑250 (MQ=255) gggACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAAC‑GCGAGGCATC‑GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCA‑GATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATAc > 1:331010/1‑248 (MQ=255) ggACAGTATCCCTGACAGCCTACAAAACGCAATTGCAGACC‑GCGAGGCATC‑GTCTTACCGAGGCACCGAGGCGTCGCATTCTTCA‑GATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCTGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGCGGTCTTTTCTCGTTACCAGCGCCGCCACTAcggcggggggacaact > 2:277614/1‑241 (MQ=25) cTGACAGCCTACAAAACGCAATTGAAGAAC‑GCGAGGCATC‑GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCA‑GATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACA‑GATGACGATCAggg > 1:55315/1‑251 (MQ=255) tGAAGAAC‑GCGAGGCATC‑GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCA‑GATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATAAC‑GATGACGGTCAGGGCGACCATCCTCGCCCTAtgctg > 2:215372/1‑251 (MQ=37) aagaaC‑GCGAGGCATC‑GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCA‑GATGGTTCAACCCTTACGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGAGATCCAAGAT‑ACGCTCAGGGCGGCAATTATCCCCCTAtgatgcc > 2:466798/1‑250 (MQ=37) gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCA‑GATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGGGATACA‑GATGACGATCAGGGGGGCAAACATCGCCTTTTGCTGCTTCATTTCTCTCTtctc > 2:65341/1‑251 (MQ=38) aaCGAGGCACCGAGGCGTCGCATTCTTCA‑GATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACA‑GATGACGATCAGGGCGCCAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTgg > 1:378184/1‑250 (MQ=255) cGAGGCGTCGCATTCTTAA‑GATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGGTACA‑GATTACGATCAGGGCGCCAATCATCCCCTTTTGCTGCTTTCTTGTTCTCTTTCCCTTTCCCTTACGGTc > 2:229718/1‑251 (MQ=255) gAGGCGTCGCATTCTTCA‑GATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGGTAAA‑GATGACGATCAGGGCGGCAAACATCGCCTTATGCTGCCTCATTTCTCTCTTCTCCtttccccttcggtc > 2:451976/1‑238 (MQ=255) aGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATTCA‑GATGACGATTAGGGCGGCAATCATCGCCTTTTTATTCTTTATTGCTCTCTTCTCCTTTGTCTTTCGGTCAGTAAAAGGCAATCTACacttgtttaacttata > 2:320759/1‑234 (MQ=255) aTGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACA‑GATTACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAG‑CATATAGGGGGCCTCGGGt > 1:179370/1‑251 (MQ=255) gATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACA‑GATTACGATCCAGGCGACAATCAACGCCTTATGCTTCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAATAAGAGGCACTCTCCATGTGTTCAGACA‑ATAGGGGGCCCCGGGGTGAt > 1:99925/1‑251 (MQ=255) aCTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATAAA‑GATGAAGATCAGGGGGGAAATCAACGCCTTATGCTGCTTCAATGCTTTCTTCTCCTTGACCCTACGGTCCGTAAGAGGCACTCTAAATGTGTTCAG‑CAAATAGGGGGCCCAGGGTTGAtagt > 1:323244/1‑247 (MQ=255) cgccgATGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGGTTCA‑AGTGAAGATCAGGGCGGAAATCATCGCCTTATGCCGCTTCCTTGCTCTCTTCTCCTTTACCTTTCGTGCCGTAAAAGGCACTCCACATTTTTTCCG‑CATTTAGGGGGCCTCGGGGTGATA‑ATAAAata > 2:449098/1‑250 (MQ=255) cgTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACC‑AATGACGATCAAGGCGACAATCATCGCCCTATGCTTCTTCATTGCCCTCTTCTCCCTGACCTTTCCGTCAGTAAGCGGCACTCCACATGTGTTCAG‑CCTATAGGGGGCCTCGGGTTGATGAGTAAA‑TATCAc > 1:459687/1‑251 (MQ=255) gTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACA‑GATGGCGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCCGTAAGAGGCACTCTACATGTGTTCAG‑CATATAGGGGGCCTCGGGGTGATA‑GTAAAATATCAcc > 2:11312/1‑250 (MQ=255) tGCTCTTATCCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGGACCTCACAGAGGTCTTTTCTCGTTCCCAGCGCCGCCACTACGGCGGTGATACA‑GATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAG‑CATATAGGAGGCCTCGGGTTGATG‑GTAAAATATCACTc < 1:124858/250‑1 (MQ=255) | GCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAAC‑GCGAGGCATC‑GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCA‑GATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACA‑GATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAG‑CATATAGGAGGCCTCGGGTTGATG‑GTAAAATATCACTC > NC_000913/16506‑16995 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |