Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,931,037 C→A 71.5% 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,0370CA71.5% 38.9 / 47.0 51A212A (GCG→GCTfucOL‑1,2‑propanediol oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (5/9);  new base A (23/13);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

CCTCGCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATCGACACCCGTCGCAGCTT  >  W3110S.gb/2930981‑2931104
                                                        |                                                                   
ccTCGCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACg                                                       >  1:3310600/1‑71 (MQ=255)
ccTCGCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACg                                                       >  1:1252615/1‑71 (MQ=255)
ccTCGCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACg                                                       >  1:2222252/1‑71 (MQ=255)
 cTCGCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTA‑‑GCCCACGCGCCACggg                                                     <  1:1375370/70‑1 (MQ=255)
    gCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACggg                                                     <  1:3225218/69‑1 (MQ=255)
    gCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACggg                                                     <  1:760848/69‑1 (MQ=255)
     cAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTaa                                                  <  1:577880/71‑1 (MQ=255)
     cAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTaa                                                  <  1:3300873/71‑1 (MQ=255)
      aGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTaa                                                  <  1:15897/70‑1 (MQ=255)
       gcgcCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAata                                                <  1:313933/71‑1 (MQ=255)
       gcgcCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAata                                                <  1:2388801/71‑1 (MQ=255)
       gcgcCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAata                                                <  1:2899788/71‑1 (MQ=255)
       gcgcCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAata                                                <  1:1227382/71‑1 (MQ=255)
       gcgcCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAata                                                <  1:2658519/71‑1 (MQ=255)
         gcCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAAtata                                              <  1:3420735/71‑1 (MQ=255)
          ccccAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAATATAc                                             <  1:898041/71‑1 (MQ=255)
          ccccAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAATATAc                                             <  1:1238623/71‑1 (MQ=255)
            ccAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAATATAccc                                           <  1:3648757/71‑1 (MQ=255)
                  aTGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAATATAcccc                                          >  1:1108046/1‑66 (MQ=255)
                  aTGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAATATAcccc                                          >  1:3374218/1‑66 (MQ=255)
                  aTGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAATATAcccc                                          >  1:407022/1‑66 (MQ=255)
                  aTGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAATATAcccc                                          >  1:952538/1‑66 (MQ=255)
                  aTGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCAAGCGCCACGGGTAATATAcc                                            >  1:3439360/1‑64 (MQ=255)
                  aTGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCAACGCGCCACGGGTAATATAcccc                                          >  1:2366300/1‑66 (MQ=255)
                   tGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTaa                                                  <  1:3107071/57‑1 (MQ=255)
                   tGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTaa                                                  <  1:2049225/57‑1 (MQ=255)
                   tGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAata                                                >  1:296791/1‑59 (MQ=255)
                   tGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAata                                                >  1:3293544/1‑59 (MQ=255)
                   tGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAata                                                >  1:3378020/1‑59 (MQ=255)
                   tGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAata                                                >  1:2406600/1‑59 (MQ=255)
                   tGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAata                                                >  1:345785/1‑59 (MQ=255)
                   tGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAata                                                >  1:471855/1‑59 (MQ=255)
                        tCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCCCGCGCCAAGGGTAATATACCCCTCAATAGCAt                                >  1:491063/1‑70 (MQ=255)
                        tCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAATATACCCCTCAATAGCAt                                >  1:2439670/1‑70 (MQ=255)
                        tCAATCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAATATACCCCTCAATAGCAt                                >  1:3722688/1‑70 (MQ=255)
                           aTCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAATATACCCCTCAATAGCAtgagt                            <  1:1453550/71‑1 (MQ=255)
                           aTCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAATATACCCCTCAATAGCAtgagt                            <  1:2732324/71‑1 (MQ=255)
                           aTCGCTTTAATGTGCAGTGCATCGGTTAGAGCCCACGCGCCACGGGTAATATACCCCTCAATAGCAtgagt                            <  1:446466/71‑1 (MQ=255)
                                tttAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACTGGTaa                                                  <  1:974880/44‑1 (MQ=255)
                                tttAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTaa                                                  <  1:988536/44‑1 (MQ=255)
                                  tAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAata                                                >  1:3684598/1‑44 (MQ=255)
                                  tAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAata                                                >  1:1484203/1‑44 (MQ=255)
                                  tAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAata                                                >  1:3242593/1‑44 (MQ=255)
                                       tgcagtgcaTCGGTTAGAGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATCGa                 <  1:2882336/70‑1 (MQ=255)
                                              caTCGGTTAGAGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATCGACACCCGTc         >  1:2779643/1‑71 (MQ=255)
                                              caTCGGTTAGAGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATCGACACCCGTc         >  1:889633/1‑71 (MQ=255)
                                                  ggTTAGAGCCCACGCGCCACGGGTAATATACCCCTCAATAg                                   >  1:237907/1‑41 (MQ=255)
                                                  ggTTAGAGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAgcg                       >  1:2825925/1‑53 (MQ=255)
                                                  ggTTAGAGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAgcg                       >  1:1026012/1‑53 (MQ=255)
                                                     tAGCGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATCGACACCCGTCGCAGCtt  >  1:2399814/1‑71 (MQ=255)
                                                     tAGCGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATCGACACCCGTCGCAGCtt  >  1:2268418/1‑71 (MQ=255)
                                                        |                                                                   
CCTCGCAGCGCCCCAGCAATGATTTCAATCGCTTTAATGTGCAGTGCATCGGTTAGCGCCCACGCGCCACGGGTAATATACCCCTCAATAGCATGAGTGAGCGCATCGACACCCGTCGCAGCTT  >  W3110S.gb/2930981‑2931104

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: