Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,926,366 | 0 | A | C | 28.6% | 77.9 / 11.4 | 42 | G30G (GGT→GGG) | mioC | flavoprotein MioC |
Reads supporting (aligned to +/- strand): ref base A (24/6); new base C (0/12); total (24/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.68e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAACG > NC_000913/3926227‑3926500 | tCTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGCCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTaaaa < 2:263419/140‑1 (MQ=255) gCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCCCCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCctct < 2:400535/140‑1 (MQ=255) aTGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCCCCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGCCCGTGCAGCGTTTC‑GGTGGTAAAC‑CCCGCCTCTTcc < 2:278827/140‑1 (MQ=255) aTGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTcc > 2:212677/1‑140 (MQ=255) aTGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTcc > 2:332661/1‑140 (MQ=255) gCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAg > 2:163115/1‑140 (MQ=255) taaTAGAAAGGAGAAAGGTTGTCCGAAATATCCCGGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACGGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCtt < 2:403813/138‑1 (MQ=255) tCAAAGAAGAGAGAAAGGTTGTCCGAAATACCCCCGCCACCGTGGGTGGAGCTGATACCCACCCAGATCCCTGAGCCAGGTAAATCTTCAAACAGCGGACGGTGCGGCGTTCG‑GGTGGTAAAC‑CCCGCCTCTTCCAGCtt < 2:550424/140‑1 (MQ=255) aaaggagaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGc > 1:734573/1‑140 (MQ=255) aaggagaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGcc < 1:319537/140‑1 (MQ=255) ggagaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAg < 1:163115/140‑1 (MQ=255) agaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAggg > 2:396482/1‑139 (MQ=255) agaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGgt > 2:155269/1‑140 (MQ=255) agaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGgt > 2:690949/1‑140 (MQ=255) gaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGgtg < 2:507473/140‑1 (MQ=255) aaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGgtgt > 1:364023/1‑140 (MQ=255) aaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGtt < 1:690949/140‑1 (MQ=255) gttggTCCGGAATATGTCACCCCCCGTGGGGGGGGGGACTACCCAGCCAGTTCACTGACGCAGGAACATCTTCAAACAGCGGGCAGCCCCGCGTTGC‑GGTGGAAAGC‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCt < 2:568041/136‑1 (MQ=255) gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTa > 1:177958/1‑140 (MQ=255) tCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGGAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCGTCCAGCTTTTCAGCCAGGTGTTCTGCTAc > 2:567276/1‑140 (MQ=255) cGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACat > 2:682673/1‑140 (MQ=255) gaaaTATCCCCGGCCCCGGGGGTGGAGCTGATACCCAGCCAGACCCCGGGGGAAGGAAAATCTTCAAACAGCGGCCGGTGCGGCTTTTG‑GGGGGAAAAC‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACata < 2:494308/138‑1 (MQ=255) gggCTAATACCCACCCAGCCCCCGGAGGCAGAAAATTCTTCAAACGGGGGCCGGTGCGGTTTTTG‑GGTGAAAAAC‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGcc < 2:734573/138‑1 (MQ=255) gCTAATACCCAGCCAGATCCCTGGGGCAGGTAAATCTTCAAACGGCGGACGGTGCAGCGTTTG‑GGTGGAAAAC‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGc < 2:290337/140‑1 (MQ=255) gccagATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTg > 1:755738/1‑140 (MQ=255) agATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGata > 2:656190/1‑140 (MQ=255) aTCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGatact > 1:777281/1‑138 (MQ=255) ccTGAGGCAGGAAAATCTTCTAACAGCGGCCCGTGCAGCTTTTC‑GGTGGAAAAC‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGc < 1:567276/140‑1 (MQ=255) gAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCAt > 1:31928/1‑140 (MQ=255) gAGGCAGGTAAATCCTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCAt > 1:280448/1‑140 (MQ=255) gAGGCAGGAAAATCTTCAAACAGCGGACCGTGCAGCTTTTG‑GGTGGAAAAC‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCAt < 1:159613/140‑1 (MQ=255) gaaaaTCTTCTAACAGGGGCCCGTGCGGCTTTTG‑GGTGGAAAAC‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc < 2:280448/138‑1 (MQ=255) gTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc > 1:50209/1‑140 (MQ=255) aatttttcaAACAGCGGCCGGTGCAGCTTTTG‑GGGGGAAAAC‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCa < 2:364023/135‑1 (MQ=255) aaaTTTTCAAACAGCGGCCCGGGCAG‑GTTTCGGGTGGTAAAC‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCa < 1:458026/140‑1 (MQ=255) aaaTCTTCTAACAGGGGCCCGTGCAGCGTTTG‑GGTGGTAAAC‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCa < 1:261508/140‑1 (MQ=255) tcttcaAACAGGGGCCGGTCCGGCTTTTG‑GGGGGAAAAC‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCt < 2:50209/140‑1 (MQ=255) tAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTttatt > 2:36572/1‑140 (MQ=255) cgggcagTGGAGCGTTT‑‑GGGGGTAAAACCCCGCCTCTTCCAGTTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGt < 2:755738/134‑1 (MQ=255) ggcccGTGCAGCTTTTC‑GGTGGAAAAC‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTg < 1:396482/137‑1 (MQ=255) ggACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTg > 1:115948/1‑140 (MQ=255) ggACCGTGCAGCGTTTC‑GCTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGGGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTg > 1:458473/1‑140 (MQ=255) cccGTGCAGCGTTTC‑GGTGGTAAAC‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGc < 1:598262/139‑1 (MQ=255) tGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTAt > 1:104807/1‑140 (MQ=255) tGCAGCGTTTC‑GGCGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTAt > 2:285469/1‑140 (MQ=255) tttC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGc > 1:305971/1‑140 (MQ=255) gggaaaac‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTCCATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGaa < 2:618927/136‑1 (MQ=255) gtAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAAcg > 2:375577/1‑140 (MQ=255) | TCTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTC‑GGTGGTAAAA‑CCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAACG > NC_000913/3926227‑3926500 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |