| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,177,748 | 0 | C | A | 27.3% | 39.1 / 12.0 | 34 | E55* (GAG→TAG) | nudF | ADP‑sugar pyrophosphatase |
| Reads supporting (aligned to +/- strand): ref base C (15/9); new base A (0/9); total (15/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.39e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCA‑GTAGCCAAGGGGTTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCC > NC_000913/3177505‑3177969 | cccccGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCA‑GTAGCCAAGGGGTTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAAtt > 2:190368/1‑251 (MQ=255)cccccGGGCTTGCCAGAAAACTTAACCCCGTTTTGGTCCGTTTGACTACCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCA‑GTAGCCAAGGGGTTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAAtt < 1:7355/251‑1 (MQ=255) cGGGCTTGCAAGAAAACTTAACCCCGTTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCA‑GTAGCCAAGGGGTTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTccc < 1:363811/251‑1 (MQ=255) aCTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTACAATCGCTTCGCGACGGGCAACATCTTCCACACTTTAACCCTCTTCAATCATCCCGGCAACCATCACCA‑GTAGCCAAGAGGTTT‑CGCTGGTGTCGTACGCGGCAATCCGTATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GGCGCACCCCATGat > 1:76794/1‑250 (MQ=255) cacCGGGTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCA‑GTAGCCAAGGGGTTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGGAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCAttt > 2:43485/1‑251 (MQ=255) tttcactaccacccccGCCTCTCCAATCGCTTCGCGACGGGCACCATCTTCCACACTTCCCCCCTCTTCAATCATCCCGGCACCCATCTCCA‑GTAGCCAAGGGGTTG‑CGGGGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGCCCGCGCTAAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGa < 2:380058/237‑1 (MQ=255) cacacccacccccGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCCCTCTTTCCCCCTCTTCAATCATCCCGGCAACCATCTCCA‑GTAGCCAGGGGTTTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TAAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGAt < 1:202790/241‑1 (MQ=255) ccttccacccccGCCTCTTCAATCGCTCCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCACCCCGGCAACCATCTCCA‑GTAGCCAAGGGGTTT‑CGCTGGGGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCATTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTCCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg < 1:281701/242‑1 (MQ=255) cAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCA‑GTAGCCAAGGGGTTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGCATAGACGCTGACGa > 1:52238/1‑250 (MQ=255) ccccGCCTCTTCAACCGCTTCGCGACGGGCAACATCTCCCCCACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCA‑GTAGCCAGGGGGTTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCCCGCACTGGG‑TCAAAGGGTAGCAA‑GACTGGGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATc < 1:121213/250‑1 (MQ=255) tcaagcacTTCGCGACGGTCAACATCTGCCTCGCTTCCACCCTCTACAATCACCCCGGCAACCCTCTCAA‑ATAGC‑AAGGGGTTTTCTGTGGTGACTTACGCGGCAATCAGAATCGTCGCAATCTGCACAAC‑TCCATCACGTTCTGGG‑CAAAAGGGTAGCAA‑GACTGCGGG‑GGGACCGCTCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGAt < 1:330704/243‑1 (MQ=255) cctccgcgACGGCCAACATCTTCCACACTTCCCCCCTCTTCAATCACCCCGGCAACCACCCCCA‑GTAGCCAAGGGGTTC‑CGCTGGTGTCGTACGCGGCAACCCGAATCTGCTCAACCACCACAAC‑TTCACCACGCACTGGG‑TAAAAGGGTAGCAA‑GCCTGCGGC‑GGGACCGCGCTAAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGa < 1:273924/247‑1 (MQ=255) gttcgtgaggggCCACATCTTCCACACTTTCACCCCCCTCAATCATCCCGGCAGCCATCTCCA‑GTGGCCAAGGGGTTG‑CGCTGGGGTCGACCGCGGCAATCCGATTCTGCACAATCAGCCCAAC‑TCCACCACGCACTGGG‑TCAAAGGGTAGCAT‑GACTGCGGC‑GTGACCGCGCTAAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaa < 2:179867/242‑1 (MQ=255) cTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCA‑GTAGCCAAGGGGTTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGCATAGACGATGACGAAATCTAGAAAGATCTAATGat > 2:27120/1‑250 (MQ=255) cATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCA‑GTAGCCAAGGGGTTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATa > 2:87409/1‑250 (MQ=255) tctccccacattCAATCCTCCCGCAAACCATC‑CCATGTAGCCAGGCGTTTT‑CGGTGGTGGCGTACGCGACAACCCGAAGCGGCTCAAGCAG‑ACAACCTACCTCTCGCACGGGG‑TAAAAGGGTACCAAAGA‑TGCGGC‑GGGACCGCGCT‑AAAAATTTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCTCGATAAAGTGTTTCTCg < 2:97765/240‑1 (MQ=255) tctTCAATCATCCCGGCAACCATCTCCA‑GTAGCCAAGGGGTTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaataat > 2:125067/1‑251 (MQ=255) tcaacccccccggacaccacccccA‑GTAGCAAGGGGGTTT‑CGCTGGTGTCGTACGCGCCAATCCAAATCTGCCAAACCCGCCAAAC‑TTCATCCCCCACTGGT‑TCAAAGGGAAGAAA‑GCGGGGGGC‑GTGACCGCGCTAAAAAATTTCCC‑GCCGCACCTCATGACTCAGTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAGAAAAAGCCGCGATAAAGTGTTTCTCGTGCaataat < 2:291746/226‑1 (MQ=255) cGGCAACCATCTCCA‑GTAGCTATGGGGTTT‑CGTTGGTGGCGTACGCGGCAACCCGAATCTGCCAAACCATCACATC‑TTCCTCACGC‑CTGGGGTCAAAGGGTAGCAA‑GCCTGCGGC‑GTGACCGCTCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGttt < 1:212456/251‑1 (MQ=255) aCCATCTCCA‑GTAGCCAAGGGGTTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGcc > 2:153076/1‑251 (MQ=255) ctcCA‑GTAGCCAAGGGGTTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGAACGCGCTCAAAACTTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATATAAGCTGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCaactg > 2:72753/1‑248 (MQ=255) ccA‑GTAGCCAAGGGGTTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGt > 2:184755/1‑250 (MQ=255) gtagccaggtgttttggggggTGGCGTACGAGGCAATCCTAATCTGCACAACCAGCACAAC‑TTCCTCCGGCACTGGG‑TCAAAGGGTCGAAA‑TCCT‑CGGCCGGGACCCCCCTAAAAATTTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACg < 2:76794/232‑1 (MQ=255) ggTGTCGTACGCGTCAATCCGATTCTCCTCAATCACCACTTC‑TTCATCACCCGCGGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGTCCGGGCTAAAAATTTTCCC‑GCCGCACCTCATGACTCATTTTCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTGCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTaa < 1:292584/251‑1 (MQ=255) gtCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGc > 1:145114/1‑250 (MQ=255) ggCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATg > 1:29350/1‑251 (MQ=255) gAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGATGTCTGGCTTAAGCATTGTTAATGTCCTGGc > 2:47948/1‑250 (MQ=255) cccgccttggttaaaggggTAGAAA‑GACTGCGGC‑GTGCCCGCGCAAAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCGCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCa < 2:204055/235‑1 (MQ=255) tGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGt > 1:55118/1‑251 (MQ=255) ggtacAAAGGGTAGCAA‑GACTGGGGC‑GTGCCCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGt < 1:331266/246‑1 (MQ=255) ggtaaaggggTAGCAG‑GACGGCGGC‑GTGCCCGCGCTAAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGttt < 1:306419/244‑1 (MQ=255) ggtaaaggggTAGCAA‑GCCTGCGGC‑GTGCCCGCGCAAAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGttt < 2:29350/244‑1 (MQ=255) gggTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGcc < 1:43806/251‑1 (MQ=255) aa‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGCTCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCg > 1:17574/1‑125 (MQ=255) | CCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCA‑GTAGCCAAGGGGTTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAAC‑TTCATCACGCACTGGG‑TCAAAGGGTAGCAA‑GACTGCGGC‑GTGACCGCGCTCAAAAATTTCCC‑GCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCC > NC_000913/3177505‑3177969 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |