Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,272,348 C→T 40.5% 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,3480CT40.5% 49.6 / 84.5 74G89S (GGT→AGT) garRtartronate semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base C (20/24);  new base T (14/16);  total (34/40)
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.79e-01

AGCGCTTCGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAGGGGAGTTTGG  >  W3110S.gb/3272284‑3272414
                                                                |                                                                  
aGCGCTTCGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACTTGGCtt                                                              <  1:3831791/71‑1 (MQ=255)
 gcgcTTCGCTGATTTCACGGCTTGCCAGCGGTGCGAAAACTCATATCGATCAATACCGTACCTGGCTTcgc                                                           <  1:3414188/71‑1 (MQ=255)
  cgcTTCGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTcg                                                            >  1:1563862/1‑71 (MQ=255)
   gcTTCGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACTTGGCTTcgc                                                           <  1:955133/71‑1 (MQ=255)
   gcTTCGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACTTGGCTTcgc                                                           <  1:1937233/71‑1 (MQ=255)
   gcTTCGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTcg                                                            >  1:3009037/1‑70 (MQ=255)
       cGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACTTGGCTTCgcgc                                                         >  1:2799022/1‑69 (MQ=255)
       cGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCt                                                       >  1:4235784/1‑71 (MQ=255)
         cTGATTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCt                                                               <  1:1855581/61‑1 (MQ=255)
            aTTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACTTGGc                                                                >  1:3623648/1‑57 (MQ=255)
            aTTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTcgcg                                                          <  1:2025281/63‑1 (MQ=255)
            aTTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCaat                                                  <  1:4286505/71‑1 (MQ=255)
            aTTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCaat                                                  <  1:3846851/71‑1 (MQ=255)
             tttCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCaata                                                 <  1:2084455/71‑1 (MQ=255)
                 aCGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGc                                             <  1:24885/71‑1 (MQ=255)
                  cGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGcc                                            <  1:942345/71‑1 (MQ=255)
                    gCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCAt                                          >  1:4464487/1‑71 (MQ=255)
                       tGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCAtt                                         <  1:4191695/69‑1 (MQ=255)
                       tGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCAtt                                         <  1:2853476/69‑1 (MQ=255)
                       tGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCAtt                                         <  1:4403066/69‑1 (MQ=255)
                        gCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATtct                                       <  1:303/70‑1 (MQ=255)
                         ccAGCGGTGCGATAGAACTCATATCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATtctc                                      <  1:2340183/70‑1 (MQ=255)
                              ggTGCGATAGAACTCATATCGATCAATACCGTACCTGGCtt                                                              >  1:1428847/1‑41 (MQ=255)
                              ggTGCGATAGAACTCATATCGATCAATACCGTACCTGGCtt                                                              >  1:3534902/1‑41 (MQ=255)
                                  cGATAGAACTCATATCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGcc                            <  1:980483/71‑1 (MQ=255)
                                  cGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCa                                           <  1:4228658/56‑1 (MQ=255)
                                  cGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGcc                            <  1:1365391/71‑1 (MQ=255)
                                    aTAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGc                                             <  1:542907/52‑1 (MQ=255)
                                    aTACAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCAtt                                         >  1:4455537/1‑56 (MQ=255)
                                     tAGAACTCATATCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCAcc                         <  1:3164881/71‑1 (MQ=255)
                                     tAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCaata                                                 >  1:590898/1‑47 (MQ=255)
                                     tAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCaaaa                                                 >  1:2971525/1‑45 (MQ=255)
                                     tAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCaa                                                   >  1:3871551/1‑45 (MQ=255)
                                      aGAATTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACct                        >  1:425994/1‑71 (MQ=255)
                                      aGAACTCATATCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACct                        >  1:213590/1‑71 (MQ=255)
                                           tCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCTATAATGCCATTCTCACCCAGCGCCAc                          <  1:2629840/64‑1 (MQ=255)
                                             atatCGATCAATACCGTACTTGGCTTCGCGCCTTCaataa                                                >  1:2307593/1‑40 (MQ=255)
                                             atatCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAgcg                              >  1:3701452/1‑58 (MQ=255)
                                             atatCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTca                  <  1:4469871/70‑1 (MQ=255)
                                              tatCGATCAATACCGTACTTGGCTTCGCGCCTTCa                                                    <  1:766385/35‑1 (MQ=255)
                                              tatCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTca                  <  1:3359703/69‑1 (MQ=255)
                                              tatCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTca                  <  1:2733210/69‑1 (MQ=255)
                                               atCGATCAATACCGTACCTGGCTTCGCGCCTTCaa                                                   >  1:531623/1‑35 (MQ=255)
                                               atCGATCAATACCGTACCTGGCTTCGCGCCTTCaa                                                   <  1:677157/35‑1 (MQ=255)
                                                tCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTcac                 >  1:2763330/1‑68 (MQ=255)
                                                tCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt               >  1:1917008/1‑70 (MQ=255)
                                                tCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt               >  1:4073580/1‑70 (MQ=255)
                                                tCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt               <  1:1231451/70‑1 (MQ=255)
                                                tCGATCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt               <  1:220741/70‑1 (MQ=255)
                                                tCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCAt                                          >  1:180334/1‑43 (MQ=255)
                                                tCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt               <  1:3926279/70‑1 (MQ=255)
                                                tCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt               >  1:327459/1‑70 (MQ=255)
                                                tCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt               >  1:1260960/1‑70 (MQ=255)
                                                tCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt               <  1:4135641/70‑1 (MQ=255)
                                                tCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt               <  1:443558/70‑1 (MQ=255)
                                                tCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATg              <  1:1036198/71‑1 (MQ=255)
                                                tCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATg              <  1:2546328/71‑1 (MQ=255)
                                                 cGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt               >  1:3867814/1‑69 (MQ=255)
                                                  gATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCttt                    <  1:2201046/63‑1 (MQ=255)
                                                   aTCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt               <  1:856454/67‑1 (MQ=255)
                                                   aTCAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACAt               <  1:818781/67‑1 (MQ=255)
                                                   aTCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAgg           <  1:1801794/71‑1 (MQ=255)
                                                     cAATACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCAACCAGCGCAACCTCTTTCACATGAggg          >  1:950629/1‑70 (MQ=255)
                                                     cAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCa                           >  1:4607606/1‑53 (MQ=255)
                                                     aaaTACCGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAggg          >  1:3011130/2‑70 (MQ=255)
                                                        tACCGTACTTGGCTTCGCGCCTTTAATAATGCCATTCTCACCCAGCGCCAc                          >  1:1572644/1‑51 (MQ=255)
                                                          ccGTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCAc                          >  1:2120231/1‑49 (MQ=255)
                                                          ccGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCa                           >  1:2809553/1‑48 (MQ=255)
                                                          ccGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCAc                          >  1:2931494/1‑49 (MQ=255)
                                                          ccGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAGGGGAGttt    <  1:4029368/71‑1 (MQ=255)
                                                            gTACTTGGCTTCGCGCCTTCAATAATGCCATTCTCa                                     >  1:1246207/1‑36 (MQ=255)
                                                            gTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAGGGGAGTTTgg  <  1:4083854/71‑1 (MQ=255)
                                                            gTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAGGGGAGTTTg   >  1:4021327/1‑70 (MQ=255)
                                                             tACTTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAGGGGAGTTTgg  >  1:847428/1‑70 (MQ=255)
                                                                |                                                                  
AGCGCTTCGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAACTCATATCGATCAATACCGTACCTGGCTTCGCGCCTTCAATAATGCCATTCTCACCCAGCGCCACCTCTTTCACATGAGGGGAGTTTGG  >  W3110S.gb/3272284‑3272414

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