Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,628,458 G→A 7.3% G250G (GGC→GGT)  gtaB ← UTP‑‑glucose‑1‑phosphate uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,628,4580GA 7.3% 157.6 / 3.2 55G250G (GGC→GGT) gtaBUTP‑‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (27/24);  new base A (2/2);  total (29/26)
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.70e-01

ACTGAGTATCATGCTTTTAATGTTATAAACCTAATGCTTTAATAAATCGTGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGTTTTTAAATAATCGAAAATTTCAGGCGTTAATACATAACG  >  USA300TCH1516_ALE/2628319‑2628592
                                                                                                                                           |                                                                                                                                      
aCTGAGTATCATGCTTTTAATGTTATAAACCTAATGCTTTAATAAATCGTGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGc                                                                                                                                       <  1:149871/141‑1 (MQ=255)
  tGAGTATCATGCTTTTAATGTTATAAACCTAATGCTTTAATAAATCGTGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCt                                                                                                                                     <  1:84020/141‑1 (MQ=255)
          aggCTTTTAATGGTAGAAACCTAGTGCTGTAATAAAGCGGGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTc                                                                                                                             <  1:152898/139‑1 (MQ=255)
               tttAATGTTATAAACCTAATGCTTTAATAAATCGTGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATg                                                                                                                        <  1:347467/141‑1 (MQ=255)
               tttAATGTTATAAACCAAATGCTTTAATAAATCGTGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATg                                                                                                                        >  1:363360/1‑141 (MQ=255)
                  aaTGTTATAAACCTAATGCTTTAATAAATCGTGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGcct                                                                                                                     >  2:180020/1‑139 (MQ=255)
                   aTGTTATAAACCTAATGCTTTAATAAATCGTGTTAATTCTTCTCGCATACTGTCATCTTTCCATGCATATTCTCTGGTAGTTTTAACGACGCCTAATTTTTCCCAACCGTCATACCGTTCGCCTTCGAAGTCATAGGcgtg                                                                                                                    >  2:17489/1‑138 (MQ=255)
                     tttATAAACATAATGCTTTAATAAATCGTGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATaca                                                                                                                  <  2:235304/140‑1 (MQ=255)
                        ataAACCTAATGCTTTAATAAATCGTGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCACCTTCGAAGTCATATGCATACACtt                                                                                                               >  2:487004/1‑141 (MQ=255)
                        ataAACCTAATGCTTTAATAAATCGTGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCACCTTCGAAGTCATATGCATACACtt                                                                                                               >  2:186026/1‑141 (MQ=255)
                           aaCCTAATGCTTTAATAAATCGTGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGt                                                                                                            <  2:88698/141‑1 (MQ=255)
                             ccTACTGCTTTAATAAATCGTGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTa                                                                                                          <  2:549856/141‑1 (MQ=255)
                                          taaATCGTGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCCTTATTCATAc                                                                                             >  1:209384/1‑141 (MQ=255)
                                           aaaTCGTGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTGGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTAGCCTTCGAAGTCATATGCATACACTTGATTATCATTCTTCATACg                                                                                            >  1:313398/1‑141 (MQ=255)
                                            aaTCGTTTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGt                                                                                           <  1:52441/141‑1 (MQ=255)
                                              tCGTGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTc                                                                                         >  1:363414/1‑141 (MQ=255)
                                              acGTGTTAATTCTTCTTGCATACTGTCATCTTTTAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTc                                                                                         <  2:149885/140‑1 (MQ=255)
                                               cGTGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCa                                                                                        >  1:202733/1‑141 (MQ=255)
                                                 tgtTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAATTCATATGCATACACTTGGTTATCATTATTCATACGTTCAAt                                                                                      <  2:196332/141‑1 (MQ=255)
                                                     aaTTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCa                                                                                  <  2:422532/141‑1 (MQ=255)
                                                         cttctCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTg                                                                              >  1:359179/1‑141 (MQ=255)
                                                            ctcGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTa                                                                           >  2:555477/1‑141 (MQ=255)
                                                               gCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACt                                                                        >  2:139355/1‑141 (MQ=255)
                                                                cATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTg                                                                       >  1:148104/1‑141 (MQ=255)
                                                                      ttCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAAGTTTTCGCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTc                                                                 <  1:346/140‑1 (MQ=255)
                                                                      gTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTCCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTc                                                                 >  2:68624/1‑141 (MQ=255)
                                                                        cttCTTTGCATGCATATTCTATGTTTGGTGTAACGAAGCTTAATTTTTCTCCAACGTCAGAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGt                                                               <  2:275156/139‑1 (MQ=255)
                                                                              tCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCACACGTTCAATCCCATCTGTTAACTGAATTTCGTTACCTg                                                         >  2:112754/1‑141 (MQ=255)
                                                                              tCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAAGTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTCCAATCGCATCTGTTAACTGAATTTGGTTACCTg                                                         >  2:284895/1‑141 (MQ=255)
                                                                                  tGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGcc                                                     <  1:555477/141‑1 (MQ=255)
                                                                                  tGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGcc                                                     >  2:26737/1‑141 (MQ=255)
                                                                                   gCATATTATATGGTATCTTTGTCTAACCCTATTTTTCTTCCACCGTCACTACGTCTGCCTTCGCAGTCATATGCATCTCCTTGCTCCTCCTTATTCCTCCGCTCCCGCGCATCTGTTAACTGAATTTTCTCACCTGCGCct                                                    <  2:14023/141‑1 (MQ=255)
                                                                                      tatTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCcttct                                                 >  1:295153/1‑141 (MQ=255)
                                                                                        ttCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCttt                                               <  2:202733/141‑1 (MQ=255)
                                                                                          cTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGTCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCttttt                                             <  2:400443/141‑1 (MQ=255)
                                                                                              ggTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAACCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGt                                         >  2:162491/1‑141 (MQ=255)
                                                                                               gTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGtt                                        <  1:139355/141‑1 (MQ=255)
                                                                                                tAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCACCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGttt                                       <  1:186026/141‑1 (MQ=255)
                                                                                                tAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCACCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGttt                                       <  1:487004/141‑1 (MQ=255)
                                                                                                 aGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGCAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGtttt                                      <  1:416282/141‑1 (MQ=255)
                                                                                                   tttcAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGTTTTTa                                    <  2:494159/141‑1 (MQ=255)
                                                                                                       aaCGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGTTTTTAaata                                >  1:154655/1‑141 (MQ=255)
                                                                                                       aaCGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGTTTTTAaata                                >  1:23820/1‑141 (MQ=255)
                                                                                                        aCGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTCGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGTTTTTAaataa                               >  1:322781/1‑141 (MQ=255)
                                                                                                              ccTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCTTACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTCCGTTTTTCATTAATCGaac                         >  2:395097/1‑140 (MQ=255)
                                                                                                                 aaTTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTCTTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGTTTTTAAATAATCGAAAAtt                      >  1:274239/1‑141 (MQ=255)
                                                                                                                 aaTTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGTTTTTAAATAATCGAAAAtt                      <  1:526537/141‑1 (MQ=255)
                                                                                                                  aTTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGTTTTTAAATAATCGAAAAttt                     >  2:234023/1‑141 (MQ=255)
                                                                                                                   tttttCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGTTTTTAAATAATCGAAAATTTc                    <  1:530822/141‑1 (MQ=255)
                                                                                                                      ttCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGTTTTTAAATAATCGAACATTTCAgg                 >  1:159068/1‑141 (MQ=255)
                                                                                                                          ccAACGTCATAACGTCCGCCTTCGAAGTCATATGCATACACTTGGTTATCGTTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGTTTTTAAATAATCGAAAATTTCAGGCGtt             <  2:416699/141‑1 (MQ=255)
                                                                                                                             aCGTCATAACGTTCGCCTTCGAAGGCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCCGCGCCTTCTTTTTGCGTTTTTAAATAATCTAAAATTTCAGGCGTTatt          >  1:195514/1‑139 (MQ=255)
                                                                                                                                  aTAACGTTCGCCTTCGAAGTCATCTGCATACACTTGGTTATCATTTTTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTCCGCCTTCTTTTTGCGTTGTTAAATAATCGAAAATTGCGGGGGTTAATACATa     >  1:113038/1‑141 (MQ=255)
                                                                                                                                  aTAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGTTTTTAAATAATCGAAAATTTCAGGCGTTAATACATa     <  1:85739/141‑1 (MQ=255)
                                                                                                                                     aCGTTCGCCTTCGAAGTCATAAGCATACGATTGATTATCATTATTCATCCGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGTTTTTAATTAATCGAAAATTTCAGGCGTTAATACATAACg  >  1:296302/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                      
ACTGAGTATCATGCTTTTAATGTTATAAACCTAATGCTTTAATAAATCGTGTTAATTCTTCTCGCATACTGTCATCTTTCAATGCATATTCTATGGTAGTTTTAACGAAGCCTAATTTTTCTCCAACGTCATAACGTTCGCCTTCGAAGTCATATGCATACACTTGGTTATCATTATTCATACGTTCAATCGCATCTGTTAACTGAATTTCGTTACCTGCGCCTTCTTTTTGCGTTTTTAAATAATCGAAAATTTCAGGCGTTAATACATAACG  >  USA300TCH1516_ALE/2628319‑2628592

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: ‑