Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 115,765 T→A 12.7% intergenic (+32/‑119) btr → / → yxeP_1 HTH‑type transcriptional activator Btr/putative hydrolase YxeP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE115,7650TA12.7% 670.8 / 71.4 220intergenic (+32/‑119)btr/yxeP_1HTH‑type transcriptional activator Btr/putative hydrolase YxeP
Reads supporting (aligned to +/- strand):  ref base T (91/101);  new base A (22/6);  total (113/107)
Fisher's exact test for biased strand distribution p-value = 2.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 225 total aligned reads displayed.

TTAAAATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTAATCAT  >  USA300TCH1516_ALE/115675‑115855
                                                                                           |                                                                                          
 taaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATa                                                                                           >  2:3308352/1‑91 (MQ=255)
 taaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATa                                                                                           >  1:1565641/1‑91 (MQ=255)
     aTGATTACGAATGGGATTAACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAAttt                                                                                       >  2:3063947/1‑91 (MQ=255)
        aTTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgtg                                                                                    >  2:1630767/1‑91 (MQ=255)
        aTTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgtg                                                                                    >  1:2478646/1‑91 (MQ=255)
        aTTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTATGTATATCAAAATTtgtg                                                                                    <  1:1082539/91‑1 (MQ=37)
          tACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAt                                                                                  >  1:2083565/1‑91 (MQ=255)
           aCGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtt                                                                                 <  1:1205158/91‑1 (MQ=255)
           aCGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTATGTATATCAAAATTTGTGAtt                                                                                 <  1:2852620/91‑1 (MQ=25)
            cGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAttt                                                                                >  1:1338398/1‑91 (MQ=255)
            cGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTATGTATATCAAAATTTGTGAttt                                                                                <  1:1762629/91‑1 (MQ=25)
             gAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtttt                                                                               >  1:164839/1‑91 (MQ=255)
             gAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTATGTATATCAAAATTTGTGAtttg                                                                               <  1:2266492/91‑2 (MQ=25)
              aaTGGTATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAttttt                                                                              <  2:3143489/91‑1 (MQ=255)
               ctGGGATTCACTTAAGTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGCTATACATAATTTGTGAtttttt                                                                             <  2:3853980/90‑1 (MQ=255)
                 gggATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTtat                                                                           <  2:2943428/91‑1 (MQ=255)
                   gATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTatat                                                                         <  1:1987642/91‑1 (MQ=255)
                    aTTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATAtt                                                                        <  2:1344044/91‑1 (MQ=255)
                      tCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCa                                                                      <  2:525566/91‑1 (MQ=255)
                      tCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCa                                                                      <  1:690772/91‑1 (MQ=255)
                      tCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCa                                                                      >  1:2406748/1‑91 (MQ=255)
                       cACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCaa                                                                     <  2:263199/91‑1 (MQ=255)
                       cACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCaa                                                                     >  2:45953/1‑91 (MQ=255)
                        aCTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCaaa                                                                    <  2:1388093/91‑1 (MQ=255)
                         cTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAg                                                                   >  1:2973576/1‑91 (MQ=255)
                         cTTAATTGAAGTAAAACGTCACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAg                                                                   >  2:3393416/1‑91 (MQ=255)
                          ttAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATACACATAATTTGTGATTTTTTGTATTCAAAgt                                                                  <  2:3809001/91‑2 (MQ=255)
                           tAATTGAAGTAATACGTTACCCAAGTTCATAAAATGCTCACAAATCACATATTGTGATATACATAATTTGTGATTTTTTATATTCAAAGCt                                                                 <  2:2053375/91‑1 (MQ=255)
                           tAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCt                                                                 <  1:2017433/91‑1 (MQ=255)
                            aaTTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTa                                                                >  1:752812/1‑91 (MQ=255)
                              ttGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTaaa                                                              >  2:3333822/1‑91 (MQ=255)
                                gAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAt                                                            <  1:1061508/91‑1 (MQ=255)
                                 aaGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAtt                                                           <  2:1697878/91‑1 (MQ=255)
                                 aaGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAtt                                                           >  2:1878296/1‑91 (MQ=255)
                                 aaGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAtt                                                           >  1:2938174/1‑91 (MQ=255)
                                 aaGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAtt                                                           >  1:641914/1‑91 (MQ=255)
                                      aaaCGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaa                                                      >  1:845764/1‑91 (MQ=255)
                                         cGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaaaaa                                                   <  2:792544/91‑1 (MQ=255)
                                         cGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaaaaa                                                   >  2:1725461/1‑91 (MQ=255)
                                         cGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaaaaa                                                   >  1:3644793/1‑91 (MQ=255)
                                          gTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAt                                                  >  2:1109687/1‑91 (MQ=255)
                                          gTTACCCAAGTTCATAAAATGATCACAAATCACAAATGCGGATATACATAATTTGGGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAt                                                  >  2:1556423/1‑91 (MQ=255)
                                           ttACGCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAtt                                                 <  1:3234926/91‑1 (MQ=255)
                                           ttACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAtt                                                 <  1:3038096/91‑1 (MQ=255)
                                           ttACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATAGACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAtt                                                 <  1:3881515/91‑1 (MQ=255)
                                            tACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTa                                                <  2:1641924/91‑1 (MQ=255)
                                             aCCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                               >  2:1844872/1‑91 (MQ=255)
                                             aCCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                               >  1:763347/1‑91 (MQ=255)
                                               ccAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATg                                             <  1:766789/91‑1 (MQ=255)
                                                cAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATgg                                            <  2:1969698/91‑1 (MQ=255)
                                                 aaGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGt                                           >  2:2652217/1‑91 (MQ=255)
                                                  aGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGtt                                          <  2:404771/91‑1 (MQ=255)
                                                   gTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTa                                         <  2:2184353/91‑1 (MQ=255)
                                                   gTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTa                                         >  1:3101092/1‑91 (MQ=255)
                                                   gTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTa                                         >  2:3676136/1‑91 (MQ=255)
                                                    ttCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTaa                                        >  2:2009034/1‑91 (MQ=255)
                                                    ttCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTaa                                        >  1:665185/1‑91 (MQ=255)
                                                       aTAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAt                                     >  1:2667986/1‑91 (MQ=255)
                                                       aTAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAt                                     >  1:2667391/1‑91 (MQ=255)
                                                        tAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtc                                    >  2:626233/1‑91 (MQ=255)
                                                        tAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtc                                    >  1:3770566/1‑91 (MQ=255)
                                                          aaaTGATCACGAATCACTAATTTTGATATACATAATTTGAGATTTTTTATATTCAAAGCTACAATTGCAAAAAATTAATGGTTAACAtctc                                  <  2:541553/91‑1 (MQ=255)
                                                          aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                                  >  1:3433707/1‑91 (MQ=255)
                                                          aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                                  >  1:2381862/1‑91 (MQ=255)
                                                          aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                                  >  1:778988/1‑91 (MQ=255)
                                                          aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                                  <  2:28660/91‑1 (MQ=255)
                                                           aaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAAAGAAAAAAATTAATGGTTAACATctct                                 <  2:2372992/91‑1 (MQ=255)
                                                            aTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtg                                >  1:252784/1‑91 (MQ=255)
                                                               atcacaaatcacaaATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttg                             <  1:2255082/91‑1 (MQ=255)
                                                                tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                            >  2:1279180/1‑91 (MQ=255)
                                                                tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                            >  2:1867202/1‑91 (MQ=255)
                                                                taaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                            >  1:1843626/5‑91 (MQ=21)
                                                                taaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                            >  1:876616/5‑91 (MQ=21)
                                                                taaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                            >  1:1728891/5‑91 (MQ=21)
                                                                taaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                            >  2:2889556/5‑91 (MQ=21)
                                                                 aaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                           >  2:132486/4‑91 (MQ=21)
                                                                 aaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                           >  1:280449/4‑91 (MQ=21)
                                                                 aaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                           >  1:3800803/4‑91 (MQ=21)
                                                                  acaaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGttgttg                          >  1:2600937/1‑91 (MQ=255)
                                                                  aaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGttgttg                          >  2:2064120/3‑91 (MQ=21)
                                                                   caaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTgg                         <  1:1319543/91‑1 (MQ=255)
                                                                   caaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTgg                         >  2:2366776/1‑91 (MQ=255)
                                                                      atcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCaa                      <  2:1427710/91‑1 (MQ=255)
                                                                      atcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCaa                      <  2:1427893/91‑1 (MQ=255)
                                                                      atcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCaa                      >  2:495963/1‑91 (MQ=25)
                                                                      atcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCaa                      >  2:3497553/1‑91 (MQ=25)
                                                                       tcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAat                     >  1:1751342/1‑91 (MQ=37)
                                                                         acaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAatat                   <  1:3665898/91‑1 (MQ=255)
                                                                          caaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAAtata                  <  2:1518161/91‑1 (MQ=255)
                                                                          caaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAAtata                  <  1:3497972/91‑1 (MQ=255)
                                                                          caaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAAtata                  <  1:3333258/91‑1 (MQ=255)
                                                                           aaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataa                 <  2:3720570/91‑1 (MQ=255)
                                                                               tttGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaataa             <  1:85758/91‑1 (MQ=255)
                                                                                 tGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATg           >  1:1188506/1‑91 (MQ=255)
                                                                                  gATATACATAATTTGTGATTTTTTATATGCAAAGCTAAAATTGCAAAAAAGTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGa          <  2:1648539/91‑1 (MQ=255)
                                                                                    tataCATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGAtt        <  2:2735497/91‑1 (MQ=255)
                                                                                    tataCATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGAtt        >  1:444857/1‑91 (MQ=255)
                                                                                      taCATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTaa      <  2:2535966/91‑1 (MQ=255)
                                                                                      taCATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTaa      <  1:559792/91‑1 (MQ=255)
                                                                                       aCATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTAAt     <  2:3632706/91‑1 (MQ=255)
                                                                                           |                                                                                          
TTAAAATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTAATCAT  >  USA300TCH1516_ALE/115675‑115855

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