Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,931,037 C→A 42.7% A212A (GCG→GCT fucO ← L‑1,2‑propanediol oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,931,0370CA42.7% 53.1 / 34.5 42A212A (GCG→GCTfucOL‑1,2‑propanediol oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (8/16);  new base A (10/8);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAACCGATCCTCGCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATCGACACCCGTCGC  >  W3110S.gb/2930971‑2931099
                                                                  |                                                              
aGCAACCGATCCTCGCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAgcg                                                               >  1:591560/1‑68 (MQ=255)
    aCCGATCCTCGCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCAc                                                          >  1:1594699/1‑69 (MQ=255)
    aCCGATCCTCGCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCAc                                                          >  1:2930313/1‑69 (MQ=255)
       gATCCTCGCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCAcgc                                                        >  1:85722/1‑68 (MQ=255)
         tCCTCGCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGcc                                                     <  1:2882866/69‑1 (MQ=255)
          ccTCGCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCa                                                    >  1:459576/1‑69 (MQ=255)
              gCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCAcg                                                         >  1:2084911/1‑60 (MQ=255)
                 gcgcCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTaa                                             <  1:3516863/69‑1 (MQ=255)
                 gcgcCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTaa                                             <  1:2724370/69‑1 (MQ=255)
                  cgcCCCAGCAATGATTTCTATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTaa                                             <  1:1952901/68‑1 (MQ=255)
                  cgcCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTaa                                             <  1:1188905/68‑1 (MQ=255)
                   gcCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAata                                           <  1:2733541/69‑1 (MQ=255)
                   gcCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAata                                           <  1:2784013/69‑1 (MQ=255)
                     cccAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAAtata                                         <  1:1898950/69‑1 (MQ=255)
                             tGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAata                                           <  1:1411923/59‑1 (MQ=255)
                             tGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACgg                                                 >  1:1699058/1‑53 (MQ=255)
                               aTTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTaa                                             <  1:3057277/55‑1 (MQ=255)
                               aTTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTaa                                             <  1:3040376/55‑1 (MQ=255)
                               aTTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTaa                                             <  1:2959468/55‑1 (MQ=255)
                               aTTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTaa                                             <  1:3601739/55‑1 (MQ=255)
                               aTTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAATATACCCCTCAAt                                >  1:3146100/1‑68 (MQ=255)
                               aTTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAATATACCCCTCAAt                                >  1:1452102/1‑68 (MQ=255)
                               aTTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAATATACCCCTCAAt                                >  1:525965/1‑68 (MQ=255)
                                   cAATCGCTTTAATGTGCAGTGCATCGGTTAGAGccc                                                            >  1:2090179/1‑36 (MQ=255)
                                   cAATCGCTTTAATGTGCAGTGCATCGGTTAGAGccc                                                            >  1:1093452/1‑36 (MQ=255)
                                   cAATCGCTTTAATGTGCAGTGCATCGGTTAGAGccc                                                            >  1:2486390/1‑36 (MQ=255)
                                      tCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACgcgc                                                      <  1:1821991/39‑1 (MQ=255)
                                      tCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACgcgc                                                      <  1:2980484/39‑1 (MQ=255)
                                        gCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAATATACCCCTCAATAGCAtgagtg                      >  1:1840548/1‑69 (MQ=255)
                                          tttAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAata                                           <  1:3238426/46‑1 (MQ=255)
                                             aaTGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTaa                                             <  1:674685/41‑1 (MQ=255)
                                             aaTGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTaa                                             <  1:1133895/41‑1 (MQ=255)
                                             aaTGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTaa                                             <  1:2625540/41‑1 (MQ=255)
                                             aaTGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTaa                                             <  1:2248439/41‑1 (MQ=255)
                                              aTGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCa                >  1:2395928/1‑69 (MQ=255)
                                                 tgcagtgcaTCGGTTAGAGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATc              <  1:2403558/68‑1 (MQ=255)
                                                   cagtgcaTCGGTTAGCGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATCGa            <  1:1536258/68‑1 (MQ=255)
                                                       gcaTCGGTTAGCGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATCGACAccc       >  1:2912833/1‑69 (MQ=255)
                                                       gcaTCGGTTAGCGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATCGACAccc       <  1:53344/69‑1 (MQ=255)
                                                       gcaTCGGTTAGCGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATCGACAccc       <  1:604085/69‑1 (MQ=255)
                                                       gcaTCGGTTAGCGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATAGACAccc       >  1:213360/1‑69 (MQ=255)
                                                            ggTTAGAGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATCGACACCCGTCGc  >  1:2909490/1‑69 (MQ=255)
                                                            ggTTAGAGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATCGACACCCGTCGc  >  1:369495/1‑69 (MQ=255)
                                                                  |                                                              
AGCAACCGATCCTCGCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATCGACACCCGTCGC  >  W3110S.gb/2930971‑2931099

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: