| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,093,139 | 0 | C | A | 28.9% | 65.6 / 11.0 | 38 | V100V (GTG→GTT) | rhaM | L‑rhamnose mutarotase |
| Reads supporting (aligned to +/- strand): ref base C (10/17); new base A (11/0); total (21/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.03e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCAGGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGC > NC_000913/4092932‑4093361 | cAGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCAGGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTTAAAAACTTCTTTCAACTCg > 1:63827/1‑223 (MQ=255) aaTGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCAGGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTGAAAAACTTCTTGCAGCTCGCCACT‑CACACGGCTGTTTTCCCCGTTAg > 2:447263/1‑228 (MQ=255) aTTCATTATTTTTGTGATTTACTCACAATGTATCAGGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTTCAGCTCGCCACT‑CACCCGGCTGTGATCACGGTTAGCGGGGATAAAATCGGTCATAGA‑TTTCCACCAC‑AGTTGGcc > 1:231651/1‑227 (MQ=255) aTTCATTATTTTTGTGATTTACTCACAAGGTATAAGGGAAATATAACTTACCGGGAGATGTAATGTGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCa < 2:635160/228‑1 (MQ=255) aTCAGGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTTAAAAACCTCTTTCAACCTGCTACT‑CACCCGTCTGTTATTCCGGTTATCGGGGCTAACATCGGGCTCATA‑TTTCTCCCCA‑CAATGGCAAAAATCA‑TGTCTGTCAAACGCCTTcccacgt > 2:651894/1‑223 (MQ=255) gggAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑Gttctt < 1:447996/228‑1 (MQ=255) aaaTATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGCCCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAAAATCGGGCATATA‑TTTTCACCAACCGTTGGCACAAACCA‑‑TGC‑GGCAACCGCATTCCaagcgtcctcag > 2:277399/1‑216 (MQ=255) atatAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGc > 2:405843/1‑114 (MQ=255) atatAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAAAATTCTTGCAGGTCGCTAAT‑CACCGGGCTGTTATCCGGGTTAGCGGGGATAACATAGGTCACATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGGAAACGCATTCCAAA‑GTTTTTCTGa > 1:19642/1‑227 (MQ=255) atatAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGCTCGCTACT‑CACCGGGCt > 2:124830/1‑134 (MQ=255) tataACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCa < 1:406505/114‑1 (MQ=255) tataACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTg < 1:125013/134‑1 (MQ=255) cTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGCTCGCTACTCCCCCGGCCTGTATTCCGGG‑TAGCGGGCATCACATCGGTCATATA‑TTTTCACCCA‑CACTTGGAAACATTA‑ATGCTGGCAAACGCATTTCAGC‑CTTCTTCAGAt > 2:679930/1‑222 (MQ=255) ttACCGGGATATGTAATAGGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGAtt < 1:681043/222‑1 (MQ=255) tACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATGTTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAAACTTCTTGCAGCTCGCTACT‑CAc > 1:145501/1‑120 (MQ=255) aCCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATGTTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CAcc < 2:145290/120‑1 (MQ=255) aCCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAAACTTCTTTCAGCTCCCTACT‑CACCCGGCTGTTTTCCCGGGTAGCGGGCATAACAACGGTCATATA‑TTTCCACCAC‑AGTTGGGAAACATCA‑ATGCTGGCAACCGCCATCCAGA‑GTGCTTCAGAATCAATatct > 2:648551/1‑228 (MQ=255) ggAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTGCAGCTCGCTACT‑CACCCGGCTGTTATCCGGGGTAGCgg > 2:634728/1‑138 (MQ=255) gagaTGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCggg < 1:635739/138‑1 (MQ=255) atGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTCAAACAATTCTTTTAGCTCGCCACT‑CACCGGGGAGTTAt > 2:13266/1‑115 (MQ=255) tGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATc < 1:13286/115‑1 (MQ=255) gTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAAACTTTTTTCAGCCCGCTACT‑CACCGGGGTGTTTTCCCGGGTAGCGGGCATAACATCGGGCCTATATTTTCCACCAC‑GCTTGGCAAAC‑TCAGGGGCTGG‑AAACGCATTCCAACGGTATTTCAGTTTC‑ACATCTATCATGGCAAAcagcag > 2:639581/1‑228 (MQ=255) tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAg > 1:389260/1‑121 (MQ=255) ttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTTTTGCAGCTCGCCACT‑CTACCGGCTGTTATTCCGGTTAACGGGCCTAACATCGGTCATATA‑TTTTCACCAA‑CCTTTGCAAACATCA‑GTGCCGGCAACCGACTTCCTGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAAcag > 2:84964/1‑223 (MQ=255) ttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTAc > 2:363512/1‑96 (MQ=255) ttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGc < 2:388619/121‑1 (MQ=255) ttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCCGGCTGGTATTCGGGGTAGCGGGCATAACATCGGGCATATA‑TTTTCACCAA‑CGTTGGCAAACATCA‑ATGCTTGCAACCGCCTTCCATC‑GTTCTTCATATTCAATCTCTACCATGGccacaagc > 2:327137/1‑217 (MQ=255) tAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACt < 1:364116/96‑1 (MQ=255) tctcTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGt < 2:63735/228‑1 (MQ=255) ctGGCACATTGGGCAATTCCGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTga < 2:231295/227‑1 (MQ=255) ggCAATTACGGCAGGTAAAACACTTCTTGCAGCCCCCTACT‑CACCCGGCTGTTTTCCGGGGTAGCGGGCATAAAAACGGTCATATA‑TTTCCCCCCA‑CGTTGGCAAACATCA‑ATGCTGGCAAACGCATTCCACC‑GTTCTTTCGATTTAAACTCTACCAAGGCAAACTACAGGATACGCGCTTCGTCGAGATAGATGGc > 2:385959/1‑197 (MQ=255) gCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCGCTTTGTCGAGATAGATGGCg < 1:386599/197‑1 (MQ=255) aaTTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCGCtt < 1:398274/178‑1 (MQ=255) ttACGGCAGGGAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCCCTTTGTCGAGATAGATGGCGTAGTTATTCGGAACGTg > 1:403947/1‑210 (MQ=255) tACGGCAGGTAAACCACTTCTTGCAGCTCGCTACT‑CACCGTGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTga < 2:403289/210‑1 (MQ=255) tACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTAcgcg > 1:4930/1‑172 (MQ=255) aCGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACgcgc < 2:4924/172‑1 (MQ=255) tAAAACACTTCTTGCAGCTCGCTCCT‑CACCGGGTTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGc < 1:239555/228‑1 (MQ=255) | CAGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCAGGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGC > NC_000913/4092932‑4093361 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |