Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,022,395 C→T 27.1% G646D (GGC→GAC)  fadB ← fused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase/3‑hydroxyacyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,022,3950CT27.1% 21.4 / 5.6 19G646D (GGC→GAC) fadBfused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase/3‑hydroxyacyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (6/7);  new base T (3/2);  total (9/10)
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.93e-01

CCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGATAATGCCTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTTACCT  >  CP009273/4022251‑4022541
                                                                                                                                                 |                                                                                                                                                  
ccGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGGACGGCGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGGGGCGGTAATGCCTtc                                                                                                                                                >  2:361530/1‑149 (MQ=255)
     gTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGATAATGTCTTCCTCCa                                                                                                                                           >  2:233528/1‑149 (MQ=255)
                            cGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGATAATGCCTTCCTCCAGACAGCGCACCACTTCGTTGAcc                                                                                                                    >  2:417961/1‑149 (MQ=255)
                                    tttGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGATAATGCCTTCCTCCAGACAGGGCACCACTTCGTTGACCATCGGGat                                                                                                            >  2:237380/1‑149 (MQ=255)
                                          gcaCCGAGGGGGTCCCGCCAGCGGAACGCCCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGATAATGCCTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGat                                                                                                            <  1:53880/141‑1 (MQ=255)
                                            aCCGAGGGTGTCCAGCCAGCGGAACGCGCGGCCGGGGAACGGTGTGAAGTCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGATAATGCCTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGc                                                                                                    <  2:405149/149‑1 (MQ=255)
                                                           ccagcGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGATAATGCCTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCgctca                                                                                           >  1:38369/1‑139 (MQ=255)
                                                                   aCGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGATCATGCCTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTg                                                                             >  2:39873/1‑149 (MQ=255)
                                                                                            gCCCAGGCCGTAGACCAGCGCCATATCCCCCTCCGCCGGAGTGGCGATAATGCCTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCggctggc                                                          <  1:160198/143‑1 (MQ=255)
                                                                                             cccAGGCCTGAAACCAGCGCCACTTCCGCTTCCGCCGGAGTGGCGATAATGCCTCCCTCCAGATAGCGCACCACTTCTTTGACCATCGTTATCCTCATGCGGGCGGTAATCTCTTCTTCGCTGAAATCGCGCTTCggctggct                                                         <  1:1990/143‑1 (MQ=255)
                                                                                                 ggCCGTAGACCAGCGCCATCTCCGCTTCCGCGCGGATGGACTTTCTGACTGCGTTTTCACAGCTCTCCCCTTCGTTGTCCATCGGTAGCTTCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACt                                                     <  1:272056/143‑1 (MQ=255)
                                                                                                     gTAGACCAGCGCCACAGCCGATTACGCCGGAGTGGCGATAATGTATTCCTCAATACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCcagca                                           <  2:99624/149‑1 (MQ=255)
                                                                                                              cgcCCTGTCCGCTTCCGCCGGACTGGCGATAATGCCTTCCGCCAGACAGCGCACCACTTCGTCGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGt                                        <  1:326798/143‑1 (MQ=255)
                                                                                                               gcCATATCCGCTTCCGCCGGAGTGGCGATAATGTCTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATCATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTc                                       <  1:204126/143‑1 (MQ=255)
                                                                                                                               ccgGAGTGGCGATAATGCCTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGtcttcttct                 <  2:227615/149‑1 (MQ=255)
                                                                                                                                        gCGATAATGCCTTCCTCCAGACAGCACACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCGTCAACGGCGGCGTCTTCTTCTTTCTTCgg         <  2:263781/149‑1 (MQ=255)
                                                                                                                                          gATAATGTCTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCCGCAGGTCTTCAACGGCGGCGTCTTCttctttct             >  1:79537/1‑143 (MQ=255)
                                                                                                                                          gATAATGTCTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCttctttct             >  1:269003/1‑143 (MQ=255)
                                                                                                                                               tGCCTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATACGGGCGATACTCTCTGCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCATCACCGGCGGCGACTTCTCCTTTCTTCGGCTTACCt  >  2:338972/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGATAATGCCTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTTACCT  >  CP009273/4022251‑4022541

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: