Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,558,086 G→A 58.1% A168A (GCG→GCA uxuB → D‑mannonate oxidoreductase, NAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,558,0860GA58.1% 9.9 / 36.4 31A168A (GCG→GCAuxuBD‑mannonate oxidoreductase, NAD‑binding
Reads supporting (aligned to +/- strand):  ref base G (8/5);  new base A (10/8);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

CTGATCAAGCACGATCTGGAAAACCCGACTGCGCCGAAGTCCGCGATTGGTTACATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTG  >  W3110S.gb/4558021‑4558150
                                                                 |                                                                
cTGATCAAGCACGATCTGGAAAACCCGACTGCGCCGAAGTCCGCGATTGGTTACATCGTCGAAGCATtgcg                                                             <  1:2615037/71‑1 (MQ=255)
   aTCAAGAACGATCTGGAAAACCCGACTGCGCCGAAGTCCGCGATTGGTTACATCGTCGAAGCGTtgcgtct                                                          >  1:135887/1‑71 (MQ=255)
      aaGCACGATCTGGAAAACCCGACTGCGCCGAAGTCCGCGATTGGTTACATCGTCGAAGCATtgcgtctgcg                                                       >  1:1444911/1‑71 (MQ=255)
                              gcgcCGAAGTCCGCGATTGGTTACATCGTCGAAGCATTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGtt                               >  1:2630356/1‑71 (MQ=255)
                               cgcCGAAGTCCGCGATTGGTTACATCGTCGAAGCATTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGtt                               <  1:2410529/70‑1 (MQ=255)
                                 ccGAAGTCCGCGATTGGTTACATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGa                                      <  1:837178/61‑1 (MQ=255)
                                    aaGTCCGCGATTGGTTACATCGTCGAAGCGTTGCGTCTGCGTCGTGaaaa                                              <  1:1266921/50‑1 (MQ=255)
                                      gTCCGCGATTGGTTACATCGTCGAAGCATTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGt                         <  1:1236235/69‑1 (MQ=255)
                                      gTCCGCGATTGGTTACATCGTCGAAGCATTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGt                         <  1:670400/69‑1 (MQ=255)
                                      gTCCGCGATTGGTTACATCGTCGAAGCATTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGa                       >  1:1272108/1‑71 (MQ=255)
                                       tCCGCGATTGGTTACATCGTCGAAGCATTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGt                         <  1:2563088/68‑1 (MQ=255)
                                       tCCGCGATTGGTTACATCGTCGAAGCATTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGa                       <  1:2647429/70‑1 (MQ=255)
                                         cgcgATTGGTTACATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGtt                               >  1:1303352/1‑60 (MQ=255)
                                         cgcgATTGGTTACATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGt                    >  1:94396/1‑71 (MQ=255)
                                         cgcgATTGGTTACATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGt                    <  1:2242476/71‑1 (MQ=255)
                                         cgcgATTGGTTACATCGTCGAAGCATTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATg                     >  1:2174658/1‑70 (MQ=255)
                                         cgcgATTGGTTACATCGTCGAAGCATTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATg                     >  1:1684861/1‑70 (MQ=255)
                                         cgcgATTGGTTACATCGTCGAAGCATTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGt                    >  1:2744832/1‑71 (MQ=255)
                                             aTTGGTTACATCGTCGAAGCATTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTg                <  1:1295503/71‑1 (MQ=255)
                                                ggTTACATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACgg                          >  1:668230/1‑58 (MQ=255)
                                                ggTTACATCGTCGAAGCATTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGa                       >  1:2390725/1‑61 (MQ=255)
                                                    aCATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGa                       >  1:2401715/1‑57 (MQ=255)
                                                    aCATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGa                       >  1:1372143/1‑57 (MQ=255)
                                                    aCATCGTCGAAGCATTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATaa          <  1:2622213/70‑1 (MQ=255)
                                                       tcgtcgAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGa                       <  1:1709041/54‑1 (MQ=255)
                                                       tcgtcgAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGa                       <  1:1537723/54‑1 (MQ=255)
                                                         gtcgAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATAAcgtgcg    >  1:238852/1‑71 (MQ=255)
                                                         gtcgAAGCATTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGt                                >  1:2242090/1‑43 (MQ=255)
                                                         gtcgAAGCATTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGt                                >  1:1005276/1‑43 (MQ=255)
                                                         gtcgAAGCATTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGa                       >  1:1484506/1‑52 (MQ=255)
                                                            gAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATAAcgtgcgtg  >  1:1301756/1‑70 (MQ=255)
                                                                 |                                                                
CTGATCAAGCACGATCTGGAAAACCCGACTGCGCCGAAGTCCGCGATTGGTTACATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTG  >  W3110S.gb/4558021‑4558150

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