Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,049,714 A→C 5.3% M458L (ATG→CTG)  menD → 2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,049,7140AC 5.3% 646.0 / 23.3 190M458L (ATG→CTG) menD2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylate synthase
Reads supporting (aligned to +/- strand):  ref base A (101/79);  new base C (4/6);  total (105/85)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

Only 100 of 201 total aligned reads displayed.

GTGCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGAACGGTTGT  >  USA300TCH1516_ALE/1049624‑1049804
                                                                                          |                                                                                          
gtgCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGa                                                                                            >  1:44259/1‑91 (MQ=255)
gtgCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGa                                                                                            >  1:547242/1‑91 (MQ=255)
gtgCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGa                                                                                            <  2:1990186/91‑1 (MQ=255)
 tgCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGac                                                                                           >  2:3847037/1‑90 (MQ=255)
  gCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATg                                                                                          >  2:2080263/1‑91 (MQ=255)
     tAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAat                                                                                       <  1:611250/91‑1 (MQ=255)
      aaaCGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCATCTGAata                                                                                      <  2:1897113/91‑1 (MQ=255)
       aaCGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAatat                                                                                     <  1:1244731/91‑1 (MQ=255)
         cGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTg                                                                                   <  1:2460800/91‑1 (MQ=255)
          gAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGt                                                                                  <  2:1151017/91‑1 (MQ=255)
          gAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGt                                                                                  >  2:2704299/1‑91 (MQ=255)
          gAATAACATTATTGATAGGTGATTTATCATTTTACCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGt                                                                                  >  2:2185254/1‑91 (MQ=255)
           aaTAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGta                                                                                 <  1:2327089/91‑1 (MQ=255)
            aTAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGtat                                                                                <  1:2579762/91‑1 (MQ=255)
             tAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGtatt                                                                               >  1:392263/1‑91 (MQ=255)
                 attattGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTg                                                                           >  1:1027204/1‑91 (MQ=255)
                 attattGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTg                                                                           <  2:3765589/91‑1 (MQ=255)
                 attattGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTg                                                                           <  2:1733449/91‑1 (MQ=255)
                   tattGATAGGTGATTTATCATTTTATCATGATATGAAGGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGaa                                                                         <  1:62329/91‑1 (MQ=255)
                    attGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGaac                                                                        <  1:621867/91‑1 (MQ=255)
                    attGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGaac                                                                        >  2:1353027/1‑91 (MQ=255)
                    attGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGaac                                                                        >  2:1418475/1‑91 (MQ=255)
                      tGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGaacaa                                                                      <  2:2017813/91‑1 (MQ=255)
                      tGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGaacaa                                                                      <  1:108564/91‑1 (MQ=255)
                        aTAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACg                                                                    >  2:3432702/1‑91 (MQ=255)
                        aTAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACg                                                                    >  2:122970/1‑91 (MQ=255)
                         tAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGa                                                                   <  2:641906/91‑1 (MQ=255)
                          aGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGAt                                                                  >  2:1583544/1‑91 (MQ=255)
                          aGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGAt                                                                  >  2:2492142/1‑91 (MQ=255)
                          aGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGAt                                                                  <  1:451219/91‑1 (MQ=255)
                             tGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGAtggt                                                               <  2:3279936/91‑1 (MQ=255)
                              gATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGAtggtg                                                              <  2:2051302/91‑1 (MQ=255)
                               aTTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGAtggtgg                                                             <  2:3152164/91‑1 (MQ=255)
                               aTTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGAtggtgg                                                             <  2:844786/91‑1 (MQ=255)
                                tttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGc                                                            <  2:2588234/91‑1 (MQ=255)
                                tttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGc                                                            >  2:2627192/1‑91 (MQ=255)
                                tttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGc                                                            <  1:2976743/91‑1 (MQ=255)
                                 ttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGc                                                            >  1:313364/1‑90 (MQ=255)
                                 ttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCg                                                           >  1:1171624/1‑91 (MQ=255)
                                 ttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCg                                                           >  1:2488423/1‑91 (MQ=255)
                                 ttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCg                                                           <  2:1545172/91‑1 (MQ=255)
                                 ttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCg                                                           <  2:298766/91‑1 (MQ=255)
                                 ttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAAAAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCg                                                           >  1:3882256/1‑91 (MQ=255)
                                  tatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCgg                                                          <  1:1819793/91‑1 (MQ=255)
                                      attttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACGACGATGGTGGCGGTAtt                                                      >  2:458003/1‑91 (MQ=255)
                                      attttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtt                                                      >  2:457909/1‑91 (MQ=255)
                                        tttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtttt                                                    >  1:3438523/1‑91 (MQ=255)
                                        tttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtttt                                                    <  2:2169157/91‑1 (MQ=255)
                                        tttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtttt                                                    >  1:3439063/1‑91 (MQ=255)
                                          tatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtttttt                                                  >  1:406789/1‑91 (MQ=255)
                                          tatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtttttt                                                  <  2:158838/91‑1 (MQ=255)
                                            tcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCa                                                >  1:1429592/1‑91 (MQ=255)
                                             catcatatGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATATTTTCat                                               >  2:3443642/5‑91 (MQ=255)
                                             catGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCat                                               >  2:1218011/1‑91 (MQ=255)
                                              atGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCata                                              >  1:3372044/1‑91 (MQ=255)
                                                gATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtttttt                                                  >  2:2732766/1‑85 (MQ=255)
                                                gATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATAtt                                            >  2:1672125/1‑91 (MQ=255)
                                                 atatGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATAttt                                           >  2:545212/1‑91 (MQ=255)
                                                   atGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTAc                                         <  1:1486504/91‑1 (MQ=255)
                                                    tGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTAcc                                        >  1:2666263/1‑91 (MQ=255)
                                                        tGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACaa                                    >  1:2401993/1‑91 (MQ=255)
                                                         ggACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACaaa                                   <  2:1090090/91‑1 (MQ=255)
                                                          gACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACaaaa                                  <  1:226378/91‑1 (MQ=255)
                                                          gACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACaaaa                                  >  1:1387251/1‑91 (MQ=255)
                                                             tATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAaaaga                               >  1:1454527/1‑91 (MQ=255)
                                                             tATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAaaaga                               <  2:712065/91‑1 (MQ=255)
                                                             tATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAaaaga                               <  2:2421885/91‑1 (MQ=255)
                                                              aTTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAaaagaa                              >  2:3105186/1‑91 (MQ=255)
                                                              aTTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAaaagaa                              >  2:2190193/1‑91 (MQ=255)
                                                                  aTGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAaaag                                >  2:2286913/1‑85 (MQ=255)
                                                                    gTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCa                        >  1:157916/1‑91 (MQ=255)
                                                                    gTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCa                        >  1:2856641/1‑91 (MQ=255)
                                                                     tCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCaa                       <  1:3460359/91‑1 (MQ=255)
                                                                     tCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCaa                       <  1:837637/91‑1 (MQ=255)
                                                                      cAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACaaa                                   >  1:2733455/1‑78 (MQ=255)
                                                                      cAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACaaa                                   >  1:1666557/1‑78 (MQ=255)
                                                                      cAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAAc                      <  2:1836343/91‑1 (MQ=255)
                                                                      cAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAAc                      >  1:3822527/1‑91 (MQ=255)
                                                                      cAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAAc                      >  1:2542469/1‑91 (MQ=255)
                                                                       aaaaTTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGTTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACt                     <  1:3275019/91‑1 (MQ=255)
                                                                       aaaaTTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACt                     >  1:1540240/1‑91 (MQ=255)
                                                                       aaaaTTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACt                     <  2:349647/91‑1 (MQ=255)
                                                                        aaaTTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTg                    <  1:3336465/91‑1 (MQ=255)
                                                                          aTTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGAc                  >  1:3848533/1‑91 (MQ=255)
                                                                           ttAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACt                 <  2:1711474/91‑1 (MQ=255)
                                                                              aataatATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAtt              <  1:2237739/91‑1 (MQ=255)
                                                                              aataatATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAtg              >  2:3849700/1‑90 (MQ=255)
                                                                               ataataTTCGGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAttt             <  1:3713272/91‑1 (MQ=255)
                                                                               ataataTTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAttt             >  1:2306727/1‑91 (MQ=255)
                                                                                taataTTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAtttt            >  1:1946116/1‑91 (MQ=255)
                                                                                taataTTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAtttt            >  2:3686857/1‑91 (MQ=255)
                                                                                caataTTCATCTGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAtttt            >  1:2906081/2‑91 (MQ=255)
                                                                                 aataTTCATCTGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTg           >  2:2887323/1‑91 (MQ=255)
                                                                                  atatTCATCTGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGa          >  2:3551210/1‑91 (MQ=255)
                                                                                  atatTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGa          <  2:3618206/91‑1 (MQ=255)
                                                                                   tatTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGaa         <  1:363662/91‑1 (MQ=255)
                                                                                     ttCAGATGAGTATTGTATTATTGAACAACGATGGTGGCGGGATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGAACg       <  1:725718/91‑1 (MQ=255)
                                                                                     ttCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGAACg       <  1:2349967/91‑1 (MQ=255)
                                                                                          aTGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTTTTTTTTTTTTCAAAAAGAAAGTGCAACTGACTATTTTGATCGgttgt  >  2:3734475/1‑91 (MQ=255)
                                                                                          aTGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGAACGgttgt  <  1:3857385/91‑1 (MQ=255)
                                                                                          |                                                                                          
GTGCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGAACGGTTGT  >  USA300TCH1516_ALE/1049624‑1049804

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