| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,538,109 | 0 | A | C | 47.9% | 62.4 / 13.7 | 48 | G219G (GGT→GGG) | nanM | N‑acetylneuraminic acid mutarotase |
| Reads supporting (aligned to +/- strand): ref base A (18/7); new base C (0/23); total (18/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.04e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG > NC_000913/4537964‑4538253 | ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCttt > 2:402774/1‑151 (MQ=255) aaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTcac < 2:457068/151‑1 (MQ=255) cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAAc < 1:90666/151‑1 (MQ=255) cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAAc < 2:40437/151‑1 (MQ=255) ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAAcc < 1:585270/150‑1 (MQ=255) cAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAAccg > 2:613131/1‑150 (MQ=255) aTCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGc > 2:253427/1‑150 (MQ=255) tCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTAAAATAAAGTTCAAATACGGCATCCGTTCGAAACCCTGTTTTGTCTTCGCC‑ATTAATAAGCCAGTTTTTATCCCCTTTATTCACAGCCGCCGCACCAGCCg < 2:374712/151‑1 (MQ=255) cTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGt > 1:161201/1‑151 (MQ=255) gTGATGAGACGGGAGCAAGCTTATTCCTTTTTAAATTATTCCCGGTAAAATCAAGTTAAAATACGGCATCCTTTCGCAATCCTGTTTTGGCTCCGCC‑ATTAATAAGCCAGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTc < 2:95022/150‑1 (MQ=255) agaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCa > 2:527559/1‑151 (MQ=255) cGGGAGCAAGCTTATTCCATTTTAAATTATTCCCGGTAAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTGGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAggg < 2:120427/151‑1 (MQ=255) gggAGCAAGCTTTTTCCATTTTAAATTATTCCCGGTGAAATCAATTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCGGTTTTTATCCCCTTTTTTCCCAACCGCCGCACCAGCCGTTCCGTACCAggg < 2:205357/150‑1 (MQ=255) gaggaaaGCTTCTCCCTTTTTAATGTATTCCCGGTAAATTCAATTTAAAAACCGCCATCCTTTCCCTACCCTGGTTCGGCTTCGCC‑TTTAATATGCCGGTTTTTATCCCTTTTTTTTACAACCGCCGCACCAGCCGTTCCGTACCAGGGc < 1:570346/145‑1 (MQ=255) aGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg > 1:363116/1‑150 (MQ=255) tttttCCTTTTAAAATTATTCCCGTTAAAATAAATTAAAAATA‑GGCCACCCTTCCAAATTCTTGTTTTGCCTCCCCAATTAATAGGCCGGTTTTTACCCCTTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAg < 1:626961/148‑1 (MQ=255) ttttAAATTATTCCCGGTGAAATCAATTTAAAAGACGGCATCCGTTCGCAATCCTGTTTGGGCTCCGCC‑ATTAATAACCCGGTTTTTATCCCCTTTTTTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTc < 1:598430/151‑1 (MQ=255) ttttAAATTATTCCCGGTGAAATCAAGTAAAAATACGGCATCCGTTCGAAATCCTGTTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTc < 2:347309/151‑1 (MQ=255) attattCCCGGTAAAATCAAGTTAAAATACGGCACCCGTTCGAAATCCTGGTTTGGCTTCGCA‑ATTAAAAAGCCAGGTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc < 2:587568/151‑1 (MQ=255) attattACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCAGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg > 1:677130/1‑150 (MQ=255) ttattaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCGCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg > 2:189394/1‑149 (MQ=255) tattcCCGGTAAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGTTTGGGCTTCGCC‑ATTAATAAGCCAGTTTTTATCACCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg < 2:115649/151‑1 (MQ=255) ccGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg > 2:93420/1‑151 (MQ=255) ccGGTAAAATCAAGTTAAAATACGGCATCCGTCCGCAATCCTGTTTGGGCTCCGCA‑ATTAAAAAGCCAGTTTTTATCCCTTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg < 1:109869/151‑1 (MQ=255) taaaataaatTTAAAATACGCCTCCCGTTCCCAACCCTGTTTGGGTTCCCCC‑ATAAAAAACCCAGTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg < 1:527553/148‑1 (MQ=255) taaaaTCAATTTAAAATCCGGCATCCGTTCGCAATCCTGTTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg < 2:507337/148‑1 (MQ=255) taaaaTCAAGTTAAAATCCGGCATCCGTCCGCAACCCTGTTTTGGCTTCCCC‑ATTAATAAGCCAGGTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg < 2:361809/148‑1 (MQ=255) taaaaTCAAGTTAAAATACGGCACCCGTTCGCAATCCTGGTTGGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg < 2:590705/148‑1 (MQ=255) aaaacaaaGTTCAATTCCGGCACCCGTCCCCAATCCTGTTTTGGCTCCCCC‑ATTAATAAGCCAGGTTTTACCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGa < 2:129607/145‑1 (MQ=255) aaaaTCAAGTTAAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGa < 2:184410/150‑1 (MQ=255) aaaaTCAAGTTAAAATACGGCATCCGTTCGCAACCCTGTTTTGGCTCCGCC‑ATTAATAAGCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGa < 2:622924/150‑1 (MQ=255) aaaaTAAAGTAAAAATCCGGCTTCCGTCCCCATTCTGTGTTTGGCTCCGCC‑ATAAAAAACCCGGTTTTTACCCCTTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGa < 1:183574/150‑1 (MQ=255) aaTCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt > 1:455341/1‑150 (MQ=255) ttaaaaTACGGCTTCCTTTCCCAATCCGGTTTTGGTTCCGCC‑ATAAATAGCCCAGGTTTTACCCCTTTTATCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac < 1:253428/148‑1 (MQ=255) ttaaaaTACGGCATCCGTTCGAAATCCGGTTTGGGCTCCGCC‑ATTAATAAGCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAg < 2:292764/146‑1 (MQ=255) aaaTACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac < 2:189239/151‑1 (MQ=255) aaaTACGGCATCCGTTCGCAATCCGGGTTGGGCTTCGCC‑ATAAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac < 2:268530/151‑1 (MQ=255) cgcaaTCCTTCCGAAATCCTGTTTTGGCTTCCCC‑ATTAATAGGCCGGTTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa < 1:574640/147‑1 (MQ=255) gCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa > 2:320347/1‑149 (MQ=255) gCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa > 1:120625/1‑149 (MQ=255) aTCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata > 2:172653/1‑151 (MQ=255) ggTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt > 2:575918/1‑151 (MQ=255) ggTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt > 2:245819/1‑151 (MQ=255) aaGCCAGTTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAg < 2:645309/151‑1 (MQ=255) ggTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc < 2:292935/150‑1 (MQ=255) gTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt > 2:612039/1‑151 (MQ=255) ttttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAtt > 1:151002/1‑151 (MQ=255) tttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg > 2:450589/1‑151 (MQ=255) | CCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG > NC_000913/4537964‑4538253 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |