Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,744,473 C→T 82.9% A67T (GCT→ACT)  apt ← Adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,744,4730CT82.9% 113.8 / 20.0 58A67T (GCT→ACT) aptAdenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (4/6);  new base T (27/21);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-01

ACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTG  >  USA300TCH1516_ALE/1744333‑1744612
                                                                                                                                            |                                                                                                                                           
acGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTACAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGc                                                                                                                                             <  1:397982/141‑1 (MQ=255)
acGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGc                                                                                                                                             <  2:195125/141‑1 (MQ=255)
    tGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATTACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACa                                                                                                                                         >  1:189792/1‑141 (MQ=255)
       ccTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACTATATTCTAGGTCATCCTCATAGCGAATGACTTAACGAGGTTATTTTCCTTCTTTTCTAACAGGTGGAAAGGCATTATCCATAGAATTAGTTACAGGa                                                                                                                                      >  2:846/1‑141 (MQ=255)
         tGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACa                                                                                                                                    >  1:522033/1‑141 (MQ=255)
            tttAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGcc                                                                                                                                 >  2:85268/1‑141 (MQ=255)
            tttAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGcc                                                                                                                                 >  2:103678/1‑141 (MQ=255)
               aatgtcATCTTTGTGCATTGTTAAAACATTTGTACCAGATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAAt                                                                                                                              <  2:548808/136‑1 (MQ=255)
                        tttGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGcc                                                                                                                     <  2:494620/141‑1 (MQ=255)
                          tgtgCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCAc                                                                                                                   >  1:309060/1‑141 (MQ=255)
                           gtgCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCAcg                                                                                                                  >  1:160001/1‑141 (MQ=255)
                              cATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACgcgc                                                                                                               >  2:328667/1‑141 (MQ=255)
                                    tAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCCCTGTACACTGAATAAATTACAGGACAGCCAATGATAAAGCCACGCGCTTCAgg                                                                                                         >  1:441912/1‑141 (MQ=255)
                                                 ccATATTCCAGGTCATCCTCATAACGCATGACTTCACGACGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATACAGCCACGCGCTTCAGGTCCTACAACGata                                                                                            <  2:558343/141‑1 (MQ=255)
                                                   atatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                                                                                          <  1:226772/141‑1 (MQ=255)
                                                   atatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCCCTTCAGGTCCTACAACGATATc                                                                                          <  1:297288/141‑1 (MQ=255)
                                                      ttCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAAc                                                                                       <  1:285697/141‑1 (MQ=255)
                                                       tCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCGAAGCCAATCCCCATTGAATAAGTTAGAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCaatc                                                                                      >  1:438703/1‑139 (MQ=255)
                                                        cTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATGTCAACac                                                                                     >  2:201913/1‑140 (MQ=255)
                                                        cTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAt                                                                                     >  1:123589/1‑141 (MQ=255)
                                                        cTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATACCAACac                                                                                     >  1:244174/1‑140 (MQ=255)
                                                          tggTCATACTCATAACGGATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACGGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtct                                                                                   <  2:249173/140‑1 (MQ=255)
                                                           ggTCATCCTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCCAACAGGTGCAAAGCCAACCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtctc                                                                                  <  2:45766/141‑1 (MQ=255)
                                                           ggTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtctc                                                                                  <  1:187224/141‑1 (MQ=255)
                                                             tCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGAAAAGCCAGTCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTg                                                                                >  2:310210/1‑141 (MQ=255)
                                                                  cTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCttt                                                                           <  2:557830/141‑1 (MQ=255)
                                                                       aaCGAATTACTTCACGACGTAATTTCCCTTCTTTACTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCCTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGt                                                                      <  2:353104/141‑1 (MQ=255)
                                                                            aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCCTTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACCACTCCATCTTTTGCGTATTcc                                                                 >  2:155929/1‑140 (MQ=255)
                                                                            aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGCCTTTTGcaaatttg                                                                 >  2:454510/1‑134 (MQ=255)
                                                                              gccgtCACGAGGTAATTGCCCTGCGTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTGCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTAc                                                               <  2:17163/137‑1 (MQ=255)
                                                                              gACGTCACGAGGTAATTCCCCTTTTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTAc                                                               <  2:56469/141‑1 (MQ=255)
                                                                               aCTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACACCGGTCTCAACATCTCTATCTTTTGCGTATTCTACa                                                              >  1:137328/1‑141 (MQ=255)
                                                                                 ttCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATGGAATAAGTTACAGGACAGCCAATGGTAAAGCCACGGGCTGGAGGCGCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAt                                                            >  2:303830/1‑141 (MQ=255)
                                                                                  tCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAtt                                                           <  2:307404/141‑1 (MQ=255)
                                                                                         gTAATTTCCCTTCCTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATACAGCCACGCGCTTCAGGTCATACAACGCTATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTg                                                    <  2:448992/141‑1 (MQ=255)
                                                                                          tAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCTTATTCTACAATATTATCTGt                                                   >  1:38103/1‑141 (MQ=255)
                                                                                           aaTTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGtt                                                  <  1:103678/141‑1 (MQ=255)
                                                                                           aaTTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGtt                                                  <  1:85268/141‑1 (MQ=255)
                                                                                            aTTTCCCTTCTTGTCTAACAGGGGCAAAGCCAAGCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTg                                                 >  2:604087/1‑141 (MQ=255)
                                                                                               tCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTGCAGGTCCTACAACGATATCACCATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCAt                                              >  1:412112/1‑141 (MQ=255)
                                                                                                cccTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATa                                             >  1:481917/1‑141 (MQ=255)
                                                                                                   ttcttTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGcc                                          <  1:321310/141‑1 (MQ=255)
                                                                                                    tcttTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATCGCCa                                         <  2:160001/141‑1 (MQ=255)
                                                                                                     cttTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGCTAAAGCCCCGCGCTTCAGGCTCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCGGTTGCATACCCaa                                        >  1:258345/1‑140 (MQ=255)
                                                                                                      ttttCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCata                                       <  1:89948/141‑1 (MQ=255)
                                                                                                             aCAGGTGCAAAGACAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCa                                <  2:441912/141‑1 (MQ=255)
                                                                                                                  tGCACAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCAtt                           >  2:406978/1‑141 (MQ=255)
                                                                                                                    cAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCCTATGCTTCCCCATTAt                         >  1:26988/1‑141 (MQ=255)
                                                                                                                    cAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTAt                         >  1:256753/1‑141 (MQ=255)
                                                                                                                      atgCCAACCCCCAGTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCGCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATcc                       <  2:123589/139‑1 (MQ=255)
                                                                                                                      aaGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTTTTGCATAGCCATATGCTTCACCATTATcc                       >  1:206563/1‑141 (MQ=255)
                                                                                                                      aaGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATcc                       <  1:587434/141‑1 (MQ=255)
                                                                                                                          cAATCCCCATTGGATAAGTTATAGGACCGCCAATGATATAGCCACGCGCTTCAGGTCCTACAACGATATCAACACCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATaa                   <  2:308021/141‑1 (MQ=255)
                                                                                                                            aTCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATCGAATTATTCTTCACCATTATCCATAAtt                 <  2:131065/141‑1 (MQ=255)
                                                                                                                              ccccATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTgt               >  1:435537/1‑141 (MQ=255)
                                                                                                                                    tGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTttagtact         >  1:148513/1‑139 (MQ=255)
                                                                                                                                     gAAAAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAata        <  1:258937/141‑1 (MQ=255)
                                                                                                                                        tAAGTTACAGGACAACCAAGGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATcc     <  1:310210/141‑1 (MQ=255)
                                                                                                                                           gCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATTCTTCACCATTATCCGTAATTGTAGTAATATCCTTg  >  1:560191/1‑141 (MQ=255)
                                                                                                                                           gCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTg  >  2:589309/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                           
ACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTG  >  USA300TCH1516_ALE/1744333‑1744612

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