Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,028,910 G→C 11.1% Q688E (CAG→GAG)  fadB ← dodecenoyl‑CoA delta‑isomerase, enoyl‑CoA hydratase, 3‑hydroxybutyryl‑CoA epimerase, 3‑hydroxyacyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,028,9100GC11.1% 93.4 / 5.9 37Q688E (CAG→GAG) fadBdodecenoyl‑CoA delta‑isomerase, enoyl‑CoA hydratase, 3‑hydroxybutyryl‑CoA epimerase, 3‑hydroxyacyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (20/12);  new base C (2/2);  total (23/14)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

CATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGATAATGCCTTCCTCCAGACAGCG  >  NC_000913/4028770‑4029051
                                                                                                                                            |                                                                                                                                             
cATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCCACCGGCCGGGCTGGCTAAACCAGCGGAGAGTACGGTTCGGTATGGCGCGCCTTATTACGCCGAGCGTTCGGCACTTCGGACAGCGGGCGGGGGTGCTGATATTGCtgtg                                                                                                                                             >  1:165648/1‑143 (MQ=255)
          ccTTCAGCCGGTTTCAGGTGTCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAg                                                                                                                                   <  2:868681/143‑1 (MQ=255)
                 ccGTTCTCAGTTCGCCACCCGGACGGGCTGGCTCAACCGGACCAGAGTCCGTTTCGTTATGACGCGCTTTATTACGCAGTCCTCCCGGCACTCCATACCGCGGGCCGTGGTGCTGATATTGCTGTGCCATATCGAGGTAtttt                                                                                                                            <  2:143146/143‑1 (MQ=255)
                  cGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGGCGCGGGGTATTGCGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGGAttttt                                                                                                                           >  1:623945/1‑142 (MQ=255)
                          ggTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTTTGCCCTCTCGAGCTCTTTTGCCCTcccc                                                                                                                   >  1:598135/1‑142 (MQ=255)
                           gTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGGGCTGATATTGCTGTGCCATCTCGAGGTACTTTGGGCTACCgt                                                                                                                  >  2:836519/1‑142 (MQ=255)
                            tccccgaCCGGACGTGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAg                                                                                                                 <  2:212832/137‑1 (MQ=255)
                            tCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAg                                                                                                                 >  1:102016/1‑143 (MQ=255)
                                     ggACGGGCTGACTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTccag                                                                                                        <  1:540489/143‑1 (MQ=255)
                                                 tCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGAAGTTTGCGCTGCCGAGGGTGTCCAGCCAGCGGAgcgc                                                                                            >  1:452262/1‑143 (MQ=255)
                                                  cAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAAcgcg                                                                                           <  2:551062/143‑1 (MQ=255)
                                                        gaggaTAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGCATTTTGCGCTCCCGAGGGTGTGCAGCCACCGCAGCGCGccgccg                                                                                     >  2:124977/1‑143 (MQ=255)
                                                                tACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGACCgg                                                                             >  2:256776/1‑143 (MQ=255)
                                                                tACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTCTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACgg                                                                             >  1:812857/1‑143 (MQ=255)
                                                                gaCGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTCTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACgg                                                                             >  1:805600/2‑143 (MQ=255)
                                                                 aCGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAgg                                                                                                                <  2:661987/107‑1 (MQ=255)
                                                                 aCGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAgg                                                                                                                >  1:661987/1‑107 (MQ=255)
                                                                   ggTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATGGAAGTATTTTGCGCTACCGACGAGGGGCAGCCCGCGGAACGCGACGCCGGGGAACGGAgg                                                                          >  2:892788/1‑143 (MQ=255)
                                                                   ggTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGacg                                                                          >  1:466622/1‑141 (MQ=255)
                                                                    gTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTCTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAggg                                                                         <  2:812857/143‑1 (MQ=255)
                                                                    gTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTCTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAggg                                                                         <  2:805600/143‑1 (MQ=255)
                                                                      tCGTGATGACGCGGTTGATGACGCGGGCCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGaa                                                                       <  2:461699/143‑1 (MQ=255)
                                                                      ccGTTATGGCGCGCTTTATGACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGaa                                                                       <  2:828639/142‑1 (MQ=255)
                                                                              aCGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCTAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGc                                                               <  1:502868/143‑1 (MQ=255)
                                                                                gcgcTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCCGCCCGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCg                                                             >  1:36000/1‑143 (MQ=255)
                                                                                          ccGCAGCCCTCCCGGCACTTCGTACAGCGGGCCGAGGTGCCGAGATTGCTGGGCCATATCGGGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAgcg                                                   <  2:763177/142‑1 (MQ=255)
                                                                                          aggcggACCTTCCGGCCCTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAgcg                                                   <  2:682538/138‑1 (MQ=255)
                                                                                              aGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGACGCCCCGCCCGCAGACCAGAGccct                                               >  2:176952/1‑141 (MQ=255)
                                                                                                          aCTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGACCGCGCCGCCGTGGAACGGAGGGCAGCCCAGGCCGTAGTCCCGCGCCATATCCGCTTccgc                                   >  1:252967/1‑143 (MQ=255)
                                                                                                            ttCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTGTTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCGGCGGCATGACCGCGTccgccg                                 >  2:474147/1‑143 (MQ=255)
                                                                                                            ttCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATACCCGCTTccgccg                                 >  2:435765/1‑143 (MQ=255)
                                                                                                                      gggCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCCGCCCGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCCGGCCGTCGCCTCGCGCCGTATGTGCTTCCGCCGGTgcggcgtt                       >  1:232280/1‑141 (MQ=255)
                                                                                                                         ccGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTCGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTCCCGCCGGAGTGGCGATAAt                    >  2:228092/1‑143 (MQ=255)
                                                                                                                             ggTGCTGATCTTGCTGTGCCATATCGAGGTACTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGATAATGCCt                <  1:193802/143‑1 (MQ=255)
                                                                                                                                gCTGATATTGCTGTGCCATATCGAGGAATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGATAATGCCTtcc             >  2:528564/1‑143 (MQ=255)
                                                                                                                                 cTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGAGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGAGAATGCCTtccg            >  1:167080/1‑142 (MQ=255)
                                                                                                                                     tatTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGATAATGCCTTCCTCCAg        >  1:633765/1‑143 (MQ=255)
                                                                                                                                         gCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGATAATGCCTTCCTccagccag    >  1:11018/1‑143 (MQ=255)
                                                                                                                                           tgtgCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGATAATGCCTTCCTCCAGACAgcg  >  1:560380/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                             
CATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATATTGCTGTGCCATATCGAGGTATTTTGCGCTACCGAGGGTGTCCAGCCAGCGGAACGCGCCGCCGTGGAACGGAGGGAAGCCCAGGCCGTAGACCAGCGCCATATCCGCTTCCGCCGGAGTGGCGATAATGCCTTCCTCCAGACAGCG  >  NC_000913/4028770‑4029051

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: