Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,791,255 G→A 45.1% G443S (GGC→AGC)  gabD → succinate‑semialdehyde dehydrogenase I, NADP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,791,2550GA45.1% 8.4 / 34.7 31G443S (GGC→AGC) gabDsuccinate‑semialdehyde dehydrogenase I, NADP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (12/5);  new base A (10/4);  total (22/9)
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

CGTGATTTAAGCCGCGTCTTCCGCGTGGGCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCGGCATTATTTCCAATGAAGTGGCCCCGTTCGGCGGCATCAAAGCCTCGGGTCTGGGTCGTGA  >  W3110S.gb/2791189‑2791316
                                                                  |                                                             
cGTGATTTAAGCCGCGTCTTCCGCGTGGGCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCGGCat                                                           <  1:156431/71‑1 (MQ=255)
  tGATTTAAGCCGCGTCTTCCGCGTGGGCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCGGCatta                                                         >  1:1792980/1‑71 (MQ=255)
               gTCTTCCGCGTGGGCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCAGCATTATTTCCAATGAAg                                             >  1:491780/1‑70 (MQ=255)
               gTCTTCCGCGTGGGCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCAGCATTATTTCCAATGAAg                                             >  1:3970782/1‑70 (MQ=255)
               gTCTTCCGCGTGGGCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCAGCATTATTTCCAATGAAg                                             >  1:1910320/1‑70 (MQ=255)
               gTCTTCCGCGTGGGCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCAGCATTATTTCCAATGAAg                                             >  1:214142/1‑70 (MQ=255)
                 cTTCCGCGTGGGCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCAGCATTATTTCCAATg                                                >  1:1804328/1‑65 (MQ=255)
                 cTTCCGCGTGGGCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCAGCATTATTTCCAATGaa                                              >  1:947669/1‑67 (MQ=255)
                 cTTCCGCGTGGGCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCAGCATTATTTCCAATGaa                                              >  1:2322060/1‑67 (MQ=255)
                  ttCCGCGTGGGCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCGGCATTATTTCCAAt                                                 >  1:1496445/1‑63 (MQ=255)
                  ttCCGCGTGGGCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCGGCATTATTTCCAATGaa                                              >  1:3428722/1‑66 (MQ=255)
                    ccGCGTGGGCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCGGCATTATTTCCAATGAAGTGGc                                         <  1:706537/69‑1 (MQ=255)
                        gTGGGCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCGGCATTATTTCCAATGAAGTGGCCCCGtt                                   >  1:855190/1‑71 (MQ=255)
                           ggCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCAGCATTATTTCCAAt                                                 >  1:140962/1‑54 (MQ=255)
                           ggCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCAGCATTATTTCCAAt                                                 >  1:1448623/1‑54 (MQ=255)
                           ggCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCAGCATTATTTCCAAt                                                 >  1:1560682/1‑54 (MQ=255)
                             cGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCGGCATTATTTCCAATGAAg                                             >  1:2173038/1‑56 (MQ=255)
                             cGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCGGCATTATTTCCAATGAAg                                             >  1:107449/1‑56 (MQ=255)
                                            cGGCATCGTCGGCATCAATACCGGCATTATTTCCAATGAAg                                             >  1:3888236/1‑41 (MQ=255)
                                            cGGCATCGTCGGCATCAATACCAGCATTATTTCCAATCAAGTGGCCCCGTTCGGCGGCATCAAAGCCTCg                <  1:4026653/70‑1 (MQ=255)
                                             ggCATCGTCGGCATCAATACCGGCATTATTTCCAATGAAGTGGCCCCGTTCGGCGGCATCAAAGCCTCggg              <  1:2002874/71‑1 (MQ=255)
                                                   gtcgGCATCAATACCGGCATTATTTCCAATGAAGTGGCCCCGTTCGGCGGCATCAAAGCCTCGGGTCTgg         >  1:767039/1‑70 (MQ=255)
                                                   gtcgGCATCAATACCGGCATTATTTCCAATGAAGTGGCCCCGGTCGGCGGCATCAAAGCCTCGGGTCTggg        >  1:56846/1‑71 (MQ=255)
                                                      ggCATCAATACCGGCATTATTTCCAATGAAGTGGCCCCGTTCGGCGGCATCAAAGCCTCGGGTCTGGGTCg     >  1:3980308/1‑71 (MQ=255)
                                                          tCAATACCGGCATTATTTCCAATGAAGTGGCCCCGTTCGGCGGCATCAAAGCCTCGGGTCTGGGTCGTg   >  1:906657/1‑69 (MQ=255)
                                                          tCAATACCGGCATTATTTCCAATGAAGTGGCCCCGTTCGGCGGCATCAAAGCCTCGGGTCTGGGTCGTGa  <  1:1660186/70‑1 (MQ=255)
                                                          tCAATACCGGCATTATTTCCAATGAAGTGGCCCCGTTCGGCAGCATCAAAGCCTCGGGTCTGGGTCGTg   >  1:3957144/1‑69 (MQ=255)
                                                             aTACCAGCATTATTTCCAATGAAGTGGCCCCGTTCGGCGGCATCAAAGCCTCGGGTCTGGGTCGTGa  <  1:848804/67‑1 (MQ=255)
                                                              tACCGGCATTATTTCCAATGAAGTGGCCCCGTTCGGCGGCATCAAAGcc                   <  1:3184921/49‑1 (MQ=255)
                                                              tACCAGCATTATTTCCAATGAAGTGGCCCCGTTcggc                               <  1:369263/37‑1 (MQ=255)
                                                              tACCAGCATTATTTCCAATGAAGTGGCCCCGTTcggc                               <  1:4004859/37‑1 (MQ=255)
                                                                  |                                                             
CGTGATTTAAGCCGCGTCTTCCGCGTGGGCGAAGCGCTGGAGTACGGCATCGTCGGCATCAATACCGGCATTATTTCCAATGAAGTGGCCCCGTTCGGCGGCATCAAAGCCTCGGGTCTGGGTCGTGA  >  W3110S.gb/2791189‑2791316

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: