Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,151,925 A→C 21.2% P207P (CCA→CCC fabD → malonyl‑CoA‑[acyl‑carrier‑protein] transacylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,151,9250AC21.2% 43.1 / 17.7 33P207P (CCA→CCCfabDmalonyl‑CoA‑[acyl‑carrier‑protein] transacylase
Reads supporting (aligned to +/- strand):  ref base A (17/9);  new base C (2/5);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGCGCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGAAT  >  W3110S.gb/1151861‑1151994
                                                                |                                                                     
cgggcgCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGcc                                                                 <  1:130251/71‑1 (MQ=255)
    cgcAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACa                                                             <  1:3724683/71‑1 (MQ=255)
    cgcAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACa                                                             <  1:4671555/71‑1 (MQ=255)
          aCGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACca                                                                       >  1:2519904/1‑55 (MQ=255)
                gcTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGAAGCCGACAAACTGGCAGTAg                                                 <  1:3317719/71‑1 (MQ=255)
                gcTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGAAGCCGACAAACTGGCAGTAg                                                 <  1:3311985/71‑1 (MQ=255)
                   gCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACa                                                             >  1:4169102/1‑56 (MQ=255)
                         aCCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCg                                         <  1:250927/70‑1 (MQ=255)
                           cAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaaa                                      <  1:2442649/71‑1 (MQ=255)
                           cAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaa                                       >  1:2853910/1‑70 (MQ=255)
                              tGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATc                                   <  1:180685/71‑1 (MQ=255)
                              tGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATc                                   <  1:3049954/71‑1 (MQ=255)
                                aGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc                                 >  1:1617685/1‑71 (MQ=255)
                                aGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc                                 >  1:4212783/1‑71 (MQ=255)
                                   gTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCtt                              <  1:1744139/71‑1 (MQ=255)
                                    tACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCttt                             <  1:1434667/71‑1 (MQ=255)
                                      ccGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTaa                           >  1:4191669/1‑71 (MQ=255)
                                      ccGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTaa                           >  1:1911013/1‑71 (MQ=255)
                                      ccGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTaa                           >  1:412501/1‑71 (MQ=255)
                                        gTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACg                         <  1:386327/71‑1 (MQ=255)
                                        gTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGc                                          >  1:2919793/1‑54 (MQ=255)
                                              cTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc                                 >  1:2390446/1‑57 (MQ=255)
                                               tgtgCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCaa                   >  1:2056495/1‑70 (MQ=255)
                                                gtgCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCaa                   >  1:3939416/1‑69 (MQ=255)
                                                 tgCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAg                >  1:1799611/1‑71 (MQ=255)
                                                     cTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCttt                             <  1:3442940/54‑1 (MQ=255)
                                                     cTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCttt                             <  1:656823/54‑1 (MQ=255)
                                                     cTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTcc            >  1:4452597/1‑71 (MQ=255)
                                                     cTGAAGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTcc            >  1:1615507/1‑71 (MQ=255)
                                                        atgaAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGt               >  1:4400683/1‑65 (MQ=255)
                                                        atgaAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtgt         >  1:432010/1‑71 (MQ=255)
                                                         tgaAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc                                 >  1:4122144/1‑46 (MQ=255)
                                                          gaAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtgttg       >  1:4434840/1‑71 (MQ=255)
                                                            aaCCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAAt                                    >  1:741028/1‑40 (MQ=255)
                                                               cagcagCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGAAt  <  1:3753915/71‑1 (MQ=255)
                                                                |                                                                     
CGGGCGCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGAAT  >  W3110S.gb/1151861‑1151994

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: