Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,608,719 | 0 | G | A | 26.7% | 76.9 / 12.1 | 45 | intergenic (+70/+32) | zntA/tusA | Zn(2(+))/Cd(2(+))/Pb(2(+)) exporting P‑type ATPase/sulfur transfer protein TusA |
Reads supporting (aligned to +/- strand): ref base G (24/9); new base A (0/12); total (24/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.02e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAA‑CCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTG > NC_000913/3608587‑3608856 | ggaGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAA‑CCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAggc > 1:657709/3‑140 (MQ=255) gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAA‑CCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAg > 1:238219/1‑140 (MQ=255) gacgggCTGGTGACAGCGAATGCGTAAAGATTGTGGCGCGGGAGAAAAGGCAAC‑CCGATCGAAACATGGAGCGCGATGGGCCCCCTCCCCCC‑TCGGGGGGGGGGGTTGGGG‑TGGGGGAAAAAGGCGGCATCGAAGCCAAt < 2:1031039/135‑1 (MQ=255) aacaggCTGGTGAAAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGAAAC‑CCGACCGCAACATTGAGCGCGATCGGCCCCCCCGCCCC‑TCGGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAAt < 2:501508/135‑1 (MQ=38) aCGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAA‑CCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAAt < 2:1174650/140‑1 (MQ=255) gTGCTGGTGACAGCGATTGGTTTAAGATTGTTGCGCGGGAGATAAGGAAAC‑CCGATCGCAACTTTGGGCGCGATCGGCCCCCTCCCCCC‑TTTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAg < 2:1314777/140‑1 (MQ=37) gTGCTGGTGACAGCGAATGCGTAAAGATTGTTGCGCAGGAGAAAAGGAAAC‑CCGATCGCAACATTGAGCGCGATCGGCCCCCCCCCCCC‑TGGGGGGAGGGGGTTGGGG‑TGGGGGGAAAAGGCGGCATCGAAGCCAATCAg < 2:106448/140‑1 (MQ=255) tGCGGGGGAAAGGGAATGCGTTAAGATTGTAGCGGAGGAGATAGGGAAAA‑CCGAGCGCAAATTGGGGCGC‑AGCGGCCCCCCCCCCCCGTGGGGGGAGGGGGTTTGGG‑GGAGGGAAAAAGGCGGCATCGAAGCCAATCAGc < 2:1255155/140‑1 (MQ=255) tGGTGACAGCGAATGCGTTAAGATTGTGGCGCAGGAGAAAAGGAAAA‑CCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGGGGGGAAAAGGCGGCATCGAAGCCAATCAGcccc < 2:1217043/140‑1 (MQ=39) aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAA‑CCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc > 1:1130401/1‑140 (MQ=255) cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAA‑CCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATcctgtc > 1:1031039/1‑135 (MQ=255) cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAA‑CCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATcctgtc > 1:501508/1‑135 (MQ=255) ttAAGTTTGTTGCGAAGGAGAAAAGGCAAC‑CCGATCGAAACTTGGAGCGCGATCGGCCCCCCCCCCCC‑TGGGGGGGGGGGTTTGGGG‑TGGGGGAAAAAGGCGGTATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg < 2:1353630/140‑1 (MQ=255) gttgCGCAGGAGATAAGGCAAA‑CCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt > 1:658235/1‑140 (MQ=255) gttgCGCAGGAGATAAGGCAAA‑CCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt > 1:1389437/1‑140 (MQ=255) ttgCGCAGGAGATAAGGCAAA‑CCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGTATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa > 1:1353630/1‑140 (MQ=255) ttgCGCAGGAGATAAGGCAAA‑CCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa > 1:672123/1‑140 (MQ=255) gggggAAATAAGGAAAC‑CCGATCGAAACATTGGGGGCGATCGGCCCCCCCCCCCC‑TGGGGGGAGAGGGTTGGGG‑TGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa < 2:639823/137‑1 (MQ=255) gcAGGAGAAAAGGCAAA‑CGGATCGCAACATGGAGCGCGATCGGCCCCCCCCCCCC‑TCGGGGGAGGGGGTTGGGG‑GGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa < 2:54171/140‑1 (MQ=255) cAGGAGATAAGGCAAA‑CCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGAt > 1:54171/1‑140 (MQ=255) aGGCAAC‑CCGATCGCAACATTGAGCGCGATCGGCCCCCTCCCCCC‑TTGGGGGAGGGGGTTAGGG‑TGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTc < 2:583441/140‑1 (MQ=255) agCAAACCCGATCGAAACATTGAGCGCGATCGGCCCCCTCCCCCC‑TCGGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTc < 2:518944/139‑1 (MQ=38) gaaacccGATCGCAAATTTGGGCGGGATGGGCCCCCCCCCCCT‑TTGGGGGAGGGGTTTGGGG‑GGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCa < 2:81183/135‑1 (MQ=255) aaa‑CCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATc > 1:48838/1‑140 (MQ=255) a‑CCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCg > 1:518944/1‑140 (MQ=255) ccGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt > 1:1245485/1‑140 (MQ=255) cgcgATCGGTCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt > 1:1352157/1‑140 (MQ=255) ttCGGCCCCCTCCCCCC‑TGGGGGGGGGGGGTTGGGG‑GGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg < 2:1245485/139‑1 (MQ=255) tgggtcccccTCG‑CCC‑TGGGGGGGGGGGGTTAGGGTTGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg < 2:1261656/134‑1 (MQ=255) aTCGGCCCCCCCGCCCC‑TCGGGGGAGGGGGTTAGGG‑TGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg < 2:1174505/140‑1 (MQ=255) aTCGGCCCCCCCCCCCC‑TCGGGGGAGGGGGTTAGGG‑TGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg < 2:1389437/140‑1 (MQ=255) cccccccGCCCC‑TCGGGGGGGGGGGTTGGGG‑TGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCa < 1:770259/139‑1 (MQ=255) ccccccccc‑TTGGGGGAGAGGGTTAGGG‑GGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa < 1:926820/135‑1 (MQ=255) gCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt > 1:1000248/1‑140 (MQ=255) tctGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa > 1:130093/1‑140 (MQ=255) gtaggggggTTAGGG‑GGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCTCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAccc < 2:227007/134‑1 (MQ=255) gggAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAccc > 2:542043/1‑140 (MQ=255) agagGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAgg > 1:444008/1‑140 (MQ=255) gagGGTTAGGG‑TGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGa < 2:658235/140‑1 (MQ=255) gagGGTTAGGG‑TGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGa < 2:48838/140‑1 (MQ=255) agGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa > 1:132084/1‑140 (MQ=255) ggTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa < 2:1237122/127‑1 (MQ=255) ggTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa > 1:1237122/1‑127 (MQ=255) aGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAg > 2:513820/1‑140 (MQ=255) aGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAg > 1:1101301/1‑140 (MQ=255) aGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCATCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAg > 1:42267/1‑140 (MQ=255) aGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTACATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAg > 2:1633965/1‑140 (MQ=255) gggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTg > 1:1213912/1‑140 (MQ=255) | GCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAA‑CCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCC‑TCTGGGGAGAGGGTTAGGG‑TGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTG > NC_000913/3608587‑3608856 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |