Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,485,771 A→C 8.4% intergenic (‑89/+292) mqo1 ← / ← lctP_2 putative malate:quinone oxidoreductase 1/L‑lactate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,485,7710AC 8.4% 161.8 / 5.3 64intergenic (‑89/+292)mqo1/lctP_2putative malate:quinone oxidoreductase 1/L‑lactate permease
Reads supporting (aligned to +/- strand):  ref base A (25/33);  new base C (3/3);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

ATACCGCCACCAATTAAGATGACATCTGTTTTGCTATGTTGTGTTGTCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTAATTATATACTTTATGTGCGATAGCACAAAGACAAAACCACAAACCTTAACAGCTATCAACCAATT  >  USA300TCH1516_ALE/2485633‑2485909
                                                                                                                                           |                                                                                                                                           
ataCCGCCACCAATTAAGATGACATCTGTTTTGCTATGTTGTGTTGTCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCaa                                                                                                                                            <  1:147770/140‑1 (MQ=255)
      ccaccaATTAAGATGACATCTGTTTTGCTATGTTGTGTTGTCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCCATAGGtt                                                                                                                                      >  1:573247/1‑140 (MQ=255)
      ccaccaATTAAGATGACATCTGTTTTGCTATGTTGTGTTGTCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCCATAGGtt                                                                                                                                      >  1:78917/1‑140 (MQ=255)
        accaATTAAGATGACATCTGTTTTGCTATGTTGTGTTGTCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTcc                                                                                                                                    <  1:625603/140‑1 (MQ=255)
            aTTAAGTTGACACCTGTTTTGCTTTATTGTTCTGTCATAACAACTACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATATGTTCCCATg                                                                                                                                <  2:603869/140‑1 (MQ=255)
            aTTAAGATGACATCTGTTTTGCTATGTTGTGTTGTCATAACGAATACAGTCACCATATTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATAACACTTTAAATCAATAAGTGTATCAATAGGTTCCAATg                                                                                                                                >  1:62995/1‑140 (MQ=255)
             ttAAGATGACATCTGTTTTGCTATGTTGTGTTGTCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATgg                                                                                                                               <  1:215004/140‑1 (MQ=255)
             ttAAGATGACATCTGTTTTGCTATGTTGTGTTGTCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATgg                                                                                                                               <  1:472898/140‑1 (MQ=255)
              tAAGATGACATCTGTTTTGCTATGTTGTGTTGTCATAACTAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGt                                                                                                                              <  1:620179/140‑1 (MQ=255)
              tAAGATGACATCTGTTTTGCTATGTTGTGTTGTCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGt                                                                                                                              <  1:362097/140‑1 (MQ=255)
              tAAGATGACATCTGTTTTGCTATGTTGTGATGTCATAACAAATACAGTCTCACTTTTTAAAATTAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAAGAGGTTCCCATGGt                                                                                                                              <  2:435407/140‑1 (MQ=255)
               aaGATGACATCTGTTTTGCTATGTTGTGTTGTCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTg                                                                                                                             <  2:545988/140‑1 (MQ=255)
                 gATGACATCTGTTTTGCTATGTTGTGTTGTCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGAt                                                                                                                          <  2:455488/140‑1 (MQ=255)
                   tGACATCTGTTTTGCTATGTTGTGTTGTCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGc                                                                                                                        <  1:5142/140‑1 (MQ=255)
                       aTCTGTTTTGTTATGTTGTCTTGTCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTGTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTa                                                                                                                    <  1:314502/140‑1 (MQ=255)
                        tCTGTTTTGCTATGTTGTGTTGTCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTAt                                                                                                                   <  2:531043/140‑1 (MQ=255)
                         cTGTTTTGCTATGTTGTGTTGTCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAATTTAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTAtt                                                                                                                  <  1:352125/140‑1 (MQ=255)
                            ttttGCTATGTTGTGTTGTCATAACAACTACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAACTCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACt                                                                                                               <  1:377524/140‑1 (MQ=255)
                              ttGCTATGTTGTGTTGTCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTAt                                                                                                             <  1:503175/140‑1 (MQ=255)
                               tGCTATGTTGTGTTGTCATAACAAATACAGTCCCCCTTTTTAAAATTAATAATGTAACGTTCTATAAATACTACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATCACTATc                                                                                                            <  1:358092/140‑1 (MQ=255)
                               tGCTATGTCGTTTTGTTATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATc                                                                                                            <  2:430600/140‑1 (MQ=255)
                                      ttttgttatCATAAAAAATACAGTCCACCTTTTTAAAATGAATAATGGAACGTTTGAGAAATATGACATATTAATACAACTTTAAACACATAAGTGTCTCAATAGATTTCCATGGTGATGCTCTATTAATATTAACTATg                                                                                                     >  2:158253/4‑140 (MQ=255)
                                         tgttgtCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGATATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGtat                                                                                                  <  1:639350/140‑1 (MQ=255)
                                         agttgtCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGATAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCAATGGTGATGCTCTATTCATATCTACGATGtat                                                                                                  >  1:592715/2‑140 (MQ=255)
                                          gttgtCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGtata                                                                                                 >  2:352439/1‑140 (MQ=255)
                                          gttgtCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACCATGaata                                                                                                 >  2:326791/1‑140 (MQ=255)
                                          gttgtCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACAAATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATAGTGATGCTCTATTACTTACTACTATGtata                                                                                                 >  1:283863/1‑140 (MQ=255)
                                           ttgtCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTatat                                                                                                >  1:545365/1‑140 (MQ=255)
                                           cggtCATAACAAATAGAGTCCCTCTTTTTAAAATGAATAATGGCACGTTCGAGAAATATGACATATTAATACCGCATTAAATCATTAAGTGTCTCATGAGGTTCGCGTGGGGATGCTTTATTCATATCTACGATGgatag                                                                                                >  2:7487/3‑139 (MQ=255)
                                             gtCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCCATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTa                                                                                              <  2:78917/140‑1 (MQ=255)
                                             gtCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCCATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTa                                                                                              <  2:573247/140‑1 (MQ=255)
                                                aTAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGTAATACGACATATTAATTTCTTTTTTTTTTTTTTTGTGTTTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAata                                                                                           <  1:39054/140‑1 (MQ=255)
                                                  aaCAATAACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATAtt                                                                                         <  2:27577/140‑1 (MQ=255)
                                                       aTACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTtata                                                                                    <  1:584258/140‑1 (MQ=255)
                                                       aTACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTtata                                                                                    >  2:230043/1‑140 (MQ=255)
                                                        tACAGTCCCCCTTTTTAAAATTAATAATCGAACGTTCGAGAAATTCGACATATTAATACCACTTTAAATCAATAAGTGTCTCCATCGGTTCCCATGTTGATGCTCTATTACTATCTACTATGTATATTAATATTTTataa                                                                                   <  2:317096/140‑1 (MQ=255)
                                                          cAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTTACTCAATAGGTTCCCATGGTGATGCTCTATTAATATCTAATATGTATATTAATATTTTataata                                                                                 >  1:41185/1‑140 (MQ=255)
                                                                 ccTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCCCTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAac                                                                          <  1:550992/140‑1 (MQ=255)
                                                                     tttAATATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAactact                                                                      <  2:11045/140‑1 (MQ=255)
                                                                        aaaaTGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTCCTTTa                                                                   >  1:571810/1‑140 (MQ=255)
                                                                        aaaaTGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTa                                                                   >  2:500170/1‑140 (MQ=255)
                                                                        aaaaTGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTAATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTa                                                                   >  1:275900/1‑140 (MQ=255)
                                                                         aaaTGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTaa                                                                  >  2:38323/1‑140 (MQ=255)
                                                                          aaTGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTAAt                                                                 >  1:272854/1‑140 (MQ=255)
                                                                            tGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTCAATATTTGAGATCTCCAATGGTTATATATTAAGATACTCTAATACGATTAACTAGTTATATTAATAGTTTATAATAACGAAACTACTTTAATta                                                               >  1:319434/1‑140 (MQ=255)
                                                                                   tGGAACCTTCGAGAAATACGACCTATTAAGACCACTTTAAATCAATAACTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGACACTACTTTAATTATATACtt                                                        <  2:409264/140‑1 (MQ=255)
                                                                                    ggAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTATCCATAGGATCTCCTGATGTCGCTCTTTCACTATCTCCTATCGAAAATATTCCTTTCAAATAACTAATCTACATGAATcatatttttt                                                       >  2:386534/1‑135 (MQ=255)
                                                                                        cGTTCGAGAAATTCGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTAAAAAAACGAAACTACTTTAATTATATACTTTAtgt                                                   <  1:17690/140‑1 (MQ=255)
                                                                                             gagaAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTAGTTATATACTTTATGTGCgtt                                              >  2:587297/1‑138 (MQ=255)
                                                                                             gagaAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTAATTTAATTATATACTTTATGTGCGAt                                              >  2:362036/1‑140 (MQ=255)
                                                                                              agaAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTAATTATATACTTTATGTGCGATa                                             <  1:628553/140‑1 (MQ=255)
                                                                                                aaaTACGATATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCTATGGTTATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTAATTATATACTTTATGTGCGATAGc                                           <  2:590184/140‑1 (MQ=255)
                                                                                                   tACGCCATATTAATATCACTTTAAATCAATAAGTGTTTCAATAGGTTCCCATGGTGATGCTCTATTACTATTTACTTTTTATATTAATATTTTATAATAACGAAACTACTTTAATTATATACTTTATGTGCGATAGcaca                                        <  1:226716/140‑1 (MQ=255)
                                                                                                   tACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTCCAAGTTTGAACTCCTTCTTAAGATTAACTATGTATATTAATCTGGTCTAAGAACAAAACTACTATAATTATATACTTTATGTGCAATAGcaca                                        >  1:656804/1‑140 (MQ=255)
                                                                                                           atTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTAATTATATACTTTATGTGCGATAGCACAAAGACaaa                                <  1:38323/140‑1 (MQ=255)
                                                                                                           atTAATACCACTTTAAATCAATAAGTATCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTAATTATATACTTTATGTGCGATAGCACAAAGACaaa                                <  1:37963/140‑1 (MQ=255)
                                                                                                                tACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTAATTATATACTTTATGTGCGATAGCACCAGGACAAAACcac                           >  2:382211/1‑140 (MQ=255)
                                                                                                                     cTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTAATTATATACTTTATGTGCGATAGAACAAAGACAAAACCACAAAcc                      >  1:665197/1‑140 (MQ=255)
                                                                                                                            tCAATAAGTGTCTCAATAGGTTCCCATGGTGATTCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTAATTATATACTTTATGTGCGATAGCACAAAGACAAAACCACAAACCTTAACAg               <  1:87584/140‑1 (MQ=255)
                                                                                                                             cAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTAATTATATACTTTATGTGCGATAGCACAAAGACAAAACCACAAACCTTCACCGc              >  2:322911/1‑140 (MQ=255)
                                                                                                                             cAATAAGTGTCCCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTAACTATATACTTTATGTGCGATAGCACATAGTTAATTCTTCAAACCTTAACAGc              <  2:473680/140‑1 (MQ=255)
                                                                                                                                      tctcAATAGGTTCCCATGGTGATGCTCTATTACTATATACTATGTATATTAATATTTTATAATAACGAAACTACTTTAATTATATACTTTATGTGCGATAGCACAAAGACAAAACCACGAACCTTAACAGCTATCAACCa     >  1:461779/1‑140 (MQ=255)
                                                                                                                                        tcAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTAATTATATAGTTTATGTACGATAGCACAAAGACAAAACCACAAACCTTAACAGCTATCAACCaat   >  1:492299/1‑140 (MQ=255)
                                                                                                                                          cAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTAATTATATACTTTATGTGCGATAGCACAAAGACAAAACCACAAACCTTAACAGCTATCAACCaatt  <  1:278006/140‑1 (MQ=255)
                                                                                                                                          cAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAAGAACGAAACTACTTTAATTATATACTTTATGTGCGATAGCACAAAGACAAAACCACAAACCTTAGAAGCTATCAACCagtt  >  2:447432/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                           
ATACCGCCACCAATTAAGATGACATCTGTTTTGCTATGTTGTGTTGTCATAACAAATACAGTCCCCCTTTTTAAAATGAATAATGGAACGTTCGAGAAATACGACATATTAATACCACTTTAAATCAATAAGTGTCTCAATAGGTTCCCATGGTGATGCTCTATTACTATCTACTATGTATATTAATATTTTATAATAACGAAACTACTTTAATTATATACTTTATGTGCGATAGCACAAAGACAAAACCACAAACCTTAACAGCTATCAACCAATT  >  USA300TCH1516_ALE/2485633‑2485909

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: