Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,373,809 A→T 15.4% intergenic (‑34/+79) pbuG ← / ← topB Guanine/hypoxanthine permease PbuG/DNA topoisomerase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,373,8090AT15.4% 137.3 / 13.2 62intergenic (‑34/+79)pbuG/topBGuanine/hypoxanthine permease PbuG/DNA topoisomerase 3
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base T (7/3);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCTTTAAACGGATTATTGTCTAAACCTTCA  >  USA300TCH1516_ALE/2373674‑2373942
                                                                                                                                       |                                                                                                                                     
aaCTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGttt                                                                                                                                   <  1:160234/140‑1 (MQ=255)
  cggCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTCTCGCACTGCAAATAGTTTTTCACGAATGGATGCTCCTTTTTGTAAAATACCTTAATAGTTTaa                                                                                                                                 <  1:100140/138‑1 (MQ=255)
    gCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAata                                                                                                                               <  1:589589/140‑1 (MQ=255)
      tAAAATATATGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATc                                                                                                                             <  2:233639/140‑1 (MQ=255)
          atatAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTc                                                                                                                         >  1:35124/1‑140 (MQ=255)
          atatAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTc                                                                                                                         >  1:367697/1‑140 (MQ=255)
              aGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTtaa                                                                                                                     >  1:571178/1‑140 (MQ=255)
              aGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGATCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTtaa                                                                                                                     >  1:666867/1‑140 (MQ=255)
                  cATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATaaa                                                                                                                 >  1:228412/1‑140 (MQ=255)
                  cATAGATAAGAAAGTTGTGATACAGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACATTAATAGTTTAATATCATTCTtaataca                                                                                                                 >  2:392450/1‑138 (MQ=255)
                       ataataAAGTTGTGATACCGCCTAAGATTTCTATTTTAAAGTTTGTTTCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACa                                                                                                            <  2:472976/140‑1 (MQ=255)
                        tAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATATCTTAATAGTTTAATATCATTCTTAATaaaataaaa                                                                                                           >  2:338860/1‑135 (MQ=255)
                          aGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTACTTGTTTAATATCTTTCTTCATAATAAACAATa                                                                                                         >  2:344101/1‑140 (MQ=255)
                            aaaGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAACATACCTGAATTTTTTACTATAATTCTGATTAATAAACAGTaac                                                                                                       >  2:244011/1‑138 (MQ=255)
                               gTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCaaa                                                                                                    >  2:434505/1‑140 (MQ=255)
                               gTTGTAATACCGACTAAGATTTCTATTTTATAGTTTGTTCCATATATATCGAAATGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCaaa                                                                                                    <  2:392371/140‑1 (MQ=255)
                                ttGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCaaaa                                                                                                   >  2:619200/1‑140 (MQ=255)
                                    gATACAGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACTAAAAGCGa                                                                                               >  2:580314/1‑140 (MQ=255)
                                         caccTAATATTTCGATTTTAAAGTTTGTTCCATATTTATCCAAATGGAAAATAACAATCACGAATGGCTGCTCCGCGATTTAAAATACCTTAATAGTTTAATATAATTCTTAATAAAAAACAATACCAAAAGCGAACttt                                                                                          <  2:315321/138‑1 (MQ=255)
                                           ccTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTTAAATACCTTAATAGTTTAATATCATTCTTAATAAAAACCAATACCAAAAGCGTCCTTTTc                                                                                        >  1:589258/1‑140 (MQ=255)
                                           ccTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTATAATACCTTAATTGTTTAATTTCATTCTTTATTAAAATCAATACCAAAAGCATTCTTTTc                                                                                        >  2:410627/1‑140 (MQ=255)
                                            cggaGAATTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGGTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCGTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCg                                                                                       <  2:593930/137‑1 (MQ=255)
                                             tAAGATTTCTCTTTTAAAGTTTGTGCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGt                                                                                      <  1:37963/140‑1 (MQ=255)
                                               aGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATTGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGttt                                                                                    <  1:675147/140‑1 (MQ=255)
                                               aGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGttt                                                                                    <  1:434505/140‑1 (MQ=255)
                                                    tctcTTTTAAAGTTTGTTACATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATTGTTTAATATCATTCTTAATAAAAAACAATACAAAAAGAGAACTTTTTGTATATCTg                                                                               >  2:366404/1‑140 (MQ=255)
                                                        ttttAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACTAAAAGCGAACTTTTCGTTTTTTTGATTc                                                                           >  2:612515/1‑140 (MQ=255)
                                                         tttAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCa                                                                          >  1:600316/1‑140 (MQ=255)
                                                            aaaGTGTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATTGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACt                                                                       <  2:589525/140‑1 (MQ=255)
                                                                 ttGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACTTTAATTGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGtt                                                                  >  2:588408/1‑140 (MQ=255)
                                                                 ttGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATTGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCTAACTTTTCTTTTTTCTGATTCAACtttttt                                                                  >  2:119526/1‑140 (MQ=255)
                                                                      ccATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGAtt                                                             >  2:2178/1‑140 (MQ=255)
                                                                      ccATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGAtt                                                             >  1:210779/1‑140 (MQ=255)
                                                                             tATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCta                                                      >  1:253496/1‑140 (MQ=255)
                                                                                    cTGGAAATAGCTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTc                                               <  2:94033/140‑1 (MQ=255)
                                                                                     tGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCa                                              >  1:193439/1‑140 (MQ=255)
                                                                                     tGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCa                                              >  1:499657/1‑140 (MQ=255)
                                                                                      ggAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTTTTTTTTTTTAATAGTTTAATATCATTCTTAATTATTTACAATACCAAAAGCGTACTTTTCGGTTTTCTGATTCATCTTTGTTCGTTTTTTTCTTTAATTTTaa                                             >  2:31509/1‑140 (MQ=255)
                                                                                      ggAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATTCCTTAATATTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGCTTCAACTTTGTTCGATTTTATCTATAAATTCaa                                             >  1:187458/1‑140 (MQ=255)
                                                                                      ggAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTCTCTATAAATTCaa                                             >  2:453441/1‑140 (MQ=255)
                                                                                      ggAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCaa                                             >  1:231508/1‑140 (MQ=255)
                                                                                       gAAATATGTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTAAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAg                                            <  2:52714/140‑1 (MQ=255)
                                                                                       gAAATAGTTTTTGACGCTTGGATGCTACGTTTTTTAAAATAACTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAg                                            <  2:144278/140‑1 (MQ=255)
                                                                                       gAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAg                                            <  1:385886/140‑1 (MQ=255)
                                                                                       gAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAg                                            >  1:704171/1‑140 (MQ=255)
                                                                                       aaaaTAATTTTTCACGAAGAGATGATCATTTTTTTAAAATAACTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAg                                            <  2:338711/139‑1 (MQ=255)
                                                                                        aaaTAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAATTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATTTATAATTTCaatt                                           >  2:329976/1‑138 (MQ=255)
                                                                                          aTAGTTTTTCACGAATGGATGCTCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAATTTCAAGTTCt                                        >  2:324764/1‑140 (MQ=255)
                                                                                           tAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCt                                        <  1:2178/140‑1 (MQ=255)
                                                                                               ttttCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTaa                                    <  1:303350/140‑1 (MQ=255)
                                                                                                     cGAATGGATGCTCCTTTTTTTAAAATACTTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGAGAACTTTTCGTTTTTCTTATTCAACTTTGTTCCATTCTATCTATAAACTCAAGTTCTTTTAtctatt                              >  2:457009/1‑135 (MQ=255)
                                                                                                       aaTGGATGCTCCTTTTTTTAAAATACTTTATTAGTTTAATATCATTTTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTACCTGCGTTGTTATCTATATATTCAAGTTCTTTAATGTATCtt                            >  2:311990/1‑140 (MQ=255)
                                                                                                        aTGGATGCTCCTTTTTTTAAAATACCTTAATTGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCttt                           <  1:561651/140‑1 (MQ=255)
                                                                                                           gATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCTTTaaa                        <  2:571178/140‑1 (MQ=255)
                                                                                                           gATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCTTTaaa                        <  2:666867/140‑1 (MQ=255)
                                                                                                           gATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCTTTaaa                        <  2:704171/140‑1 (MQ=255)
                                                                                                               cTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCTTTAAACGGa                    <  1:619200/140‑1 (MQ=255)
                                                                                                                      ttttAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCTTTAAACGGATTATTGt             <  2:69644/140‑1 (MQ=255)
                                                                                                                            ttttttttAATAGTTTAATATCATTCTTAATACATTTCAATACCATAAGCGTACTTTTCGTTTTTCTGAATCCACTTTTTTTTTTTTTATTTATAACTTCAACTTCTTTAATTCTTCTTTAAACGGATTATTTTCTAAAc       >  2:568235/7‑140 (MQ=255)
                                                                                                                             aTACCTTAATTGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTTTTTAATGTATCTTTAAATGGATTATTGTCTAAAcc      >  2:141795/1‑140 (MQ=255)
                                                                                                                                ccTTGAGAGTTTAATAGCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTGTCTGATTCAACTTTGTTCGATTTTATCTATCAATTCAAGTTCTTTAAGGCATCTCTAAACGGATTATTGTCTAAACCTTc   <  2:13424/140‑1 (MQ=255)
                                                                                                                                 taataTAGTTTAATACCAGTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCTTTAAACGGATTATTGTCTAAACCTTCa  <  2:517448/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
AACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCTTTAAACGGATTATTGTCTAAACCTTCA  >  USA300TCH1516_ALE/2373674‑2373942

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