Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,761,889 C→A D388Y (GAT→TAT)  fpgS ← Folylpolyglutamate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,761,8890CA96.3% 151.9 / ‑4.3 55D388Y (GAT→TAT) fpgSFolylpolyglutamate synthase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (29/23);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATCTACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGTTTTC  >  USA300TCH1516_ALE/1761754‑1762026
                                                                                                                                       |                                                                                                                                         
cttcATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATatacta                                                                                                                                      >  2:502512/1‑141 (MQ=255)
  tcATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAAt                                                                                                                                    >  2:482164/1‑141 (MQ=255)
    aTGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTg                                                                                                                                  >  2:152977/1‑141 (MQ=255)
      ggCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCTTATACTAATTGTa                                                                                                                                >  1:250558/1‑141 (MQ=255)
       gCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTAt                                                                                                                               <  1:717901/141‑1 (MQ=255)
        cTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTGACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTAtt                                                                                                                              >  1:493511/1‑141 (MQ=255)
         tATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATTTACTAATTGTAttt                                                                                                                             <  1:502512/141‑1 (MQ=255)
         tATTTTAGTTATAGACTTAGAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTAttt                                                                                                                             >  2:894177/1‑141 (MQ=255)
            tttagttttagtTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTtcat                                                                                                                          <  1:262853/141‑1 (MQ=255)
                gttttagttttaTAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCAtt                                                                                                                      >  2:729223/1‑141 (MQ=255)
                 ttttagttttaTAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAACTTTTTATAAATTCAACGTAGTCATCTACTAATTGTATTTCATCAttt                                                                                                                     >  1:661224/1‑141 (MQ=255)
                 ttttagtatcaTAATGCCTCAAAGTCTAATTTTCAATTAACTGCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCAttt                                                                                                                     <  2:265454/141‑1 (MQ=255)
                      gttttaTAATGCTTCAAAGTCCAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAAt                                                                                                                <  1:784935/141‑1 (MQ=255)
                               tGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATCCTAATTGTATTTCATCATTTTCAATActacctac                                                                                                       >  1:245860/1‑141 (MQ=255)
                                         tCTAATTTTGATTTAACTTCACTTATGAAATACAGACTCCAGGGCCTTCCTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTGTTTCATCATTTTCAATAATACCTTCAATTtcttct                                                                                             >  1:404203/1‑141 (MQ=255)
                                          cTAATTTTGCTTTACCTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCACCGTAGTCATCTACTAATTGTATTTCCTCATTTTCAATACTACCTCCAATTTcttctt                                                                                            >  2:332099/1‑141 (MQ=255)
                                          cTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTcttctt                                                                                            >  1:95142/1‑141 (MQ=255)
                                          cTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTcttctt                                                                                            >  1:126836/1‑141 (MQ=255)
                                           tAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGGTAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCttt                                                                                           >  1:193800/1‑141 (MQ=255)
                                           tAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCttt                                                                                           >  1:374806/1‑141 (MQ=255)
                                               tttGATTTAACTTCACTTTTGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGt                                                                                       <  2:492759/141‑1 (MQ=255)
                                               tttGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGGAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGt                                                                                       <  1:580537/141‑1 (MQ=255)
                                               tttGAGTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGt                                                                                       <  2:453631/141‑1 (MQ=255)
                                                       aaCTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCttt                                                                               <  2:742555/141‑1 (MQ=255)
                                                       aaCTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCttt                                                                               <  1:790472/141‑1 (MQ=255)
                                                            cACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGaa                                                                          <  2:536961/141‑1 (MQ=255)
                                                                tATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTGGGAAAATc                                                                      >  1:344403/1‑141 (MQ=255)
                                                                  tGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCaa                                                                    >  1:214601/1‑141 (MQ=255)
                                                                    aaaTACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAAt                                                                  <  2:195446/141‑1 (MQ=255)
                                                                     aaTACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAAtt                                                                 >  1:293016/1‑141 (MQ=255)
                                                                      aTACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTTCAATTTCTTCTTTGCGTAACGCTTTCGAAAAATCAAATTc                                                                >  1:827455/1‑141 (MQ=255)
                                                                          aGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTTCCTACCATTTCTTCTTTGCGTAACGCTTTCGGAAAATCCAATTCAttt                                                            >  2:233712/1‑141 (MQ=255)
                                                                           gagtACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTg                                                           <  2:50704/138‑1 (MQ=255)
                                                                            aCTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGc                                                          >  1:425451/1‑141 (MQ=255)
                                                                              tACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCAt                                                        <  2:769688/141‑1 (MQ=255)
                                                                                 cGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTCCAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATaaa                                                     >  2:442660/1‑141 (MQ=255)
                                                                                                   cACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAAc                                   >  2:315162/1‑141 (MQ=255)
                                                                                                      cTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTa                                <  1:576391/141‑1 (MQ=255)
                                                                                                       ttGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTaa                               >  1:492798/1‑141 (MQ=255)
                                                                                                          aTAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATCAAACGTATGCGCAATTAAACTTAAATg                            >  2:204832/1‑141 (MQ=255)
                                                                                                          aTAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAACCTTAACTg                            >  2:64705/1‑141 (MQ=255)
                                                                                                          aTAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATg                            >  2:364102/1‑141 (MQ=255)
                                                                                                               tttttATAAATTCAACGTAGTCATATACTAATTGTATTTGATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGa                       <  1:169559/141‑1 (MQ=255)
                                                                                                               tttttATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGa                       <  2:425451/141‑1 (MQ=255)
                                                                                                                   tataAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAACTCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATATTTGACCAt                   >  1:419456/1‑141 (MQ=255)
                                                                                                                   tataAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCAt                   >  1:435345/1‑141 (MQ=255)
                                                                                                                   tataAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACAAt                   >  1:744267/1‑141 (MQ=255)
                                                                                                                      aaaTTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCAtctc                <  2:264434/141‑1 (MQ=255)
                                                                                                                              cGTAGTCATATTCTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCg        <  2:284040/141‑1 (MQ=255)
                                                                                                                              cGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCg        >  2:133589/1‑141 (MQ=255)
                                                                                                                               gTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCgg       <  1:68746/141‑1 (MQ=255)
                                                                                                                                tAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGt      <  1:126345/141‑1 (MQ=255)
                                                                                                                                tAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGt      <  1:571398/141‑1 (MQ=255)
                                                                                                                                tAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGt      <  1:894177/141‑1 (MQ=255)
                                                                                                                                  gTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTCGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGttt    >  2:536937/1‑141 (MQ=255)
                                                                                                                                    cATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGTTTTc  <  1:341436/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                         
CTTCATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATCTACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGTTTTC  >  USA300TCH1516_ALE/1761754‑1762026

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑