Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 629,681 C→A L84I (CTA→ATA)  USA300HOU_RS02990 → 3‑hexulose‑6‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE629,6810CA98.8% 317.8 / ‑3.7 86L84I (CTA→ATA) USA300HOU_RS029903‑hexulose‑6‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (35/50);  total (36/50)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

AAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGC  >  USA300TCH1516_ALE/629594‑629766
                                                                                       |                                                                                     
aaaGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAg                                                                                    <  1:29343/91‑1 (MQ=255)
aaaGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAg                                                                                    <  1:3048518/91‑1 (MQ=255)
   gTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGgtg                                                                                 <  1:2793740/91‑1 (MQ=255)
    tATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGgtgt                                                                                <  1:3011086/91‑1 (MQ=255)
    tATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGgtgt                                                                                <  1:3012139/91‑1 (MQ=255)
     ataaGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGtt                                                                               <  1:2666424/88‑1 (MQ=255)
     aTTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGtt                                                                               <  1:291352/91‑1 (MQ=255)
     aTTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGtt                                                                               <  1:1300164/91‑1 (MQ=255)
      ttAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTg                                                                              >  1:1428569/1‑91 (MQ=255)
      ttAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTg                                                                              >  1:3205171/1‑91 (MQ=255)
       tAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGc                                                                             <  1:1113408/91‑1 (MQ=255)
       tAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGc                                                                             <  1:1113419/91‑1 (MQ=255)
        aGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCa                                                                            <  1:1288320/91‑1 (MQ=255)
        aGCAGACATAAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGATGTTGCa                                                                            >  1:1850894/1‑91 (MQ=255)
         gCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCag                                                                           <  1:3105022/91‑1 (MQ=255)
          cAGACATGATAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCaga                                                                          >  1:1857751/1‑91 (MQ=255)
               aTGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATg                                                                     <  1:1358109/91‑1 (MQ=255)
               aTGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATg                                                                     >  1:1126228/1‑91 (MQ=255)
               aTGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATg                                                                     <  1:2569930/91‑1 (MQ=255)
               aTGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATg                                                                     <  1:2569889/91‑1 (MQ=255)
               aTGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATg                                                                     >  1:3656680/1‑91 (MQ=255)
               aTGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATg                                                                     <  1:2904408/91‑1 (MQ=255)
                 gAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCa                                                                   >  1:2008910/1‑91 (MQ=255)
                  aaaaTTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCAt                                                                  >  1:1783969/1‑91 (MQ=255)
                  aaaaTTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCAt                                                                  >  1:1784126/1‑91 (MQ=255)
                  aaaaTTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCAt                                                                  <  1:73784/91‑1 (MQ=255)
                  aaaaTTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCAt                                                                  <  1:2707853/91‑1 (MQ=255)
                    aaTTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCa                                                                >  1:2177134/1‑91 (MQ=255)
                    aaTTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCa                                                                >  1:1243091/1‑91 (MQ=255)
                    aaTTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCa                                                                >  1:3706700/1‑91 (MQ=255)
                      ttATGGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAAt                                                              >  1:1969123/1‑91 (MQ=255)
                        ataGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTa                                                            <  1:1463338/91‑1 (MQ=255)
                        ataGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTa                                                            >  1:1615902/1‑91 (MQ=255)
                          aGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTaaa                                                          >  1:467861/1‑91 (MQ=255)
                          aGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTaaa                                                          >  1:1865059/1‑91 (MQ=255)
                          aGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTaaa                                                          >  1:3198118/1‑91 (MQ=255)
                          aGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTaaa                                                          >  1:231971/1‑91 (MQ=255)
                          aGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTaaa                                                          >  1:232498/1‑91 (MQ=255)
                           gATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAag                                                         <  1:3372405/91‑1 (MQ=255)
                           gATGCAGCTGATTATGAAGGTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAag                                                         <  1:1658295/91‑1 (MQ=255)
                             tGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAagca                                                       <  1:2589811/91‑1 (MQ=255)
                                   tGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTAtt                                                 >  1:765942/1‑91 (MQ=255)
                                   tGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTAtt                                                 >  1:309981/1‑91 (MQ=255)
                                   tGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTAtt                                                 <  1:3126120/91‑1 (MQ=255)
                                   tGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTAtt                                                 <  1:998343/91‑1 (MQ=255)
                                      ttATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTgaa                                              <  1:2288784/91‑1 (MQ=255)
                                       tATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTgaag                                             >  1:1465660/1‑91 (MQ=255)
                                       tATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTgaag                                             >  1:1465597/1‑91 (MQ=255)
                                           aaGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGc                                         <  1:2766394/91‑1 (MQ=255)
                                           aaGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGc                                         <  1:945494/91‑1 (MQ=255)
                                            aGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCt                                        <  1:3098566/91‑1 (MQ=255)
                                             gTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTc                                       >  1:1559482/1‑91 (MQ=255)
                                              ttAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCa                                      >  1:3617525/1‑91 (MQ=255)
                                              ttAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCa                                      >  1:3616403/1‑91 (MQ=255)
                                              ttAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCa                                      >  1:3616776/1‑91 (MQ=255)
                                               tAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCAt                                     <  1:2051734/91‑1 (MQ=255)
                                               tAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCAt                                     <  1:1353043/91‑1 (MQ=255)
                                               tAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCAt                                     <  1:994245/91‑1 (MQ=255)
                                                  ccAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATaaa                                  <  1:3056349/91‑1 (MQ=255)
                                                  ccAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATaaa                                  <  1:491917/91‑1 (MQ=255)
                                                  ccAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATaaa                                  <  1:2359869/91‑1 (MQ=255)
                                                      gCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAaata                              >  1:2161019/1‑91 (MQ=255)
                                                      gCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAaata                              >  1:2162412/1‑91 (MQ=255)
                                                       cAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAaataa                             <  1:2378169/91‑1 (MQ=255)
                                                        aaTTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAaataat                            <  1:1862525/91‑1 (MQ=255)
                                                         aTTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAataata                           <  1:2337488/91‑1 (MQ=255)
                                                         aTTAAAGTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAataata                           <  1:2337644/91‑1 (MQ=255)
                                                          ttAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAAtaataa                          >  1:3409115/1‑91 (MQ=255)
                                                             aaTTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACa                       >  1:2547315/1‑91 (MQ=255)
                                                             aaTTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACa                       >  1:1075528/1‑91 (MQ=255)
                                                               tttGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGTTATTGAAGAAGCTCATAAAAATAATAAACAAt                     <  1:218061/91‑1 (MQ=255)
                                                               tttGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAAt                     <  1:2821081/91‑1 (MQ=255)
                                                                ttGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAAtt                    >  1:2077363/1‑91 (MQ=255)
                                                                ttGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAAtt                    >  1:1451495/1‑91 (MQ=255)
                                                                 tGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTa                   <  1:40333/91‑1 (MQ=255)
                                                                 tGGCGAGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGAAGCTATTGAAGAAGGTCATAAAAATAATAAACAATTa                   <  1:3115780/91‑1 (MQ=255)
                                                                    cgcgGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTa                <  1:3656516/91‑1 (MQ=255)
                                                                     gcgGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAg               <  1:1441608/91‑1 (MQ=255)
                                                                     gcgGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAg               <  1:2986312/91‑1 (MQ=255)
                                                                        gATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTg            >  1:1904354/1‑91 (MQ=255)
                                                                          tGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGat          <  1:3092044/91‑1 (MQ=255)
                                                                            tAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGatat        <  1:1212171/91‑1 (MQ=255)
                                                                             aaTTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATg       <  1:444480/91‑1 (MQ=255)
                                                                             aaTTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATg       <  1:444633/91‑1 (MQ=255)
                                                                             aaTTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATg       <  1:1302343/91‑1 (MQ=255)
                                                                                  cAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGc  >  1:2043296/1‑91 (MQ=255)
                                                                                       |                                                                                     
AAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGC  >  USA300TCH1516_ALE/629594‑629766

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