Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,077,892 G→A 16.2% R171R (CGC→CGT glpD ← sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,077,8920GA16.2% 83.6 / 8.8 37R171R (CGC→CGTglpDsn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (18/13);  new base A (3/3);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGCCGGTATCGATATCTTCCGCTTCCACAATCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACGCACCACCATCTGGGCGTTGGCGAGTACCAGACGGGCG  >  W3110S.gb/4077826‑4077961
                                                                  |                                                                     
tttGCCGGTATCGATATCTTCCGCTTCCACAATCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGa                                                                    <  1:251116/70‑1 (MQ=255)
    ccGGTATCGATATCTTCCGCTTCCACAATCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTcc                                                                <  1:144015/70‑1 (MQ=255)
          tCGATATCTTCCGCTTCCACAATCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAg                                                         >  1:2633078/1‑71 (MQ=255)
          tCGATATCTTCCGCTTCCACAATCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAg                                                         >  1:2941163/1‑71 (MQ=255)
            gATATCTTCCGCTTCCACAATCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCa                                                       >  1:2415370/1‑71 (MQ=255)
            gATATCTTCCGCTTCCACAATCCACAGGCCGTTTTCGAGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCa                                                       >  1:3555811/1‑71 (MQ=255)
              tatCTTCCGCTTCCACAATCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACt                                                     <  1:815026/71‑1 (MQ=255)
                tCTTCCGCTTCCACAATCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACtt                                                    >  1:163556/1‑70 (MQ=255)
                  ttCCGCTTCCACAATCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTcac                                                 >  1:3604815/1‑69 (MQ=255)
                    ccGCTTCCACAATCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTcgccg                                               >  1:1650660/1‑71 (MQ=255)
                     cGCTTCCACAATCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTcgccgc                                              <  1:1452180/71‑1 (MQ=255)
                     cGCTTCCACAATCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTcgccgc                                              <  1:391246/71‑1 (MQ=255)
                          ccACAATCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCtt                                           <  1:3646184/69‑1 (MQ=255)
                           cacaATCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTCTAc                                        >  1:262292/1‑71 (MQ=255)
                             caATCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACg                                       >  1:3268436/1‑70 (MQ=255)
                                tCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACt                                                     <  1:1078907/53‑1 (MQ=255)
                                  cacaGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACGcaccac                                 >  1:1074314/1‑71 (MQ=255)
                                   acaGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACGCaccacc                                <  1:34800/71‑1 (MQ=255)
                                    caGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTa                                         >  1:2806608/1‑61 (MQ=255)
                                    caGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACGCAccacca                               <  1:887942/71‑1 (MQ=255)
                                        ccGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACGCACCACCATCTg                           <  1:1391692/71‑1 (MQ=255)
                                           ttttCGGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACGCACCACCAt                              <  1:2117664/64‑1 (MQ=255)
                                             ttCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACGCACCACCATCTGGGCg                       >  1:3555413/1‑70 (MQ=255)
                                             ttCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACGCACCACCATCTGGGCGt                      >  1:28418/1‑71 (MQ=255)
                                             ttCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACGCACCACCATCTGGGCGt                      >  1:146989/1‑71 (MQ=255)
                                             ttCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACGCACCACCATCTGGGCGt                      >  1:1449394/1‑71 (MQ=255)
                                                       gCAGAGGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACGCACCACCATCTGGGCGTTGGCGAGTAc            >  1:1019120/1‑71 (MQ=255)
                                                          gagGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACGCACCACCATCTGGGCGTTGGCGAGTACCAg         <  1:646686/71‑1 (MQ=255)
                                                           acgTGGCACGAGTACGAGTAAGCACTTCACCGCCTTTACGCACCACCATCTGGGCGTTgg                   >  1:2800053/3‑60 (MQ=255)
                                                           acgTGGCACGAGTACGAGTAAGCACTTCACCGCCTTTACGCACCACCATCTGGGCGTTgg                   >  1:65758/3‑60 (MQ=255)
                                                           acgTGGCACGAGTACGAGTAAGCACTTCACCGCCTTTACGCACCACCATCTGGGCGTTgg                   >  1:2585903/3‑60 (MQ=255)
                                                             gTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACGCACCACCATCTGGGCGTTGGCGAGTACCAGACg      >  1:2353448/1‑71 (MQ=255)
                                                             gTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACGCACCACCATCTGGGCGTTGGCGAGTACCAGACg      <  1:2609874/71‑1 (MQ=255)
                                                             gTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACGCACCACCATCTGGGCGTTGGCGAGTACCAGACg      >  1:3258209/1‑71 (MQ=255)
                                                              tGGCACGAGTACGAGTAAGCACTTCACCGCCTTTACGCACCACCATCTGGGCGTTGGCGAGTACCAGAc       <  1:191594/69‑1 (MQ=255)
                                                              tGGCACGAGTACGAGTAAGCACTTCACCGCCTTTACGCACCACCATCTGGGCGTTGGCGAGTACCAGAc       <  1:1981201/69‑1 (MQ=255)
                                                              tGGCACGAGTACGAGTAAGCACTTCACCGCCTTTACGCACCACCATCTGGGCGTTGGCGAGTACCAGAc       <  1:2824523/69‑1 (MQ=255)
                                                                  gcgAGTCCGAGTAAGCACTTCGCCGCCTTTACGCACCACCATCTGGGCGTTGGCGAGTACCAGACGGGcg  >  1:3721202/1‑70 (MQ=255)
                                                                  |                                                                     
TTTGCCGGTATCGATATCTTCCGCTTCCACAATCCACAGGCCGTTTTCGCGGCGAGCAGAGGTGGCGCGAGTCCGAGTAAGCACTTCGCCGCCTTTACGCACCACCATCTGGGCGTTGGCGAGTACCAGACGGGCG  >  W3110S.gb/4077826‑4077961

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: