| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,027,458 | 0 | A | . | 22.7% | 50.1 / 14.1 | 22 | intergenic (+230/+151) | fre/fadA | NAD(P)H‑flavin reductase/3‑ketoacyl‑CoA thiolase |
| Reads supporting (aligned to +/- strand): ref base A (10/7); new base . (5/0); total (15/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.35e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTC > NC_000913/4027260‑4027702 | aagaaCGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAT‑TGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCTCATCCGGCaagagg > 2:55047/1‑247 (MQ=11)aagaaCGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAT‑TGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTgg > 2:455701/1‑251 (MQ=11)aagaaCGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAT‑TGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTgg > 1:305515/1‑251 (MQ=11) agaaCGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAT‑TGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTggg > 1:176919/1‑250 (MQ=11) aaCGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAT‑TGGTGCCGGATCGGTAGACCGGATAAGGCGATCACGCCGCATCCGGCAAGTGGtt > 1:397517/1‑249 (MQ=11) agatctatgacgtatctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGAc > 1:58029/26‑251 (MQ=18) gatctatgacgtatctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGAcc > 1:475718/25‑251 (MQ=18) ctatgacgtatctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGAccc > 2:206459/22‑249 (MQ=18) tatgacgtatctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGaacc > 1:202749/21‑245 (MQ=18) gacgtttctggcaacccaggcccgATGCGACGCTGTCGCGACTTATCGTGCCTACAAATAGT‑CCG‑AACCGGAGGCCGGATAAGGCGGTTACGCCGCATCCGGCAATTGGTGCATGATGCCGGATGCGACGCTGGCGCGTCTTAGCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCAGCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGAcc < 2:164906/227‑1 (MQ=21) cgtatctggcagaccctgCCCGATGCGACGCTGTCGCGTCTTAGCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCTGGATAAGGCGTTTACGCCGCTTCCGGCAATTGGTGCATGATGCCTTATGCGTCGCTGGCGCGTCTTATCAGGCGTCCATTGGTGCCGGATCGGTAGACCGGGTAAGGCGGTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTGAAACACCGTCGCAATACCCTGACCCAGACCGAt < 1:131843/235‑1 (MQ=255) gtatctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATa < 1:310256/237‑1 (MQ=255) aaaccatgCCCGATGCGGCGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACAtcgtc > 1:106010/6‑251 (MQ=255) atgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGcc > 1:412508/1‑249 (MQ=255) gcccgATGCGCCGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGAc < 2:334302/247‑1 (MQ=255) aTGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAAc > 1:446678/1‑251 (MQ=255) gCGACGCGGTCGCGTCTTATCGTGCCTACAAATAGC‑CCG‑ACCCGTAGGCCGGATAAGTCGTTTACGCCGCATCCGGGAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTg < 2:202749/251‑1 (MQ=255) gCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGACAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGt < 1:344537/251‑1 (MQ=255) cgTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCg > 1:315348/1‑251 (MQ=255) tCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGtt > 1:370478/1‑251 (MQ=255) gTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAg < 1:359947/251‑1 (MQ=255) tACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTc > 1:287610/1‑250 (MQ=255) | AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTC > NC_000913/4027260‑4027702 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |