Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,272,176 T→C 31.5% D146G (GAC→GGC)  garR ← tartronate semialdehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,272,1760TC31.5% 92.2 / 64.9 73D146G (GAC→GGC) garRtartronate semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base C (15/8);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGC  >  W3110S.gb/3272109‑3272242
                                                                   |                                                                  
aCCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGt                                                                    >  1:238844/1‑68 (MQ=255)
aCCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGt                                                                    >  1:2203027/1‑68 (MQ=255)
aCCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGt                                                                    >  1:1534804/1‑68 (MQ=255)
aCCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGt                                                                    >  1:506881/1‑68 (MQ=255)
    gCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGaa                                                                >  1:1445232/1‑68 (MQ=255)
     cACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGaa                                                                >  1:844486/1‑67 (MQ=255)
         gATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGc                                                          <  1:3378562/69‑1 (MQ=255)
         gATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGCCGAAAATAGc                                                          <  1:1252945/69‑1 (MQ=255)
         gATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGCCGAAAATAGc                                                          >  1:2306037/1‑69 (MQ=255)
         gATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGCCGAAAATAGc                                                          >  1:3707727/1‑69 (MQ=255)
         gATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGCCGAAAATAGc                                                          >  1:3431012/1‑69 (MQ=255)
         gATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGCCGAAAATAGc                                                          >  1:174694/1‑69 (MQ=255)
          aTTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGcc                                                         >  1:2284718/1‑69 (MQ=255)
            ttCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCtt                                                       <  1:522267/69‑1 (MQ=255)
             tCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTg                                                      >  1:3332681/1‑69 (MQ=255)
             tCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTg                                                      >  1:1706957/1‑69 (MQ=255)
             tCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTg                                                      >  1:3019713/1‑69 (MQ=255)
                    tATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCtt                                                       <  1:3621829/61‑1 (MQ=255)
                    tATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCtt                                                       <  1:1092638/61‑1 (MQ=255)
                     aTGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCtt                                                       >  1:999765/1‑60 (MQ=255)
                     aTGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCtt                                                       >  1:2205013/1‑60 (MQ=255)
                     aTGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCtt                                                       >  1:3359752/1‑60 (MQ=255)
                     aTGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCtt                                                       <  1:3456451/60‑1 (MQ=255)
                      tGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGaaaa                                                              <  1:3313866/52‑1 (MQ=255)
                      tGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGCCGaaaa                                                              >  1:2830388/1‑52 (MQ=255)
                      tGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGCCGaaaa                                                              >  1:2561734/1‑52 (MQ=255)
                      tGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGCCGaaaa                                                              <  1:3394595/52‑1 (MQ=255)
                      tGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGCCGaaaa                                                              <  1:1317868/52‑1 (MQ=255)
                        caccacGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTc                                                    <  1:2987598/60‑1 (MQ=255)
                        caccacGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGCCGAAAATAGCCTTGTCGCCGCcca                                            <  1:327723/68‑1 (MQ=255)
                           cacGGAACCCGCCGTCGCTTTCATCAAATCATAGTATTTGTCGaa                                                                >  1:3633473/1‑45 (MQ=255)
                           cacGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGaa                                                                >  1:1506430/1‑45 (MQ=255)
                           cacGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGaa                                                                >  1:403331/1‑45 (MQ=255)
                           cacGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGaa                                                                >  1:3293416/1‑45 (MQ=255)
                           cacGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCAt                                        >  1:1005180/1‑69 (MQ=255)
                            acGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGccc                                             >  1:2352592/1‑63 (MQ=255)
                            acGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATc                                       >  1:3198427/1‑69 (MQ=255)
                              ggAACCCGCCATCGCTTTCATCAAATCATAGTATTTGCCGAAAATAGCCtt                                                       <  1:3100892/51‑1 (MQ=255)
                                aaCCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGaaaa                                                              <  1:1394571/42‑1 (MQ=255)
                                 aCCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGaaaa                                                              >  1:1652468/1‑41 (MQ=255)
                                   ccGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCtt                                                       <  1:3472648/46‑1 (MQ=255)
                                       cATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGt                            <  1:736863/69‑1 (MQ=255)
                                       cATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGt                            <  1:3147329/69‑1 (MQ=255)
                                          cGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATc                                       >  1:1605114/1‑55 (MQ=255)
                                          cGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATc                                       >  1:2343924/1‑55 (MQ=255)
                                           gCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACt                                    <  1:1144192/57‑1 (MQ=255)
                                           gCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACt                                    <  1:3313006/57‑1 (MQ=255)
                                           gCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACt                                    <  1:738690/57‑1 (MQ=255)
                                           gCTTTCATCAAATCATAGTATTTGCCGAAAATAGCCTTGTcgc                                                  <  1:665175/43‑1 (MQ=255)
                                           gCTTTCATCAAATCATAGTATTTGCCGAAAATAGCCTTGTCGCCGCccac                                           <  1:3432317/50‑1 (MQ=255)
                                           gCTTTCATCAAATCATAGTATTTGCCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCg                        >  1:3286403/1‑69 (MQ=255)
                                           gCTTTCATCAAATCATAGTATTTGCCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCg                        >  1:408128/1‑69 (MQ=255)
                                           gCTTTCATCAAATCATAGTATTTGCCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCg                        >  1:3034029/1‑69 (MQ=255)
                                           gCTTTCATCAAATCATAGTATTTGCCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCg                        >  1:3453285/1‑69 (MQ=255)
                                             tttCATCAAATCATAGTATTTGCCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTc                      <  1:2043536/69‑1 (MQ=255)
                                              ttCATCAAATCATAGTATTTGCCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCg                     >  1:983261/1‑69 (MQ=255)
                                              ttCATCAAATCATAGTATTTGCCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCg                     >  1:1693492/1‑69 (MQ=255)
                                                 atcaAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATg                  >  1:219941/1‑69 (MQ=255)
                                                   caAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGc                <  1:413406/69‑1 (MQ=255)
                                                   caAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGc                <  1:1083770/69‑1 (MQ=255)
                                                       tCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTc            >  1:2360286/1‑69 (MQ=255)
                                                       tCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTc            >  1:2629732/1‑69 (MQ=255)
                                                            gTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGAt                   <  1:3239831/57‑1 (MQ=255)
                                                            gTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTc       >  1:1516127/1‑69 (MQ=255)
                                                            gTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTc       >  1:2935952/1‑69 (MQ=255)
                                                               tttGCCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTc                      >  1:523680/1‑51 (MQ=255)
                                                               tttGCCGAAAATAGCCTTGTCGC‑‑CCCACCATCACTGACAGCGTACCGTc                      >  1:2553048/1‑49 (MQ=255)
                                                               tttGCCGAAAATAGCCTTGACGCCGCCCACCATCACTGACAGCGTACCGTc                      >  1:2074619/1‑51 (MQ=38)
                                                                ttGTCGAAAATAGCCTTGTCGCCGCCCACCATCACt                                    <  1:1193538/36‑1 (MQ=255)
                                                                ttGTCGAAAATAGCCTTGTCGCCGCCCACCATCACt                                    <  1:1101537/36‑1 (MQ=255)
                                                                ttGTCGAAAATAGCCTTGTCGCCGCCCACCATCACt                                    <  1:3329672/36‑1 (MQ=255)
                                                                 tGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTc                      <  1:330252/49‑1 (MQ=255)
                                                                 tGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTc                      <  1:2111642/49‑1 (MQ=255)
                                                                 tGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTc                      >  1:2792973/1‑49 (MQ=255)
                                                                 tGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTc                      <  1:1604005/49‑1 (MQ=255)
                                                                 tGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCAccgc  >  1:360377/1‑69 (MQ=255)
                                                                  gTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTc       <  1:133328/63‑1 (MQ=255)
                                                                   |                                                                  
ACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGC  >  W3110S.gb/3272109‑3272242

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: