| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,177,748 | 0 | C | A | 46.5% | 34.0 / 29.0 | 43 | E55* (GAG→TAG) | nudF | ADP‑ribose pyrophosphatase |
| Reads supporting (aligned to +/- strand): ref base C (15/8); new base A (0/20); total (15/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.34e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCA > NC_000913/3177603‑3177891 | tCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCaaaa > 2:436755/1‑150 (MQ=255) aCCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAAtt > 1:855508/1‑151 (MQ=255) ctctTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTcc > 2:152591/1‑151 (MQ=255) tctTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTccc > 1:549855/1‑151 (MQ=255) tCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGATCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccg < 1:550699/149‑1 (MQ=255) cAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgccg > 2:769635/1‑151 (MQ=255) cAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgccg > 1:589798/1‑151 (MQ=255) aaTCACCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgccg < 2:496030/150‑1 (MQ=255) cccGGCACCCATCCCCAGTAGCCAGGGGG‑TTT‑CGCTGGTGTCGTACGCGGCAACCCGAATCTCCTCAATCACCACAACTCCACCACGCCCTGGGTCAAGGGGTAGCAAGACTGGGGCGTGCCCGCCCTAAAAAATTTCCCGCCGCACCTc < 1:665494/150‑1 (MQ=255) cccGGCACCCATCCCCAGTAGCAAGGGGT‑TTT‑CGCGGGTGTCGTACGCGGAAACCCGAATCTGCTCAATCACCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGGGTGACCGCGCTAAAAAATTTCCCGCCGCACCt < 1:91410/149‑1 (MQ=255) cGGCAACCATCTCCAGTAGCCAAGGGG‑TTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATg > 2:520804/1‑151 (MQ=255) ggCAACCATCTCCAGTAGCCAAGGGG‑TTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATg > 2:98560/1‑150 (MQ=255) acccATCCCCAGTAGCAAGGGTT‑TTT‑CGCTGGTGTCGTACGCGGCAATCCAAATCTCCTCAATCAGCACACCTCCCCCCCGCACGGGGTCAAAGGGTAGCAAGACTGCGGCGAGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGa < 2:423875/146‑1 (MQ=255) aCCATCTCCAGTAGCCAAGGGG‑TTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTc > 2:316977/1‑150 (MQ=255) ccccgtaagcAAGGGGTTTT‑CGCTGGTGTCTTCCCGGGCAACCCAATTCTCCCAAACCAGCACACCTTCACCCCCCACTGGGTCAAAGGGTACCAAGCCTGCGGCGTGACCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCAtt < 1:498254/139‑1 (MQ=255) aGTAGCCAGGGGG‑TTT‑CGCGGGGGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCAGCACGCACGGGGGCAAAGGTTACCAAGACTGCGGCGTGACCGCGCTAAAAAATTGCCCGCCGCACCTCATGACTCATTTGCCCGt < 2:430908/151‑1 (MQ=255) gCCAAGGGG‑TTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGa > 1:619975/1‑150 (MQ=255) agcgaGGGG‑TTTGCGCGGTTGTCGTACGCGGAAAACCGTATCTTCGCAGACAACACACCTTCTTCCCGCTCTGGTTCAAAGGGTATCAAGACGGCGGCGTACCCGCGCTCAAAATTTCCCCGCCGCACCTCATGACTCATTTGCCCGtt < 1:453763/145‑1 (MQ=255) aaGGGG‑TTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTCCACCACCCACGGGGTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATa < 2:262085/150‑1 (MQ=255) aaGGGG‑TTT‑CGCTGGTCTGGTACGCGGCAACCCAAATCCGCCCAATCAGCACAACTCCAGCACGCACGGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGGCCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGa < 1:671904/147‑1 (MQ=255) ggtg‑TTT‑CGCGGGTGGCGTCCGCGGAAACCCGAATCTGCCCAATCAGCAAAACTTCACCACCCACGGGGTCAAAGGGTAGCAAATCTGCGGTGTGCCCGCAAAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGa < 1:111914/147‑1 (MQ=255) tccGCTGGTGTCGTACGCGGCAACCCGAATCTGCTCATCCAGCACAACTGCACCCCGGCCTGGGTCAAAGGTTAGCAAGCCTGTGGCGTGACCGAGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGAt < 1:604608/147‑1 (MQ=255) tGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaa > 1:741550/1‑150 (MQ=255) cgcgGAAACCCGAATCTGCTCAATCAACACACCTTCATCACGCACTGGGTCAAAGGGTAGCAAGCCTGCGGCGTCCCCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAg < 2:359702/150‑1 (MQ=255) aaTCCGAATCTGCTCAATCACCACAATTTCATCCCGCACTGGGTAAAGGGTTAGAAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTa < 1:740758/149‑1 (MQ=255) cTCAATCTGCACAGCTCCATCACGCATGTTGCCAAGGGCGAGCAATATGGCGGGGTGCCCGCGATCAAAAGTTTCCCGCCGCACCTCATGACTCATTTGCCAGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAg < 1:66842/148‑1 (MQ=255) cTCAACCAGCACAACTTCACCACCAACTGGTCCAAAGGTTAGCAGGCCTGCGGCGTGCCCGCGGTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaaaa < 2:62380/147‑1 (MQ=255) caccacaACTCCACCCCCCACTGGTTAAAAGGTTAGCAAGACTGCGGCGTGACCGCGCCAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATa < 1:330789/147‑1 (MQ=255) aGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGAGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATa > 2:362200/1‑149 (MQ=255) acctccatcaCGCACGGGATAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGAt < 1:520766/138‑1 (MQ=255) ttCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtct > 1:287227/1‑151 (MQ=255) accgcGCACTGGTAAAAAGGTTAGAAAGCCTGCGTGGTGACCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtctc < 1:692594/145‑1 (MQ=255) caCGCAGGGGTTAAAGGGGTAGCAAGACGGCGGCGTGCCCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGc < 1:606236/151‑1 (MQ=255) caCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTg < 2:653754/150‑1 (MQ=255) caCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGc > 2:146315/1‑151 (MQ=255) aaaGGGTACCAAAACGGGGGCGTGACCGCGATAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCt < 1:34305/148‑1 (MQ=255) gggTAGCAAGACTGCGGCGTGCCCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCg < 1:63300/151‑1 (MQ=255) gTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGtt < 2:525556/151‑1 (MQ=255) aaGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGc < 2:614075/151‑1 (MQ=255) cggcggcgTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCa < 1:808250/147‑1 (MQ=255) gACCGCGCCAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAAc < 1:184872/148‑1 (MQ=255) ccGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGc > 2:690790/1‑150 (MQ=255) ccGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCa < 2:287252/151‑1 (MQ=255) | TCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGG‑TTT‑CGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCA > NC_000913/3177603‑3177891 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |