Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 118,976 G→A 12.2% G82S (GGT→AGT)  ampD → N‑acetyl‑anhydromuranmyl‑L‑alanine amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb118,9760GA12.2% 135.7 / 9.6 49G82S (GGT→AGT) ampDN‑acetyl‑anhydromuranmyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base G (17/26);  new base A (3/3);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

GCACATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCA  >  W3110S.gb/118910‑119040
                                                                  |                                                                
gCACATCCTTTCTTTTCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGa                                                              <  1:385240/71‑1 (MQ=255)
gCACATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGa                                                              <  1:701949/71‑1 (MQ=255)
   caTCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGaaa                                                            <  1:2627166/70‑1 (MQ=255)
    aTCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTTaaa                                                            <  1:2694288/69‑1 (MQ=255)
       ctttctttGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTc                                                       >  1:1165722/1‑71 (MQ=255)
          tctttGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAg                                                    >  1:812388/1‑71 (MQ=255)
           ctttGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAg                                                    <  1:1050495/70‑1 (MQ=255)
              tGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGt                                                        >  1:812827/1‑63 (MQ=255)
              tGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTa                                                  <  1:943246/69‑1 (MQ=255)
              tGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTa                                                  <  1:917947/69‑1 (MQ=255)
                cTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGt                                               <  1:1442689/70‑1 (MQ=255)
                   agaTCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTACt                                           <  1:2567991/71‑1 (MQ=255)
                      tCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTc                                        <  1:2106670/71‑1 (MQ=255)
                       cGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATa                                                          >  1:2143086/1‑52 (MQ=255)
                             tttGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATAGTGAAATAGTCCAGt                                                   <  1:1499535/53‑1 (MQ=255)
                              ttGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACg                                >  1:681074/1‑71 (MQ=255)
                              ttGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATAGTGAAATAGTCCAGTATGTTCCTTTCGATAAACg                                >  1:430288/1‑71 (MQ=255)
                               tGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTc                                             >  1:1758072/1‑57 (MQ=255)
                               tGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCt                                           >  1:2722102/1‑59 (MQ=255)
                               tGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGt                               >  1:1693611/1‑71 (MQ=255)
                               tGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGt                               >  1:1405592/1‑71 (MQ=255)
                               tGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGt                               >  1:2420721/1‑71 (MQ=255)
                               tGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGt                               >  1:2371070/1‑71 (MQ=255)
                               tGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGt                               >  1:1812629/1‑71 (MQ=255)
                                 cgcgTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTg                              <  1:697516/70‑1 (MQ=255)
                                  gcgTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTgca                            <  1:675972/71‑1 (MQ=255)
                                     tctcCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGa                                      <  1:1193232/58‑1 (MQ=255)
                                     tctcCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTgcatgg                         <  1:732704/71‑1 (MQ=255)
                                       tcCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTg                              >  1:1200763/1‑64 (MQ=255)
                                       tcCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCAAGTATGTTCCTTTCg                                       >  1:786470/1‑55 (MQ=255)
                                       tcCGCTCACTGTTTGATTCGCCGTGATAGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTgcatggca                       >  1:2667650/1‑71 (MQ=255)
                                       tcCGCTCACTGTTTGATTCGCCGTGATAGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTgcatggca                       >  1:66227/1‑71 (MQ=255)
                                         cGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGt                                               <  1:1010887/45‑1 (MQ=255)
                                         cGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTgca                            >  1:1324885/1‑64 (MQ=255)
                                               cTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTg                              <  1:2811936/56‑1 (MQ=255)
                                                tGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCtt                                          <  1:2739056/43‑1 (MQ=255)
                                                tGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCtt                                          <  1:1501463/43‑1 (MQ=255)
                                                tGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCtt                                          <  1:867052/43‑1 (MQ=255)
                                                    tGATTCGCCGTGATAGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGtctc           <  1:739790/70‑1 (MQ=255)
                                                     gATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTgcatggca                       >  1:1176800/1‑57 (MQ=255)
                                                     gATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTgca                            <  1:1471146/52‑1 (MQ=255)
                                                     gATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGtctct          <  1:1182741/70‑1 (MQ=255)
                                                     gATTCGCCGTGATAGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTgcatgg                         <  1:2547492/55‑1 (MQ=255)
                                                      aTTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTgcatggc                        <  1:1850183/55‑1 (MQ=255)
                                                        tCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTgcatggcat                      >  1:1404847/1‑55 (MQ=255)
                                                           ccGTTATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTgcatgg                         <  1:1166272/49‑1 (MQ=255)
                                                            cGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCa  <  1:1829094/71‑1 (MQ=255)
                                                                aTGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTgcatg                          <  1:2602802/43‑1 (MQ=255)
                                                                aTGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTgcatg                          <  1:755462/43‑1 (MQ=255)
                                                                  |                                                                
GCACATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCA  >  W3110S.gb/118910‑119040

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: