| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,596,228 | 0 | G | C | 23.1% | 52.2 / 14.3 | 26 | A59G (GCC→GGC) ‡ | livH | branched‑chain amino acid ABC transporter permease |
| Reads supporting (aligned to +/- strand): ref base G (11/9); new base C (3/3); total (14/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AATCAGGCGCTTAGAGTTACGCACCGGGCGGTAAGCCACCCGTTCGATACTCCAGCCGTAGGCGCTGGCAATGACGATTGCGCCGACGAATCCCGCAGCTACCAGCAGCCAGCCGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAA > NC_000913/3596086‑3596367 | aaTCAGGCGCTTAGAGTTACGCACCGGGCGGTAAGCCACCCGTTCGATACTCCAGCCGTAGGCGCTGGCAATGACGATTGCGCCGACGAATCCCGCAGCTACCAGCAGCCAGCCGGTATCAATGCCCATCATCATCAGCGCGGCgatga > 6:98856/1‑149 (MQ=255) gagTTACGCACCGGGCGGTAAGCCACCCGTTCGATACTCCAGCCGTAGGCGCTGGCAATGACGATTGCGCCGACGAATCCCGCAGCTACCAGCAGCCAGCCGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATgaga < 5:98856/148‑1 (MQ=255) cagTTCGATACTCCAGCCGTAGGCGCTGGCAATGACGATTGCGCCGACGAATCCCGCAGCTACCAGCAGCCAGCCTGTATCAATGCCCATCATCATCCGCGCGGCGATGATCATAAATCAGACGTAGCt > 2:109458/3‑129 (MQ=255) aGCCGTAGGCGCTGGCAATGACGATTGCGCCGACGAATCCCGCAGCTACCAGCAGCCAGCCGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTg < 2:28879/149‑1 (MQ=255) gaATCCCGCAGCTACCAGCAGCCAGCCGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATgg < 6:4030/143‑1 (MQ=255) gaATCCCGCAGCTACCAGCAGCCAGCCGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATgg > 5:4030/1‑143 (MQ=255) gCTACCAGCAGCCAGCCGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATGGTGTAGCCGATGGCTAt < 2:84358/149‑1 (MQ=255) gcagcCAGCCGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATc > 6:85644/1‑101 (MQ=255) gcagcCAGCCGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATGGTGTAGCCGATGGCTATCAGcgcg > 2:152768/1‑149 (MQ=255) cagccGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAgg > 1:82051/1‑148 (MQ=255) cagccGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAg < 2:74876/147‑1 (MQ=255) cagccGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAg > 1:74876/1‑147 (MQ=255) agccGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGaa > 2:28342/1‑88 (MQ=255) agccGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGaa < 1:28342/88‑1 (MQ=255) cGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCata > 3:126442/1‑65 (MQ=255) cGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCata < 4:126442/65‑1 (MQ=255) atttatgatcatcgCCGCGCTGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGa > 8:99998/8‑149 (MQ=255) atttatgatcatcgCCGCGCTGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGa < 7:99998/142‑1 (MQ=255) ttatgatcatcgCCGCGCTGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATGgtg < 8:32757/102‑1 (MQ=255) ttatgatcatcgCCGCGCTGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATGgtg > 7:32757/6‑107 (MQ=255) ttatgatcatcgCCGCGCTGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATGgtg < 6:55328/102‑1 (MQ=255) ttatgatcatcgCCGCGCTGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATGgtg > 5:55328/6‑107 (MQ=255) catcatcatCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATg < 2:145203/103‑1 (MQ=255) catcatcatCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATg > 1:145203/1‑103 (MQ=255) atcatcaGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGt < 1:152768/148‑1 (MQ=255) tcaGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTaaa > 3:889/1‑148 (MQ=255) | AATCAGGCGCTTAGAGTTACGCACCGGGCGGTAAGCCACCCGTTCGATACTCCAGCCGTAGGCGCTGGCAATGACGATTGCGCCGACGAATCCCGCAGCTACCAGCAGCCAGCCGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATGGTGTAGCCGATGGCTATCAGCGCGTAGGTACTGCCCAGCGTGACGCCGTTAAA > NC_000913/3596086‑3596367 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |