Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,272,348 C→T 24.4% G89S (GGT→AGT)  garR ← tartronate semialdehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,272,3480CT24.4% 151.7 / 49.9 82G89S (GGT→AGT) garRtartronate semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base C (28/33);  new base T (11/9);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TCAGCGCTTCGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAGGGGAGTTTGGCAG  >  W3110S.gb/3272282‑3272417
                                                                  |                                                                     
tCAGCGCTTCGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGc                                                                   >  1:3562303/1‑71 (MQ=255)
tCAGCGCTTCGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGc                                                                   >  1:3357519/1‑71 (MQ=255)
        tCGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCgcgc                                                            >  1:1247336/1‑70 (MQ=255)
              aTTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACTTGGCTTCGCGCCTTCaat                                                     <  1:432846/71‑1 (MQ=255)
              aTTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACTTGGCTTCGCGCCTTCaat                                                     <  1:2088461/71‑1 (MQ=255)
              aTTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACTTGGCTTCGCGCCTTCaat                                                     <  1:376002/71‑1 (MQ=255)
              aTTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCaat                                                     <  1:2890193/71‑1 (MQ=255)
              aTTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCaat                                                     <  1:3303732/71‑1 (MQ=255)
              aTTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCaa                                                      <  1:1964692/70‑1 (MQ=255)
              aTTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCa                                                       <  1:3482035/69‑1 (MQ=255)
              aTTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCa                                                       <  1:3636314/69‑1 (MQ=255)
                ttCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACGTGGCTTCGCGCCTTCaataa                                                   >  1:1199109/1‑71 (MQ=255)
                ttCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCaataa                                                   >  1:763866/1‑71 (MQ=255)
                ttCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCaataa                                                   >  1:235336/1‑71 (MQ=255)
                ttCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCaataa                                                   >  1:376716/1‑71 (MQ=255)
                ttCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCaataa                                                   >  1:2930214/1‑71 (MQ=255)
                ttCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCaataa                                                   >  1:2073673/1‑71 (MQ=255)
                ttCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCaataa                                                   >  1:3893657/1‑71 (MQ=255)
                   aCGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGc                                                <  1:561087/71‑1 (MQ=255)
                           ccAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATtctc                                         <  1:1913225/70‑1 (MQ=255)
                              gCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCAccc                                     <  1:474016/71‑1 (MQ=255)
                              gCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCAccc                                     <  1:3827197/71‑1 (MQ=255)
                              gCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCAccc                                     <  1:1122007/71‑1 (MQ=255)
                                  tGCGATAGAACTCATATCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAgcg                                 >  1:1002290/1‑71 (MQ=255)
                                  tGCGATAGAACTCATATCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAgcg                                 >  1:2443866/1‑71 (MQ=255)
                                    cGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCaataa                                                   >  1:1118256/1‑51 (MQ=255)
                                    cGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAgcgc                                <  1:577912/70‑1 (MQ=255)
                                      aTAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATTCCATTCTCACCCAGCGCCAc                             <  1:720608/71‑1 (MQ=255)
                                          aaCTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACctct                         <  1:73371/71‑1 (MQ=255)
                                           aCTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCtt                        >  1:3982787/1‑71 (MQ=255)
                                           aCTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCtt                        >  1:2303348/1‑71 (MQ=255)
                                            cTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGc                                                <  1:2410789/46‑1 (MQ=255)
                                            cTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGc                                                <  1:168038/46‑1 (MQ=255)
                                            cTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGc                                                <  1:1439091/46‑1 (MQ=255)
                                            cTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGc                                                <  1:3104704/46‑1 (MQ=255)
                                              cATATCGATCAATACCGTACTTGGCTTCCCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTca                     >  1:3813997/1‑71 (MQ=255)
                                               atatCGATCAATACCGTACTTGGCTTCGCGCCTTCaataa                                                   >  1:3631748/1‑40 (MQ=255)
                                               atatCGATCAATACCGTACTTGGCTTCGCGCCTTCaataa                                                   >  1:661474/1‑40 (MQ=255)
                                               atatCGATCAATACCGTACTTGGCTTCGCGCCTTCaataa                                                   >  1:3021207/1‑40 (MQ=255)
                                               atatCGATCAATACCGTACTTGGCTTCGCGCCTTCaataa                                                   >  1:1102235/1‑40 (MQ=255)
                                                tatCGATCAATACCGTACTTGGCTTCGCGCCTTCaat                                                     <  1:733382/37‑1 (MQ=255)
                                                  tCGATCAATACCGTACTTGGCTTCGCGCCTTCaataa                                                   >  1:1383160/1‑37 (MQ=255)
                                                  tCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCAGCTCTTTcaca                   >  1:805876/1‑69 (MQ=255)
                                                  tCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt                  >  1:1218200/1‑70 (MQ=255)
                                                  tCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATg                 <  1:3264429/71‑1 (MQ=255)
                                                  tCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATg                 <  1:1172967/71‑1 (MQ=255)
                                                  tCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATg                 <  1:3941760/71‑1 (MQ=255)
                                                  tCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATg                 <  1:3757067/71‑1 (MQ=255)
                                                  tCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATg                 <  1:3575705/71‑1 (MQ=255)
                                                  tCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCa                                              >  1:3684322/1‑42 (MQ=255)
                                                  tCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCa                                              >  1:3807268/1‑42 (MQ=255)
                                                  tCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCa                                              >  1:1315065/1‑42 (MQ=255)
                                                  tCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCa                                              >  1:2655861/1‑42 (MQ=255)
                                                  tCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACctct                         <  1:3917165/63‑1 (MQ=255)
                                                  tCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATg                 <  1:2599987/71‑1 (MQ=255)
                                                  tCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATg                 <  1:447664/71‑1 (MQ=255)
                                                    gATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAg               >  1:3349657/1‑71 (MQ=255)
                                                    gATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATtct                                          >  1:3088252/1‑44 (MQ=255)
                                                          tACCGTACCTGGCTTCGCGCCTTCAATAATGCCAtt                                            <  1:2758919/36‑1 (MQ=255)
                                                          tACCGTACCTGGCTTCGCGCCTTCAATAATGCCAtt                                            <  1:1439841/36‑1 (MQ=255)
                                                          tACCGTACCTGGCTTCGCGCCTTCAATAATGCCAtt                                            <  1:2977367/36‑1 (MQ=255)
                                                          tACCGTACCTGGCTTCGCGCCTTCAATAATGCCATtct                                          <  1:403137/38‑1 (MQ=255)
                                                          tACCGTACCTGGCTTCGCGCCTTCAATAATGCCATtct                                          <  1:2013730/38‑1 (MQ=255)
                                                          tACCGTACCTGGCTTCGCGCCTTCAATAATGCCATtct                                          <  1:404904/38‑1 (MQ=255)
                                                          tACCGTACCTGGCTTCGCGCCTTCAATAATGCCATtct                                          <  1:2443311/38‑1 (MQ=255)
                                                          tACCGTACCTGGCTTCGCGCCTTCAATAATGCCATtct                                          <  1:1519936/38‑1 (MQ=255)
                                                          tACCGTACCTGGCTTCGCGCCTTCAATAATGCCATtct                                          <  1:731368/38‑1 (MQ=255)
                                                          tACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGcc                               <  1:259165/49‑1 (MQ=255)
                                                          tACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTcac                    <  1:320606/60‑1 (MQ=255)
                                                          tACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAGGGGAg          <  1:1220429/70‑1 (MQ=255)
                                                           aCCGTACCTGGCTTCGCGCCTTCAATAATGCCAtt                                            <  1:629590/35‑1 (MQ=255)
                                                            ccGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTcaca                   >  1:1919003/1‑59 (MQ=255)
                                                            ccGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTcaca                   >  1:1850551/1‑59 (MQ=255)
                                                            ccGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt                  >  1:1533060/1‑60 (MQ=255)
                                                            ccGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt                  >  1:2477590/1‑60 (MQ=255)
                                                            ccGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt                  >  1:705086/1‑60 (MQ=255)
                                                            ccGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt                  >  1:3725344/1‑60 (MQ=255)
                                                            ccGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt                  >  1:3609427/1‑60 (MQ=255)
                                                            ccGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt                  >  1:3318776/1‑60 (MQ=255)
                                                            ccGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAGGGGAGttt       >  1:204868/1‑71 (MQ=255)
                                                            ccGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAGGGGAGttt       >  1:2075679/1‑71 (MQ=255)
                                                            ccGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAGGGGAGttt       >  1:3666818/1‑71 (MQ=255)
                                                                 ccTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAGGGGAGTTTGgcag  >  1:586009/1‑71 (MQ=255)
                                                                  cTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt                  <  1:1538828/54‑1 (MQ=255)
                                                                  cTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAGGGGAGTTTGgcag  <  1:3653552/70‑1 (MQ=255)
                                                                  |                                                                     
TCAGCGCTTCGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAGGGGAGTTTGGCAG  >  W3110S.gb/3272282‑3272417

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: