| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,888,549 | 0 | G | A | 68.2% | 16.2 / 12.2 | 22 | P400S (CCA→TCA) | fadD | long‑chain‑fatty‑acid‑‑CoA ligase |
| Reads supporting (aligned to +/- strand): ref base G (2/5); new base A (8/7); total (10/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.81e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.50e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGCTGCATGACGACATCTTCAATCTCGTTGGGATAGACGTTAAAACCGGAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCC > NC_000913/1888317‑1888753 | tgctgcATGACGACATCTTCAATCTCGTTGGGATAGACGTTAAAACCGGCAACCAGAATCATGACTCTTTCACGATCGCCAATGCGCAGGAATACTTCTTCATCCATTACCGCGATGTCGCAGCTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGAAACAAAACTCACCCGGTTGACATGGTGATACTTCATTatcatcatc > 1:289767/1‑251 (MQ=255) cATCTTCAATCTCGTTGGGATAGACGTTAAAACCGGAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGATACTTCAGTATCATCATCATCCACCAGTTTg < 1:86893/251‑1 (MQ=255) tCAATCTCGTTGGGATAGACGTTAAAACCGGAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTgg < 1:149306/247‑1 (MQ=255) cAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGGCCTTTGACACAAAGCTGACCCGGTTGACCTGGTGATACTTCATTATCATCATCATCCACCAGTGTTGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAAGATCATATGGGTTAACg < 1:56467/251‑1 (MQ=255) aTCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGATACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATgn < 1:96451/204‑2 (MQ=255) aTTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAgcgg < 1:248416/236‑1 (MQ=255) cgATGTCGCCGGTGTGTAACCAGCCATTTGTGATGATTTCATCGGTATCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATGATCATCATCATCCACCAGTTTGGCTTGCTGCTACGGCGCCGGGAATCCGATGCTACCACTATGGTAATCTAGATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAg < 1:242141/251‑1 (MQ=255) cGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAACATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGTCGCACACTCGGTAAGGCCATAGCCTTCcagn < 1:29481/251‑2 (MQ=255) tgtgtNNACNNNCCNNTTTGATGATTTCATCGGTAGCATNGGGACGCTGCCAGTAACCCAGCATCACCTNCNGNCCTTNNACANAAAGCTCCCCNGGATGACCTGGTGATACTTCATTATCATCATCATCCACAAGTTTGGCTNCNGTCGACGGCACCGGCAAACCGATGCTACCACNATGATAATCAATATCCTATGGGTTNACGcnngancnnnnnnnn > 2:237225/1‑207 (MQ=255) tgtgtAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTAGACACAAAGCTCACCCGGTTGACCTGGTGAGACTTCATTATCATCATCATCCACCAGTTTGTCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCt < 1:237225/208‑1 (MQ=255) cATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAGCCCAGCATCACCTGCGGTCGTTTGACACAAAGCTCGCTCTGTTGACCGGGGGATACTTCATTATCATCATCATCCACCAGTTTGGCtt < 1:148595/132‑1 (MQ=255) gatTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGATNNTTCATTnnc > 2:1984/1‑99 (MQ=255) gatTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGATACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTc < 1:1984/127‑1 (MQ=255) tCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGATACTTCATTATCATCATCATCCCCCAg > 1:147308/1‑107 (MQ=255) cATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGATACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGccacn < 1:81390/251‑2 (MQ=255) gCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAACGCTCACCCGGTTGACCTGGTGATACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATACTCAATATCCTATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTtgctgc > 1:46043/1‑251 (MQ=255) cTGCCAGTAACCGAGCATCACCTGCAGTCCTTTGACACAAAGCTCACCCGGTTGAACCGGTGGTACTTCCTTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGAGCAAACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTtgctgc > 1:107567/1‑250 (MQ=255) gTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCan < 1:31407/125‑2 (MQ=255) aCCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGATACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAAc > 1:11907/1‑110 (MQ=255) gTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGATACTTCATTATCATCATCATCCACAAGTTTGGCTTCCGTCgccg > 1:158616/1‑78 (MQ=255) cacaAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGcc > 1:174640/1‑232 (MQ=255) gCTCACCCGGTTGACCTGGTGATACTTCCTTATCATCATCATCCACAAGTTTGGCTTCCGTCGACGTCACCGGCAAACCGATGCTCCCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCAGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCCACGCTCTGcc > 1:206378/1‑197 (MQ=255) | TGCTGCATGACGACATCTTCAATCTCGTTGGGATAGACGTTAAAACCGGAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCC > NC_000913/1888317‑1888753 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |