Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,245,338 A→T 8.0% intergenic (+143/‑100) fabG → / → acpP 3‑oxoacyl‑[acyl‑carrier‑protein] reductase FabG/Acyl carrier protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,245,3380AT 8.0% 784.4 / 36.9 240intergenic (+143/‑100)fabG/acpP3‑oxoacyl‑[acyl‑carrier‑protein] reductase FabG/Acyl carrier protein
Reads supporting (aligned to +/- strand):  ref base A (105/113);  new base T (4/16);  total (110/130)
Fisher's exact test for biased strand distribution p-value = 1.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

GAATTGTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTAAAGGAGG  >  USA300TCH1516_ALE/1245248‑1245428
                                                                                          |                                                                                          
gaATTGTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATaa                                                                                            >  2:377739/1‑91 (MQ=255)
 aaTTGTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATaaa                                                                                           >  1:1924844/1‑91 (MQ=255)
 aaTTGTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATaaa                                                                                           <  2:3610302/91‑1 (MQ=255)
  aTTGTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAAt                                                                                          >  1:2038663/1‑91 (MQ=255)
   ttGTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAAtt                                                                                         <  1:71103/91‑1 (MQ=255)
   ttGTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAAtt                                                                                         <  2:2849466/91‑1 (MQ=255)
     gTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGa                                                                                       <  2:1686608/91‑1 (MQ=255)
       cTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATa                                                                                     >  2:2614888/1‑91 (MQ=255)
        tGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTTATCAATTTATTGAtta                                                                                    <  1:537471/91‑3 (MQ=38)
        tGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATaa                                                                                    >  1:1077950/1‑91 (MQ=255)
         gACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTTATCAATTTATTGAttag                                                                                   <  2:617854/91‑4 (MQ=38)
         gACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTTATCAATTTATTGATTAg                                                                                   <  1:3464800/91‑1 (MQ=38)
          aCCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTTATCAATTTATTGATTAgt                                                                                  <  2:2642433/91‑1 (MQ=38)
          aCCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTTATCAATTTATTGATTAgt                                                                                  <  2:2025242/91‑1 (MQ=38)
           ccTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTTATCAATTTATTGATTAgtg                                                                                 <  2:472225/91‑1 (MQ=38)
            cTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTTATCAATTTATTGATTAgtgt                                                                                <  2:356773/91‑1 (MQ=37)
               gtcaacTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTc                                                                             >  1:1033630‑M2/7‑91 (MQ=255)
               gtcaacTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTc                                                                             >  1:2970382‑M2/7‑91 (MQ=255)
                 caacTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTa                                                                           >  1:269434‑M2/5‑91 (MQ=255)
                  aacTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTtctta                                                                          <  1:368223‑M2/88‑3 (MQ=255)
                   acTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGAGTTTCTaaa                                                                         >  1:2872264‑M2/3‑91 (MQ=255)
                     tttGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAAt                                                                       >  1:2040410/1‑91 (MQ=255)
                       tGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAAttt                                                                     <  1:2787008/91‑1 (MQ=255)
                         cGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCt                                                                   >  2:2328023/1‑91 (MQ=255)
                           ggggAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTAc                                                                 <  2:2811772/91‑1 (MQ=255)
                           ggggAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTAc                                                                 <  2:3871413/91‑1 (MQ=255)
                           ggggAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTAc                                                                 <  2:765022/91‑1 (MQ=255)
                             ggAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACtt                                                               <  1:2909321/91‑1 (MQ=255)
                              gAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTg                                                              >  1:3014965/1‑91 (MQ=255)
                               aaaTTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGt                                                             <  2:1080669/91‑1 (MQ=255)
                                  ttCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGtttt                                                          >  2:1295107/1‑91 (MQ=255)
                                  ttCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGACAAGTGTTTCTAAAATTTCTACTTGtttt                                                          >  1:1131471/1‑91 (MQ=255)
                                   tCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGttttt                                                         >  1:1549632/1‑91 (MQ=255)
                                   tCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGttttt                                                         <  1:2714120/91‑1 (MQ=255)
                                     tAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTa                                                       <  1:1200492/91‑1 (MQ=255)
                                      aaGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAg                                                      <  2:2159821/91‑1 (MQ=255)
                                       aGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATACATTGATAAGTGTTTCTAAAATTTCTACTTGTTCTTTAGa                                                     >  2:3346739/1‑91 (MQ=255)
                                       aGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGa                                                     >  1:3249660/1‑91 (MQ=255)
                                       aGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGa                                                     >  1:1144735/1‑91 (MQ=255)
                                       aGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGa                                                     >  2:3346605/1‑91 (MQ=255)
                                        gCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGaa                                                    <  2:599285/91‑1 (MQ=255)
                                         cAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAAt                                                   <  1:3865642/91‑1 (MQ=255)
                                          aaCCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAAtt                                                  >  2:3786481/1‑91 (MQ=255)
                                           aCCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAAttt                                                 <  1:2505962/91‑1 (MQ=255)
                                           aCCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAAttt                                                 >  1:1727817/1‑91 (MQ=255)
                                             cTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTaa                                               <  2:296031/91‑1 (MQ=255)
                                              tAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTaaa                                              <  1:3185303/91‑1 (MQ=255)
                                                gATAAGGTTGCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAAt                                            <  1:3787254/91‑1 (MQ=255)
                                                gATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAAt                                            >  1:2995360/1‑91 (MQ=255)
                                                gATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAAt                                            >  1:2540456/1‑91 (MQ=255)
                                                 aTAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATg                                           >  2:3817978/1‑91 (MQ=255)
                                                   aaGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTGCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATggg                                         <  2:501000/91‑1 (MQ=255)
                                                    aGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGa                                        >  1:3095426/1‑91 (MQ=255)
                                                     ggTTCCAGAATTTTTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTGTTTAGAATTTAAAATGGGaa                                       <  2:241796/91‑1 (MQ=255)
                                                       ttCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGaaaa                                     <  2:2478387/91‑1 (MQ=255)
                                                        tCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAat                                    <  1:3868647/91‑1 (MQ=255)
                                                         ccAGAATTTCTCCCTAAGAAACACTAATCAATACATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTGGAATTTAAAATGGGAAAata                                   <  1:3236393/91‑1 (MQ=255)
                                                         ccAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAata                                   >  1:151041/1‑91 (MQ=255)
                                                         ccAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAata                                   >  2:2679951/1‑91 (MQ=255)
                                                           aGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAAtata                                 <  2:95115/91‑1 (MQ=255)
                                                           aGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAAtata                                 <  2:2102186/91‑1 (MQ=255)
                                                               tttCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATAtagtag                             >  2:138248/1‑91 (MQ=255)
                                                               tttCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATAtagtag                             >  2:1759150/1‑91 (MQ=255)
                                                                 tctcCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTc                           <  2:1467510/91‑1 (MQ=255)
                                                                  ctcCCTAAGAAACACTAATCAATAAATTGATAAGTGGTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCt                          <  2:511368/91‑1 (MQ=255)
                                                                   tcCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTa                         >  1:478645/1‑91 (MQ=255)
                                                                    cccTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTAt                        >  1:68501/1‑91 (MQ=255)
                                                                     ccTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATg                       <  2:2123387/91‑1 (MQ=255)
                                                                     ccTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATg                       >  1:2858776/1‑91 (MQ=255)
                                                                       tAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGta                     <  1:2024883/91‑1 (MQ=255)
                                                                       tAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGta                     <  1:185669/91‑1 (MQ=255)
                                                                         aGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGtata                   <  2:3727134/91‑1 (MQ=255)
                                                                         aGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGtata                   <  1:361432/91‑1 (MQ=255)
                                                                          gAAACACTTATCAATTTATTGATTAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAg                  >  2:3526619/1‑91 (MQ=37)
                                                                          gAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAg                  <  1:2158686/91‑1 (MQ=255)
                                                                          gAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAg                  <  1:2158740/91‑1 (MQ=255)
                                                                            aaCACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGc                >  2:2608285/1‑91 (MQ=255)
                                                                            aaCACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGc                <  2:1279312/91‑1 (MQ=255)
                                                                             acacTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCa               <  2:3955696/91‑1 (MQ=255)
                                                                             acacTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCa               <  2:1943438/91‑1 (MQ=255)
                                                                             acacTAAACAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCa               <  1:1662335/91‑1 (MQ=255)
                                                                              cacTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAAGCAt              <  2:443106/91‑1 (MQ=255)
                                                                               acTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCAtt             >  1:1158934/1‑91 (MQ=255)
                                                                               acTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCAtt             >  1:504187/1‑91 (MQ=255)
                                                                                cTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCAttt            >  2:821744/1‑91 (MQ=255)
                                                                                 tAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCAtttt           >  2:1096847/1‑91 (MQ=255)
                                                                                  aaTCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCAttttt          <  1:2338618/91‑1 (MQ=255)
                                                                                  aaTCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCAttttt          >  2:627278/1‑91 (MQ=255)
                                                                                   aTCAATTTATTGATTAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGACAATATAGTAGTCTATGTATAGGCATTTTTa         >  2:2800382/1‑91 (MQ=39)
                                                                                    tCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTaa        >  1:2944788/1‑91 (MQ=255)
                                                                                     cAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTaaa       <  1:973170/91‑1 (MQ=255)
                                                                                      aaTAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTAAag      <  2:3162268/91‑1 (MQ=255)
                                                                                       aTAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTAAagg     >  1:245215/1‑91 (MQ=255)
                                                                                        tAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTAAagga    >  1:1426574/1‑91 (MQ=255)
                                                                                        tAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTAAagga    <  2:2575485/91‑1 (MQ=255)
                                                                                        tAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTAAagga    >  2:2241630/1‑91 (MQ=255)
                                                                                          aaTTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTAAaggagg  >  1:3593402/1‑91 (MQ=255)
                                                                                          aaTTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTAAaggagg  <  1:2286004/91‑1 (MQ=255)
                                                                                          aaTTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTAAaggagg  <  2:2019787/91‑1 (MQ=255)
                                                                                          aaTTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGGCTATGTATAGGCATTTTTAAaggagg  <  2:3378290/91‑1 (MQ=255)
                                                                                          |                                                                                          
GAATTGTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTAAAGGAGG  >  USA300TCH1516_ALE/1245248‑1245428

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