Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,475,701 T→C S352P (TCG→CCG)  yheS → ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,475,7010TC96.3% 283.5 / 1.7 82S352P (TCG→CCG) yheSABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (44/35);  total (46/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

GATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACTT  >  NZ_CP009273/3475617‑3475786
                                                                                    |                                                                                     
gATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACa                                                                                  <  2:2476372/90‑1 (MQ=255)
 aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAt                                                                                 >  1:82147/1‑90 (MQ=255)
 aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAt                                                                                 >  2:893725/1‑90 (MQ=255)
 aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAt                                                                                 <  2:2224764/90‑1 (MQ=255)
 aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAt                                                                                 >  2:1813972/1‑90 (MQ=255)
 aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAt                                                                                 >  1:84340/1‑90 (MQ=255)
 aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAt                                                                                 >  1:361781/1‑90 (MQ=255)
  tCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAtt                                                                                <  2:2960785/90‑1 (MQ=255)
  tCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAtt                                                                                >  2:1905453/1‑90 (MQ=255)
  tCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAtt                                                                                <  1:2719875/90‑1 (MQ=255)
  tCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAtt                                                                                <  1:2506685/90‑1 (MQ=255)
     cATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAAt                                                                             >  2:2909243/1‑90 (MQ=255)
       ttattCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCa                                                                           >  1:2054206/1‑90 (MQ=255)
          ttCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAAc                                                                        <  1:1049152/90‑1 (MQ=255)
          ttCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAAc                                                                        <  1:373215/90‑1 (MQ=255)
            ctcgactcgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTg                                                                      <  1:950209/90‑1 (MQ=255)
              cgactcgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGtt                                                                    >  1:1434294/1‑90 (MQ=255)
              cgactcgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGtt                                                                    >  2:2422302/1‑90 (MQ=255)
              cgactcgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGtt                                                                    >  1:1113779/1‑90 (MQ=255)
               gactcgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTa                                                                   >  1:3014302/1‑90 (MQ=255)
                  tcgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGcc                                                                >  1:2602308/1‑90 (MQ=255)
                   cgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCg                                                               >  1:104524/1‑90 (MQ=255)
                   cgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCaa                                                                          >  2:819705/1‑79 (MQ=255)
                   cgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCaa                                                                          <  1:819705/79‑1 (MQ=255)
                   cgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCg                                                               >  1:271464/1‑90 (MQ=255)
                   cgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCg                                                               >  2:607404/1‑90 (MQ=255)
                     aTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGt                                                             <  1:1959285/90‑1 (MQ=255)
                             gAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTgc                                                     >  2:2581686/1‑90 (MQ=255)
                             gAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTgc                                                     >  1:990567/1‑90 (MQ=255)
                             gAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTgc                                                     >  1:58803/1‑90 (MQ=255)
                             gAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTgc                                                     >  1:482162/1‑90 (MQ=255)
                               aCCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTgcgc                                                   >  1:68404/1‑90 (MQ=255)
                               aCCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTgcgc                                                   >  2:227013/1‑90 (MQ=255)
                                  tGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACtt                                                       >  1:1033208/1‑83 (MQ=255)
                                  tGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACtt                                                       <  2:1033208/83‑1 (MQ=255)
                                   ggTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGt                                               >  1:318083/1‑90 (MQ=255)
                                    gTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTc                                              <  1:838423/90‑1 (MQ=255)
                                        ccGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCg                                          <  2:58803/90‑1 (MQ=255)
                                        ccGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCg                                          >  2:1766293/1‑90 (MQ=255)
                                        ccGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCg                                          <  1:1013124/90‑1 (MQ=255)
                                        ccGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCg                                          <  1:2805956/90‑1 (MQ=255)
                                         cGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCgg                                         <  1:2581686/90‑1 (MQ=255)
                                         cGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCgg                                         >  2:700287/1‑90 (MQ=255)
                                          ggCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGt                                        >  2:2205418/1‑90 (MQ=255)
                                          ggCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGt                                        <  2:2926188/90‑1 (MQ=255)
                                          ggCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGt                                        >  2:579322/1‑90 (MQ=255)
                                              cgcgTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGaaa                                    >  2:2934163/1‑90 (MQ=255)
                                               gcgTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAAt                                   >  2:814996/1‑90 (MQ=255)
                                                    ttGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTc                              <  1:2909243/90‑1 (MQ=255)
                                                    ttGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTc                              <  2:1225779/90‑1 (MQ=255)
                                                    ttGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTc                              <  2:1308573/90‑1 (MQ=255)
                                                        tCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGc                          >  2:676331/1‑90 (MQ=255)
                                                         cTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCg                         >  1:2647218/1‑90 (MQ=255)
                                                          tGTTAGCCGGCAATGGCTCGGGTAAACCGGCATTAATCAAACTGGTACCCCGGGAACTTGCGCCGGTCAGCGGGGAAATTGGTCTGGGGa                        >  1:1352384/1‑90 (MQ=255)
                                                           gTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAAtt                                  >  1:1024653/1‑79 (MQ=255)
                                                           gTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAAtt                                  <  2:1024653/79‑1 (MQ=255)
                                                           gTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGaa                       >  2:874706/1‑90 (MQ=255)
                                                           gTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGaa                       >  1:225012/1‑90 (MQ=255)
                                                           gTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGaa                       >  1:1293366/1‑90 (MQ=255)
                                                           gTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGaa                       >  2:1477993/1‑90 (MQ=255)
                                                            ttAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGaaa                      <  1:700287/90‑1 (MQ=255)
                                                            ttAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGaaa                      <  1:814996/90‑1 (MQ=255)
                                                            ttAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGCTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGaaa                      <  2:509416/90‑1 (MQ=255)
                                                                gCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGa                  >  2:2546854/1‑90 (MQ=255)
                                                                 ccGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGAt                 <  1:489999/90‑1 (MQ=255)
                                                                  cGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATc                <  2:2740141/90‑1 (MQ=255)
                                                                    cAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCaa              <  1:299772/90‑1 (MQ=255)
                                                                      aTGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAAc            <  1:1503080/90‑1 (MQ=255)
                                                                         gcgcGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGt                                        >  2:570074/1‑59 (MQ=255)
                                                                         gcgcGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGt                                        <  1:570074/59‑1 (MQ=255)
                                                                         gcgcGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGaaa                      >  1:2091513/1‑77 (MQ=255)
                                                                         gcgcGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGaaa                      <  2:2091513/77‑1 (MQ=255)
                                                                         gcgcGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCg         >  1:1901047/1‑90 (MQ=255)
                                                                         gcgcGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCg         >  2:123214/1‑90 (MQ=255)
                                                                         gcgcGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCg         >  1:1751650/1‑90 (MQ=255)
                                                                         gcgcGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCg         >  2:2161989/1‑90 (MQ=255)
                                                                           gcgGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGc       <  2:2039901/90‑1 (MQ=255)
                                                                           gcgGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGc       <  1:1718951/90‑1 (MQ=255)
                                                                           gcgGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGc       <  1:227013/90‑1 (MQ=255)
                                                                                tAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAACCTCGGTTACtt  <  1:1633999/90‑1 (MQ=255)
                                                                                tAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACtt  <  2:225012/90‑1 (MQ=255)
                                                                                tAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACtt  <  1:2172592/90‑1 (MQ=255)
                                                                                    |                                                                                     
GATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACTT  >  NZ_CP009273/3475617‑3475786

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: