Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 115,758 Δ1 bp 6.1% intergenic (+25/‑126) 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,7580A. 6.1% 768.7 / 36.1 231intergenic (+25/‑126)btr/yxeP_1HTH‑type transcriptional activator Btr/putative hydrolase YxeP
Reads supporting (aligned to +/- strand):  ref base A (105/111);  new base . (9/5);  total (114/117)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 241 total aligned reads displayed.

TTAAACATTAAAATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGAT  >  USA300TCH1516_ALE/115668‑115848
                                                                                          |                                                                                           
  aaacattaaaATGATTACGAATGGGATTCACTTAATTGATGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGata                                                                                           >  2:1501803/1‑91 (MQ=255)
  aaacattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGata                                                                                           >  1:2249895/1‑91 (MQ=255)
    acattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGAtata                                                                                         <  1:533247/91‑1 (MQ=255)
    acattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGAtata                                                                                         >  1:1237128/1‑91 (MQ=255)
    acattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGAtata                                                                                         >  1:3203898/1‑91 (MQ=255)
    acattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGAtata                                                                                         >  2:272218/1‑91 (MQ=255)
     cattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATAc                                                                                       <  2:2631052/91‑1 (MQ=255)
      attaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACa                                                                                      >  2:2785057/1‑91 (MQ=255)
       ttaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACAt                                                                                     <  2:1121115/91‑1 (MQ=255)
         aaaaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATaa                                                                                   >  1:247423/1‑91 (MQ=255)
         aaaaTGATTACGAATGGGATTCACTTAACTGAAGTAAAACGTTACCCAAGTTCATAAACTGATCACAACTCACAAATCTTGATATACATaa                                                                                   >  2:1567242/1‑91 (MQ=255)
          aaaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAAt                                                                                  >  1:2646974/1‑91 (MQ=255)
           aaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAAtt                                                                                 >  1:171912/1‑91 (MQ=255)
             tGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtg                                                                               <  1:1690494/91‑1 (MQ=255)
             ggATGACGAAATGGATTCTCTTAATTGAAGTAATAGGGAACCCAAATTCATTAAATGAGCACAAATCACAAATTATGTATATCATAATTtg                                                                               <  2:310274/90‑1 (MQ=12)
               aTTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgtg                                                                             >  2:3267883/1‑91 (MQ=255)
               aTTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgtg                                                                             <  1:1177130/91‑1 (MQ=255)
                ttACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGa                                                                            <  2:3853140/91‑1 (MQ=255)
                ttACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGa                                                                            >  1:2752064/1‑91 (MQ=255)
                ttACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGa                                                                            >  1:1081252/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)
                   cGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAttt                                                                         <  1:3606987/91‑1 (MQ=255)
                    gAATGGGATTCACTTAATTGAAGGAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTGGATATACATAATTTGTGAtttt                                                                        <  2:568538/91‑1 (MQ=255)
                      aTGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCCCAAATTTTGATATACATAATTTGTGAtttttt                                                                      <  2:3884746/91‑1 (MQ=255)
                       tGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTta                                                                     >  1:1575646/1‑91 (MQ=255)
                       tGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTta                                                                     <  2:1073117/91‑1 (MQ=255)
                           aTTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATAtt                                                                 <  1:2493110/91‑1 (MQ=255)
                            ttCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTc                                                                <  1:825794/91‑1 (MQ=255)
                            ttCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTc                                                                >  1:2945141/1‑91 (MQ=255)
                              cACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCaa                                                              <  2:3082750/91‑1 (MQ=255)
                                  tAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCt                                                          <  2:918149/91‑1 (MQ=255)
                                    aTTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGGGATTTTTTATATTCAAAGCTaa                                                        <  2:696172/91‑1 (MQ=255)
                                     ttGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTaaa                                                       >  1:2443182/1‑91 (MQ=255)
                                        aaGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAtt                                                    >  2:476002/1‑91 (MQ=255)
                                         aGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTg                                                   >  2:3403880/1‑91 (MQ=255)
                                           tAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCa                                                 <  2:1021279/91‑1 (MQ=255)
                                            aaaaCGTTACCCAAGTTCAAAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaa                                                <  2:1984053/91‑1 (MQ=255)
                                                cGTTAGCCAAGTTCATAAAATGATCACAAATCACGACTTTTGATATACGTAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaaaaa                                            <  2:224176/91‑1 (MQ=255)
                                                agTTACCCAAGTTCATAAAAAGATCAAAAATCACAAATTTTGATATACATAATTTGTAATATTTAATAATAAAAGAAAAAATAaaaaaaaa                                            <  2:2049151/90‑1 (MQ=255)
                                                 gTTACCCAAGTTAATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAt                                           >  1:3048491/1‑91 (MQ=255)
                                                  ttACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAtt                                          <  2:1334061/91‑1 (MQ=255)
                                                  ttACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATATGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAtt                                          <  2:1146926/91‑1 (MQ=255)
                                                   tACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTa                                         >  1:681849/1‑91 (MQ=255)
                                                    aCCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                        >  2:456809/1‑91 (MQ=255)
                                                    aCCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                        >  1:272983/1‑91 (MQ=255)
                                                    aCCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                        <  1:3436090/91‑1 (MQ=255)
                                                    aCCCAAGTAAATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                        <  2:1234667/91‑1 (MQ=255)
                                                      ccAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATg                                      >  2:2585775/1‑91 (MQ=255)
                                                       cAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTATAATTGCAAAAAATTAATgg                                     <  2:3461263/91‑1 (MQ=255)
                                                         aGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGtt                                   >  1:2720526/1‑91 (MQ=255)
                                                         aGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGtt                                   <  1:3211623/91‑1 (MQ=255)
                                                          gTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTa                                  <  2:2218184/91‑1 (MQ=255)
                                                            tCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAAc                                >  2:2546658/1‑91 (MQ=255)
                                                               tAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtc                             >  1:1061980/1‑91 (MQ=255)
                                                               tAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtc                             >  2:75741/1‑91 (MQ=255)
                                                               tAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtc                             >  1:3124568/1‑91 (MQ=255)
                                                                aaaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtct                            <  1:957417/91‑1 (MQ=255)
                                                                aaaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtct                            <  2:2350910/91‑1 (MQ=255)
                                                                aaaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtct                            <  2:3701227/91‑1 (MQ=255)
                                                                 aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                           <  1:564515/91‑1 (MQ=255)
                                                                 aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                           <  1:2865283/91‑1 (MQ=255)
                                                                 aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                           >  1:2851765/1‑91 (MQ=255)
                                                                  aaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATctct                          <  1:936213/91‑1 (MQ=255)
                                                                  aaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATctct                          <  1:1875445/91‑1 (MQ=255)
                                                                  aaTGATCACAAATCACAAAGTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATctct                          <  1:3903004/91‑1 (MQ=255)
                                                                       tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAGAATTAATTGTTAACATCTCcgttgt                     >  2:1257186/1‑91 (MQ=255)
                                                                       tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                     >  2:3135198/1‑91 (MQ=255)
                                                                       tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                     <  2:3061715/91‑1 (MQ=255)
                                                                       tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                     >  1:2020922/1‑91 (MQ=255)
                                                                       tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                     >  1:1626305/1‑91 (MQ=255)
                                                                       tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                     <  2:2278904/91‑1 (MQ=255)
                                                                       taaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                     >  2:2368792/5‑91 (MQ=21)
                                                                        aaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                    >  2:2926623/4‑91 (MQ=21)
                                                                         aaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGttgttg                   >  1:3346561/3‑91 (MQ=21)
                                                                          caaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTgg                  <  1:1761096/91‑1 (MQ=255)
                                                                           aaatcacaaatTTTGATATACATAATTTGTGGTTTTTTATATTCTAAGCTAAAGTTGGAAAAAATTACTGGTTAACATCTCTGTTGTTGGc                 <  1:3155341/91‑1 (MQ=255)
                                                                           aaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGc                 >  2:1593482/1‑91 (MQ=255)
                                                                            aatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCa                >  2:1814365/1‑91 (MQ=25)
                                                                              tcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAat              <  1:1532003/91‑1 (MQ=255)
                                                                              tcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGAAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAat              <  2:2436113/91‑1 (MQ=255)
                                                                              tcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAATAATTAATGGTTAACATCTCTGTTGTTGGCAat              >  2:1275382/1‑91 (MQ=25)
                                                                                acaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAatat            >  1:1027658/1‑91 (MQ=255)
                                                                                acaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAatat            >  2:360846/1‑91 (MQ=255)
                                                                                acaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAGTAATGGTTAACATCTCTGTTGTTGGCAatat            <  1:3533177/91‑1 (MQ=255)
                                                                                 caaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAAtata           <  1:2369162/91‑1 (MQ=255)
                                                                                    atTTTGATATACATAATTTGTGATTTTTTTTTTTTTTTGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaat        >  2:3735445/1‑91 (MQ=255)
                                                                                    atTTTGATATACATAATTTGTGATTTTTTTTATTTTAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaat        >  2:2962862/1‑91 (MQ=255)
                                                                                    atTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaat        >  1:2715978/1‑91 (MQ=255)
                                                                                    atTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaat        <  1:1381099/91‑1 (MQ=255)
                                                                                       ttGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaatgat     <  1:1619652/91‑4 (MQ=255)
                                                                                        tGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATg    <  2:3945670/91‑1 (MQ=255)
                                                                                        tGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAACTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATg    <  2:3469105/91‑1 (MQ=255)
                                                                                         gATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGa   <  2:1776994/91‑1 (MQ=255)
                                                                                         gATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGa   >  1:3847680/1‑91 (MQ=255)
                                                                                         gATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGa   <  1:20111/91‑1 (MQ=255)
                                                                                         gATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTCGCAATATAAATAATGa   >  2:3789278/1‑91 (MQ=255)
                                                                                          atatTCATAATTTGTGATTTTTTATATGCAATGCTAATATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGAt  <  2:2135103/91‑1 (MQ=255)
                                                                                          atatACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGAt  <  2:2467319/91‑1 (MQ=255)
                                                                                          |                                                                                           
TTAAACATTAAAATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGAT  >  USA300TCH1516_ALE/115668‑115848

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