Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,608,041 C→T 23.3% 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,0410CT23.3% 57.6 / 12.3 30L45L (CTG→TTG) fadA3‑ketoacyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base C (9/14);  new base T (4/3);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAACGTGCGTGCAGAAGATCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTTGTGTGCAGCAGACGCTGGAGCAGGGTTTT  >  W3110S.gb/3607976‑3608106
                                                                 |                                                                 
gTAACGTGCGTGCAGAAGATCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGTTGGa                                                               <  1:2160912/70‑1 (MQ=255)
            cAGAAGATCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTc                                                  >  1:3328773/1‑71 (MQ=255)
               aagaTCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCgtcggc                                                      <  1:3534803/64‑1 (MQ=255)
               aagaTCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCggcggc                                                      <  1:3676/64‑1 (MQ=255)
               aagaTCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCggcggc                                                      <  1:3483613/64‑1 (MQ=255)
               aagaTCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTcga                                                >  1:737282/1‑70 (MQ=255)
               aagaTCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTcga                                                >  1:170776/1‑70 (MQ=255)
                 gaTCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTcgacga                                             <  1:1880121/71‑1 (MQ=255)
                  aTCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGTTGGAAGCGGCGGCCCTcgac                                               <  1:2409307/68‑1 (MQ=255)
                    ctctcCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGTTGGAAGCGGCGGCCCTcgacga                                             >  1:69941/1‑68 (MQ=255)
                    ctctcCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGTTGGAAGCGGCGGCCCTcgacga                                             >  1:2823385/1‑68 (MQ=255)
                    ctctcCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGata                                           >  1:3548611/1‑70 (MQ=255)
                       tcCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGat                                            <  1:1240702/66‑1 (MQ=255)
                          gCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTg                                    <  1:2487924/71‑1 (MQ=255)
                            tCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATAtt                                         <  1:2974019/64‑1 (MQ=255)
                            tCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATAtt                                         <  1:1223665/64‑1 (MQ=255)
                            tCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATAtt                                         <  1:563950/64‑1 (MQ=255)
                            tCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATAtt                                         <  1:330844/64‑1 (MQ=255)
                              aTTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATAttt                                        <  1:2643189/63‑1 (MQ=255)
                                 tAATGCTTAGCCTGCTGGCGCGTAACCCGGCGTTGGAAg                                                             >  1:522379/1‑39 (MQ=37)
                                 tAATGCGTAGCCTGCTGGCGCGTAACCCGGCGTTGGAAg                                                             >  1:1399357/1‑39 (MQ=255)
                                 tAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTtgt                             >  1:364680/1‑71 (MQ=255)
                                 tAAAGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGCCGATATTTACTGGGGTtgt                             >  1:359700/1‑71 (MQ=255)
                                               ctgGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTTGTGTGCAGCAGACGCt               <  1:3540529/71‑1 (MQ=255)
                                               ctgGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTTGTGTGCAGCAGACGCt               <  1:2759850/71‑1 (MQ=255)
                                                  gcgcGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTa                                       >  1:4020199/1‑44 (MQ=255)
                                                      gTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTTGTGTGCAGCAGACGCTGGAGCAg        >  1:2494911/1‑71 (MQ=255)
                                                            cGGCGTTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTTGTGTGCAGCAGACGCTGGAGCAGGGtttt  <  1:3545297/71‑1 (MQ=255)
                                                            cGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTTGTGTgca                        >  1:2323055/1‑49 (MQ=255)
                                                            cGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTTGTGTGCAGCAGACGCTGGAGCAGGGtttt  <  1:1775667/71‑1 (MQ=255)
                                                                 |                                                                 
GTAACGTGCGTGCAGAAGATCTCTCCGCTCATTTAATGCGTAGCCTGCTGGCGCGTAACCCGGCGCTGGAAGCGGCGGCCCTCGACGATATTTACTGGGGTTGTGTGCAGCAGACGCTGGAGCAGGGTTTT  >  W3110S.gb/3607976‑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: