Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,507,931 C→T 6.7% intergenic (‑356/+119) ald1 ← / ← ypcP Alanine dehydrogenase 1/5'‑3' exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,507,9310CT 6.7% 171.3 / 2.8 60intergenic (‑356/+119)ald1/ypcPAlanine dehydrogenase 1/5'‑3' exonuclease
Reads supporting (aligned to +/- strand):  ref base C (26/30);  new base T (2/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

CTCAAAAACTAATTAATCAAATATGTTTTTTAGATTAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACTTTTTAAAAATGGGATGAAATA  >  USA300TCH1516_ALE/1507791‑1508068
                                                                                                                                            |                                                                                                                                         
cTCAAAAACTAATTAATCAAATATGTTTTTTAGATTAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTTGTTAATAAAATATGTATCTATCAAAAATc                                                                                                                                           >  2:833673/1‑141 (MQ=255)
  cAAAAACTAATTAATCAAATATGTTTTTTAGATTAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCgg                                                                                                                                         <  2:78092/141‑1 (MQ=255)
  cAAAAACTAATTAATCAAATATGTTTTTTAGATTAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCgg                                                                                                                                         <  2:5744/141‑1 (MQ=255)
   aaaaaCTAATTAATCAAATATGTTTTTTAGATTAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGc                                                                                                                                        >  1:749010/1‑141 (MQ=255)
    aaaaCTAATTAATCAAATATGTTTTTTAGTTTAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAGATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTTGGCt                                                                                                                                       >  2:547884/1‑141 (MQ=255)
    aaaaCTAATTAATCAAATATGTTTTTTAGATTAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTTGGCt                                                                                                                                       >  2:158305/1‑141 (MQ=255)
    aaaaCTAATTAATCAAATATGTTTTTTAGATTAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCt                                                                                                                                       >  1:380971/1‑141 (MQ=255)
     aaaCTAATTAATCAAATATGTTTTTTAGATTAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCtt                                                                                                                                      <  1:256300/141‑1 (MQ=255)
     aaaCTAATTAATCAAATATGTTTTTTAGATTAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCtt                                                                                                                                      <  1:296827/141‑1 (MQ=255)
        cTAATTAATCAAATATGTTTTTTATATTAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCttatc                                                                                                                                   <  1:879032/141‑2 (MQ=255)
           attaattaAATATTTTTTTTAGATTAATAAAGTTGTAAATGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCttatttat                                                                                                                                <  2:115412/141‑1 (MQ=255)
           attaatCAAATATGTTTTTTAGATTAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCttatttat                                                                                                                                <  2:179107/141‑1 (MQ=255)
                     tatGTTTTTTAGATTAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGAc                                                                                                                      <  2:185129/141‑1 (MQ=255)
                          tttttAGATTAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATc                                                                                                                 <  2:416810/141‑1 (MQ=255)
                                 aTTAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATAc                                                                                                          <  1:618915/141‑1 (MQ=255)
                                  ttAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACg                                                                                                         <  2:634685/141‑1 (MQ=255)
                                           ttGTAACTGTTCTATTTTGTATTGTATGTTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATCTGTATCTATCAAAAAGCGGCTCATTTATTCACTCGTACATATCCAAATACGCAAACAGAc                                                                                                <  2:201701/141‑1 (MQ=255)
                                           ttGTAACTGTACTATTTTGTAGTGTATGGTAAATTATTTGGATGAAATATAATTCTTTATACTCTATTGATTCATCAGATATGTATCTATCAAAATTCGGCTTATTTTTTCACTCTGACATATCCAAAAAAAAAAAAAGAc                                                                                                <  2:417356/141‑1 (MQ=255)
                                              tAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTAt                                                                                             <  1:708535/141‑1 (MQ=255)
                                              tAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTAt                                                                                             <  2:847174/141‑1 (MQ=255)
                                                   gTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTa                                                                                        >  1:299187/1‑141 (MQ=255)
                                                           ttttAGTGTATGGGAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCCGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTAt                                                                                <  2:709603/138‑1 (MQ=255)
                                                           ttGTGGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTAt                                                                                <  2:259341/141‑1 (MQ=255)
                                                            tGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATc                                                                               <  2:145249/141‑1 (MQ=255)
                                                                gtgtATGGTAATTTATTTGGATGAAATATAATTCTTTATACTGTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTCCCTTGTATCGATa                                                                           >  2:262751/1‑141 (MQ=255)
                                                                gtgtATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATa                                                                           >  1:95345/1‑141 (MQ=255)
                                                                              atttGGGTGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTTTCGATGGGAGTAGGCttttt                                                             >  1:68070/1‑141 (MQ=255)
                                                                               tttGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTta                                                            >  2:895684/1‑141 (MQ=255)
                                                                               tttGAATGAAATATAATTCTTTATACTCTACTGATTAATCAAATATGTATCTAACAAAATTCGGCTTATTTATTCACTCTGACATATCTAAATCCGCACAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTta                                                            >  2:574990/1‑141 (MQ=255)
                                                                                  ggctGAAATATAATTCTTTATGCTTTATTGATTAATCAAATATGTATCTATCAAACTTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACATTGTATCGATAGAAATAGTCTTTTTATAc                                                         <  2:522676/138‑1 (MQ=255)
                                                                                    aTGAAATATAATTCTTTATCCTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTTTCGATAGAAATAGTCTTTTTATACtt                                                       <  1:455809/141‑1 (MQ=255)
                                                                                    aTGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACtt                                                       <  1:833673/141‑1 (MQ=255)
                                                                                       aaaTATAATTTTTTATACTCTATTGATTAATCAAATATGTATCGATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAAt                                                    <  2:574561/141‑1 (MQ=255)
                                                                                             aaTTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTACTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGttt                                              >  1:568010/1‑141 (MQ=255)
                                                                                              aTTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAGATACTATTTCTACCTTGTATCAATAGAAATAGTCTTTTTATACTTAATTGGTTTa                                             >  2:517544/1‑141 (MQ=255)
                                                                                                  ttgctaCTCTATTGATTAATCAAACATGTAGCTATCAAAATTCGGCTTATTTATTCACGCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAAttt                                         <  1:56865/137‑1 (MQ=255)
                                                                                                      taCTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCa                                     <  2:561704/141‑1 (MQ=255)
                                                                                                      taCTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCa                                     <  2:561557/141‑1 (MQ=255)
                                                                                                      taCTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTACCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCa                                     <  2:470996/141‑1 (MQ=255)
                                                                                                         tctATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTACTTGGTTTAATTTTTCATaa                                  >  2:329065/1‑141 (MQ=255)
                                                                                                         tctATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTACTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATaa                                  >  2:600037/1‑141 (MQ=255)
                                                                                                          ctATAGATTAATAAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATGGTCTTTATATACTTCATTGGTTTAATTTTTCATaaa                                 >  2:426230/1‑141 (MQ=255)
                                                                                                                     tCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTATTTTTTCATAAAATCTTATACtt                      >  2:884547/1‑141 (MQ=255)
                                                                                                                     tCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACtt                      >  2:561777/1‑141 (MQ=255)
                                                                                                                     tCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACtt                      >  2:208934/1‑141 (MQ=255)
                                                                                                                      cAAATATGTATCTATCAAAATTTGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACttt                     <  1:158305/141‑1 (MQ=255)
                                                                                                                      cAAATATGTATCTATCAAAATTTGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACttt                     <  1:547884/141‑1 (MQ=255)
                                                                                                                      cAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACttt                     <  2:380971/141‑1 (MQ=255)
                                                                                                                      cAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACttt                     >  1:468976/1‑141 (MQ=255)
                                                                                                                      cAAATATGTAGCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACttt                     <  2:265680/141‑1 (MQ=255)
                                                                                                                        aaTATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACttttt                   >  2:784045/1‑141 (MQ=255)
                                                                                                                        aaTATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACttttt                   <  2:917374/141‑1 (MQ=255)
                                                                                                                             gTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCACCTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACTTTTTaaaaa              >  2:428664/1‑141 (MQ=255)
                                                                                                                              tatctatcAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACTTTTTAAAAAt             >  2:581745/1‑141 (MQ=255)
                                                                                                                              tacctatcAAACTTCGGCTTATTTATTCAGTCTGGCATATCCAAATACGCAATAAGACTTTTTCTGCCTTGTATCAGTAGAAATATTATTTTTTTACTTAATAGGTTTAATTTTTCATAAAATCTTATACTTTTTcaaatt             >  2:271824/4‑139 (MQ=255)
                                                                                                                                tctatcAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACTTTTTAAAAATgg           >  1:811486/1‑141 (MQ=255)
                                                                                                                                 ctatcAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACTTTTTAAAAATggg          >  1:912486/1‑141 (MQ=255)
                                                                                                                                     cAAAATTCGGCTTATTTATTCACTCTGAAATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTAAACAAAAAAAAAATGGGATGa      <  2:697098/141‑1 (MQ=255)
                                                                                                                                         aTTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATATTTTTTAAAAATGGGATGAAatt  >  2:694722/1‑140 (MQ=255)
                                                                                                                                         aTTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACTTTTTAAAAATGGGATGAAata  >  2:5933/1‑141 (MQ=255)
                                                                                                                                         aTTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACTTTTTAAAAATGGGATGAAata  >  2:167664/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                         
CTCAAAAACTAATTAATCAAATATGTTTTTTAGATTAATAAAGTTGTAACTGTACTATTTTGTAGTGTATGGTAATTTATTTGGATGAAATATAATTCTTTATACTCTATTGATTAATCAAATATGTATCTATCAAAATTCGGCTTATTTATTCACTCTGACATATCCAAATACGCAAAAAGACTATTTCTACCTTGTATCGATAGAAATAGTCTTTTTATACTTAATTGGTTTAATTTTTCATAAAATCTTATACTTTTTAAAAATGGGATGAAATA  >  USA300TCH1516_ALE/1507791‑1508068

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: