Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 115,762:1 +T 6.1% intergenic (+29/‑122) 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,7621.T 6.1% 946.8 / 50.9 228intergenic (+29/‑122)btr/yxeP_1HTH‑type transcriptional activator Btr/putative hydrolase YxeP
Reads supporting (aligned to +/- strand):  ref base . (102/112);  new base T (9/5);  total (111/117)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 242 total aligned reads displayed.

ACATTAAAATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTAAT  >  USA300TCH1516_ALE/115672‑115852
                                                                                           |                                                                                          
acattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGAtata                                                                                             >  1:1237128/1‑91 (MQ=255)
acattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGAtata                                                                                             >  2:272218/1‑91 (MQ=255)
acattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGAcata                                                                                             >  1:798551/1‑91 (MQ=255)
 cattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATAc                                                                                           <  1:1719069/91‑1 (MQ=255)
  attaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACa                                                                                          >  2:2634698/1‑91 (MQ=255)
     aaaaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATaa                                                                                       >  1:247423/1‑91 (MQ=255)
     aaaaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATaa                                                                                       >  2:1442506/1‑91 (MQ=255)
       aaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAAtt                                                                                     >  1:559066/1‑91 (MQ=255)
       aaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAAtt                                                                                     <  1:3301088/91‑1 (MQ=255)
       aaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAAtt                                                                                     >  1:403163/1‑91 (MQ=255)
        aTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACGCAAGTTCATCAAATGATCACAAATCACAATTTTTGATATACATAAttt                                                                                    <  2:3796744/91‑1 (MQ=255)
        aTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAAttt                                                                                    >  1:2006886/1‑91 (MQ=255)
         tGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTATGTATATCAAAATTtg                                                                                   <  1:3206764/91‑1 (MQ=37)
         ggATGACGAAATGGATTCTCTTAATTGAAGTAATAGGGAACCCAAATTCATTAAATGAGCACAAATCACAAATTATGTATATCATAATTtg                                                                                   <  2:310274/90‑1 (MQ=12)
          gATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgt                                                                                  <  2:499822/91‑1 (MQ=255)
          gATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgt                                                                                  <  1:2456425/91‑1 (MQ=255)
           aTTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgtg                                                                                 <  1:1555107/91‑1 (MQ=255)
           aTTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgtg                                                                                 <  1:1177130/91‑1 (MQ=255)
           aTTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgtg                                                                                 >  2:3267883/1‑91 (MQ=255)
            ttACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGa                                                                                <  2:3853140/91‑1 (MQ=255)
            ttACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGa                                                                                >  2:2909005/1‑91 (MQ=255)
            ttACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGa                                                                                >  2:3078795/1‑91 (MQ=255)
            ttACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGa                                                                                >  1:1081252/1‑91 (MQ=255)
            ttACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGa                                                                                >  1:2752064/1‑91 (MQ=255)
             tACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAt                                                                               >  1:1814956/1‑91 (MQ=255)
              aCGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtt                                                                              >  2:2081656/1‑91 (MQ=255)
               cGAATGGTATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAttt                                                                             <  2:3521354/91‑1 (MQ=255)
                gAATGGGATTAACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtttt                                                                            <  2:389682/91‑1 (MQ=255)
                   tGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTta                                                                         >  1:1576831/1‑91 (MQ=255)
                        ttCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTc                                                                    >  1:229053/1‑91 (MQ=255)
                         tCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCa                                                                   <  2:3415212/91‑1 (MQ=255)
                              tAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCt                                                              <  2:918149/91‑1 (MQ=255)
                               aaTTGTAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTa                                                             <  2:3634993/91‑1 (MQ=255)
                               aaTTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTa                                                             >  1:385875/1‑91 (MQ=255)
                                   gAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAt                                                         <  2:1244120/91‑1 (MQ=255)
                                   gAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAt                                                         <  1:3396856/91‑1 (MQ=255)
                                    aaGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATCCATAATTTGTGATTTTTTATATTCAAAGCTAAAAtt                                                        >  2:3117801/1‑91 (MQ=255)
                                     aGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTg                                                       >  1:2916683/1‑91 (MQ=255)
                                        aaaaCGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaa                                                    <  2:112044/91‑1 (MQ=255)
                                        aaaaCGTTACCCAAGTTCAAAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaa                                                    <  2:1984053/91‑1 (MQ=255)
                                           aCGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaaaa                                                 >  1:3058862/1‑91 (MQ=255)
                                            cGTTAGCCAAGTTCATAAAATGATCACAAATCACGACTTTTGATATACGTAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaaaaa                                                <  2:224176/91‑1 (MQ=255)
                                            agTTACCCAAGTTCATAAAAAGATCAAAAATCACAAATTTTGATATACATAATTTGTAATATTTAATAATAAAAGAAAAAATAaaaaaaaa                                                <  2:2049151/90‑1 (MQ=255)
                                             gTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAt                                               >  2:1058247/1‑91 (MQ=255)
                                              ttACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAtt                                              >  1:1473219/1‑91 (MQ=255)
                                               tACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTa                                             >  1:2441559/1‑91 (MQ=255)
                                                aCCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                            >  1:272983/1‑91 (MQ=255)
                                                aCCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                            >  1:2569514/1‑91 (MQ=255)
                                                aCCCAAGTAAATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                            <  2:1234667/91‑1 (MQ=255)
                                                  ccAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATg                                          >  2:2585775/1‑91 (MQ=255)
                                                    aaGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGt                                        <  1:161970/91‑1 (MQ=255)
                                                     aGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGtt                                       >  1:2720526/1‑91 (MQ=255)
                                                      gTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTa                                      >  1:3522435/1‑91 (MQ=255)
                                                      gTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTa                                      <  2:2218184/91‑1 (MQ=255)
                                                       ttCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTaa                                     >  2:1261954/1‑91 (MQ=255)
                                                       ttCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTaa                                     <  2:1575646/91‑1 (MQ=255)
                                                        tCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAAc                                    >  2:588926/1‑91 (MQ=255)
                                                        tCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAAc                                    >  2:2546658/1‑91 (MQ=255)
                                                           tAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtc                                 >  1:1061980/1‑91 (MQ=255)
                                                           tAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtc                                 >  2:75741/1‑91 (MQ=255)
                                                            aaaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtct                                <  2:3701227/91‑1 (MQ=255)
                                                            aaaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtct                                <  2:1125695/91‑1 (MQ=255)
                                                            aaaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtct                                <  2:2350910/91‑1 (MQ=255)
                                                            aaaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtct                                <  1:957363/91‑1 (MQ=255)
                                                             aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                               >  1:2851701/1‑91 (MQ=255)
                                                             aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                               <  1:2865283/91‑1 (MQ=255)
                                                              aaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATctct                              <  1:1875445/91‑1 (MQ=255)
                                                                   tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                         >  1:2020922/1‑91 (MQ=255)
                                                                   tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                         <  2:2278904/91‑1 (MQ=255)
                                                                   tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                         >  1:1626306/1‑91 (MQ=255)
                                                                   taaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                         >  2:2408432/5‑91 (MQ=21)
                                                                    aaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                        >  2:1291983/4‑91 (MQ=21)
                                                                    aaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                        >  2:2926623/4‑91 (MQ=21)
                                                                     aaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGttgttg                       >  1:872084/3‑91 (MQ=21)
                                                                      caaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTgg                      <  1:1761096/91‑1 (MQ=255)
                                                                       aaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCTAAGCTAAAGTTGAAAAAAATTATTGGTTAACATCTCTGTTGTTGGc                     <  1:3155091/91‑1 (MQ=255)
                                                                       aaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGc                     >  2:1593482/1‑91 (MQ=255)
                                                                        aatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCa                    >  2:1814365/1‑91 (MQ=25)
                                                                         atcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCaa                   <  1:367097/91‑1 (MQ=255)
                                                                          tcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGAAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAat                  <  2:2436113/91‑1 (MQ=255)
                                                                            acaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAatat                >  1:1027658/1‑91 (MQ=255)
                                                                              aaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataa              <  1:982758/91‑1 (MQ=255)
                                                                               aatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaa             <  1:3043176/91‑1 (MQ=255)
                                                                                atTTTGATATACATAATTTGTGATTTTTTTTTTTTTTTGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaat            >  2:3735445/1‑91 (MQ=255)
                                                                                atTTTGATATACATAATTTGTGATTTTTTTTATTTTAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaat            >  2:2962862/1‑91 (MQ=255)
                                                                                atTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaat            >  2:2180779/1‑91 (MQ=255)
                                                                                atTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaat            >  1:2715978/1‑91 (MQ=255)
                                                                                   ttGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAaataat         >  1:315553/1‑91 (MQ=255)
                                                                                    tGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATg        <  1:3555955/91‑1 (MQ=255)
                                                                                    tGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATg        <  1:578906/91‑1 (MQ=255)
                                                                                    tGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATg        <  2:3945670/91‑1 (MQ=255)
                                                                                    tGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAACTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATg        <  2:3469105/91‑1 (MQ=255)
                                                                                     gATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGa       <  2:465540/91‑1 (MQ=255)
                                                                                     gATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGa       <  2:1776994/91‑1 (MQ=255)
                                                                                      atatACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGAt      >  2:3367787/1‑91 (MQ=255)
                                                                                      atatACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGAt      <  2:2467319/91‑1 (MQ=255)
                                                                                       tataCATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGAtt     <  2:535774/91‑1 (MQ=255)
                                                                                        ataCATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTa    <  1:2547892/91‑1 (MQ=255)
                                                                                        ataCATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTa    <  2:724163/91‑1 (MQ=255)
                                                                                          aCATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTAAt  >  2:3067158/1‑91 (MQ=255)
                                                                                           |                                                                                          
ACATTAAAATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTAAT  >  USA300TCH1516_ALE/115672‑115852

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