| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,807,668 | 0 | A | . | 66.7% | 8.7 / 22.1 | 15 | coding (269/291 nt) | ydiZ | protein YdiZ |
| Reads supporting (aligned to +/- strand): ref base A (4/1); new base . (4/6); total (8/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.82e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAACCGGCAAAGATCCTGAAATCACCATCACGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTAC > NC_000913/1807592‑1807742 | aTTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAACCGGCAAAGATCCTGAAATCACCATCACGACCTGggaggatgga < 2:302128/81‑5 (MQ=255) agagAACTGGTTTCGGGGCTTACCCAAAGTGCAACCGGCAAAGATCCTGAAATCACCATCACGACCTGGGAggaatggaa < 1:302129/80‑1 (MQ=255) gagaACTGGTTTCGGGGCTTACCCAAAGTGCAACCGGCAAAGATCCTGAAATCACCATCACGACCTGGGAGG‑ATGGAATa < 2:302130/80‑1 (MQ=255) aCCCAAAGTGCAACCGGCAAAGATCCTGAAATCACCATCACGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGtt > 2:302162/1‑81 (MQ=255) aCCCAAAGTGCAACCGGCAAAGATCCTGAAATCACCATCACGACCTGGGAGG‑ATGGAATAGCAACAGAAAATAAATGGttt > 1:302157/1‑81 (MQ=255) cccAAAGTGCAACCGGCAAAGATCCTGAAATCACCATCACGACCTGGGAGG‑ATGGAATAGCAACAGAAAATAAATGGtttt < 2:302131/81‑1 (MQ=255) cccAAAGTGCAACCGGCAAAGATCCTGAAATCACCATCACGACCTGGGAGG‑ATGGAATAGCAACAGAAAATAAATGGtttt > 1:302166/1‑81 (MQ=255) ccAAAGTGCAACCGGCAAAGATCCTGAAATCACCATCACGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGtttt > 2:302161/1‑81 (MQ=255) ccAAAGTGCAACCGGCAAAGATCCTGAAATCACCATCACGACCTGGGAGG‑ATGGAATAGCAACAGAAAATAAATGGttttt < 1:302132/81‑1 (MQ=255) cAACCGGCAAAGATCCTGAAATCACCATCACGACCTGGGAGG‑ATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCaat < 2:302133/80‑1 (MQ=255) ccGGCAAAGATCCTGAAATCACCATCACGACCTGGGAGG‑ATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATc < 2:435305/81‑1 (MQ=255) aaGATCCTGAAATCACCATCACGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCt > 1:302159/1‑81 (MQ=255) aTCACCATCACGACCTGGGAGG‑ATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGtt < 1:302134/81‑1 (MQ=255) tCACCATCACGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGtt < 2:302135/81‑1 (MQ=255) cATCACGACCTGGGAGG‑ATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTc > 1:302156/1‑81 (MQ=255) aTCACGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTc > 2:302167/1‑81 (MQ=255) tGGGAGG‑ATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTa > 1:302164/1‑80 (MQ=255) ggaggaATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTac < 1:302136/80‑1 (MQ=255) | ATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAACCGGCAAAGATCCTGAAATCACCATCACGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTAC > NC_000913/1807592‑1807742 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |