Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,608,041 C→T 20.0% L45L (CTG→TTG)  fadA → 3‑ketoacyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,608,0410CT20.0% 62.3 / 9.8 30L45L (CTG→TTG) fadA3‑ketoacyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base C (10/14);  new base T (4/2);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

TTTCGTAACGTGCGTGCAGAAGATCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTTGTGTGCAGCAGACGCTGGAGCAGGGTTTT  >  W3110S.gb/3607972‑3608106
                                                                     |                                                                 
tttCGTAACGTGCGTGCAGAAGATCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCt                                                                  >  1:3015844/1‑71 (MQ=255)
  tCGTAACGTGCGTGCAGAAGATCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTg                                                                 >  1:2459475/1‑70 (MQ=255)
       aCGTGCGTGCAGAAGATCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAgc                                                            <  1:1061413/70‑1 (MQ=255)
         gtgcgtgcAGAAGATCTCTCCTCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGTAAgcggc                                                         <  1:4497519/71‑1 (MQ=255)
           gcgtgcAGAAGATCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTg                                                                 <  1:2738477/61‑1 (MQ=255)
                cAGAAGATCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTc                                                  >  1:283006/1‑71 (MQ=255)
                 agaagaTCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCt                                                                  <  1:3561007/54‑1 (MQ=255)
                    agaTCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTcgac                                               >  1:2804969/1‑70 (MQ=255)
                     gaTCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGTTGGaa                                                              >  1:3868400/1‑54 (MQ=255)
                     gaTCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTg                                                                 >  1:3648978/1‑51 (MQ=255)
                     gaTCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCt                                                   >  1:2494119/1‑65 (MQ=255)
                       tctctcCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGTTGGa                                                               >  1:4011686/1‑51 (MQ=255)
                        ctctcCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGatat                                          >  1:498301/1‑71 (MQ=255)
                           tcCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTcgac                                               >  1:3893453/1‑63 (MQ=255)
                           tcCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGata                                           <  1:4245317/67‑1 (MQ=255)
                           tcCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCGGGAAGCGGCGGCCCTCGACGat                                            <  1:3116646/66‑1 (MQ=255)
                             cGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTa                                       <  1:714565/69‑1 (MQ=255)
                              gCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTg                                    <  1:4061991/71‑1 (MQ=255)
                               cTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGTTGGAAGCGGCGGCCCTCGCCGATAATTAc                                      >  1:2331639/1‑68 (MQ=255)
                                tCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGTTGGAAGCGGCGGCCCTCGACGATATTTa                                       >  1:3656793/1‑66 (MQ=255)
                                   tttAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGata                                           <  1:3564362/59‑1 (MQ=255)
                                               ccTGCTGGCGCGTAACCCGGCGTTGGAAGCGGCGGCCCTcgacga                                             <  1:714096/45‑1 (MQ=255)
                                               ccTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTtg                              <  1:3260021/60‑1 (MQ=255)
                                               ccTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTTGTGTgcagca                     <  1:1955694/69‑1 (MQ=255)
                                                   ctgGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCt                                                   >  1:3635179/1‑35 (MQ=255)
                                                   ctgGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTTGTGTGCAGCAGACGc                <  1:628907/70‑1 (MQ=255)
                                                      gcgcGTAACCCGGCGTTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTTGTGTGCAGCAGACGCTg              <  1:1456796/69‑1 (MQ=255)
                                                      gcgcGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTtg                              <  1:2422529/53‑1 (MQ=255)
                                                            aaCCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTTGTGTGCAGCAGACGCTGGAGCAggg      >  1:4158760/1‑71 (MQ=255)
                                                                cGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTTGTGTGCAGCAGACGCTGGAGCAGGGtttt  <  1:833593/71‑1 (MQ=255)
                                                                     |                                                                 
TTTCGTAACGTGCGTGCAGAAGATCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTTGTGTGCAGCAGACGCTGGAGCAGGGTTTT  >  W3110S.gb/3607972‑3608106

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: