Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 115,755 T→A 13.2% intergenic (+22/‑129) 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,7550TA13.2% 622.9 / 75.0 214intergenic (+22/‑129)btr/yxeP_1HTH‑type transcriptional activator Btr/putative hydrolase YxeP
Reads supporting (aligned to +/- strand):  ref base T (90/94);  new base A (22/6);  total (113/101)
Fisher's exact test for biased strand distribution p-value = 4.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.75e-01

Only 100 of 220 total aligned reads displayed.

CATTAAACATTAAAATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAA  >  USA300TCH1516_ALE/115666‑115844
                                                                                         |                                                                                          
cattaaacattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAAtttt                                                                                           >  2:2990040/1‑91 (MQ=255)
  ttaaacattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGa                                                                                         <  1:1044155/91‑1 (MQ=255)
    aaacattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGata                                                                                       >  1:2184353/1‑91 (MQ=255)
     aacattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGatat                                                                                      >  1:28660/1‑91 (MQ=255)
      acattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGAtata                                                                                     <  2:3119664/91‑1 (MQ=255)
       cattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATAc                                                                                   <  2:1789536/91‑1 (MQ=255)
        attaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACa                                                                                  >  2:1750198/1‑91 (MQ=255)
         ttaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACAt                                                                                 <  1:1926522/91‑1 (MQ=255)
          taaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATa                                                                                >  2:3308352/1‑91 (MQ=255)
          taaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATa                                                                                >  1:1565641/1‑91 (MQ=255)
          taaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATa                                                                                <  2:1386297/91‑1 (MQ=255)
              aTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAAttt                                                                            <  2:1442723/91‑1 (MQ=255)
              aTGATTACGAATGGGATTAACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAAttt                                                                            >  2:3063947/1‑91 (MQ=255)
               tGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtg                                                                           <  2:684201/91‑1 (MQ=255)
               tGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtg                                                                           <  1:3659418/91‑1 (MQ=255)
                 aTTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTATGTATATCAAAATTtgtg                                                                         <  1:1082539/91‑1 (MQ=37)
                   tACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAt                                                                       >  1:2886003/1‑91 (MQ=255)
                    aCGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtt                                                                      <  1:1205158/91‑1 (MQ=255)
                     cGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAttt                                                                     >  1:1338398/1‑91 (MQ=255)
                     cGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAttt                                                                     <  1:2757649/91‑1 (MQ=255)
                     cGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTATGTATATCAAAATTTGTGAttt                                                                     <  1:1762629/91‑1 (MQ=25)
                      gAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtttt                                                                    >  1:717251/1‑91 (MQ=255)
                      gAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtttt                                                                    >  2:482781/1‑91 (MQ=255)
                      gAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtttt                                                                    >  1:164839/1‑91 (MQ=255)
                      gAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTATGTATATCAAAATTTGTGAtttg                                                                    <  1:2266492/91‑2 (MQ=25)
                       aaTGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAttttt                                                                   >  2:2073376/1‑91 (MQ=255)
                        aTGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtttttt                                                                  <  2:414086/91‑1 (MQ=255)
                        aTGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTATGTATATCAAAATTTGTGAtttgtg                                                                  <  1:2940193/91‑4 (MQ=21)
                         tGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTta                                                                 >  2:3519451/1‑91 (MQ=255)
                          gggATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTtat                                                                >  2:3835650/1‑91 (MQ=255)
                          gggATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTtat                                                                <  2:2943428/91‑1 (MQ=255)
                           ggATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTtata                                                               >  1:5404/1‑91 (MQ=255)
                            gATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTatat                                                              >  1:2013145/1‑91 (MQ=255)
                             aTTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATAtt                                                             <  1:3741724/91‑1 (MQ=255)
                               tCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCa                                                           >  2:3408626/1‑91 (MQ=255)
                               tCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCa                                                           >  1:2406748/1‑91 (MQ=255)
                                cACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCaa                                                          >  1:285198/1‑91 (MQ=255)
                                  cTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAg                                                        <  2:2810558/91‑1 (MQ=255)
                                   ttAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATACACATAATTTGTGATTTTTTGTATTCAAAgt                                                       <  2:3809001/91‑2 (MQ=255)
                                    tAATTGAAGTAATACGTTACCCAAGTTCATAAAATGCTCACAAATCACATATTGTGATATACATAATTTGTGATTTTTTATATTCAAAGCt                                                      <  2:2053375/91‑1 (MQ=255)
                                     aaTTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTa                                                     >  1:752812/1‑91 (MQ=255)
                                         gAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAt                                                 <  2:3635145/91‑1 (MQ=255)
                                         gAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAt                                                 <  1:1061508/91‑1 (MQ=255)
                                          aaGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAtt                                                >  2:1878156/1‑91 (MQ=255)
                                          aaGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAtt                                                >  1:641914/1‑91 (MQ=255)
                                          aaGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAtt                                                <  2:1697878/91‑1 (MQ=255)
                                              aaaaCGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaa                                            >  1:76822/1‑91 (MQ=255)
                                               aaaCGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaa                                           >  1:845764/1‑91 (MQ=255)
                                                 aCGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaaaa                                         <  2:394274/91‑1 (MQ=255)
                                                  cGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaaaaa                                        <  2:792544/91‑1 (MQ=255)
                                                   gTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAt                                       >  2:2595678/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)
                                                    ttACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAtt                                      <  1:3375875/91‑1 (MQ=255)
                                                    ttACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAtt                                      <  2:1720244/91‑1 (MQ=255)
                                                    ttACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATAGACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAtt                                      <  1:3881515/91‑1 (MQ=255)
                                                     tACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTa                                     <  2:1641924/91‑1 (MQ=255)
                                                      aCCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                    >  1:764283/1‑91 (MQ=255)
                                                        ccAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATg                                  >  1:813746/1‑91 (MQ=255)
                                                         cAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATgg                                 <  2:1969698/91‑1 (MQ=255)
                                                          aaGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGt                                >  2:2652217/1‑91 (MQ=255)
                                                            gTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTa                              >  1:3101092/1‑91 (MQ=255)
                                                            gTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTa                              <  2:2184353/91‑1 (MQ=255)
                                                             ttCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTaa                             >  1:665185/1‑91 (MQ=255)
                                                             ttCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTaa                             >  2:2009034/1‑91 (MQ=255)
                                                                aTAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAt                          >  1:2667940/1‑91 (MQ=255)
                                                                 tAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtc                         >  1:3770566/1‑91 (MQ=255)
                                                                   aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                       <  2:28661/91‑1 (MQ=255)
                                                                   aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                       >  1:2381862/1‑91 (MQ=255)
                                                                   aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                       <  2:540903/91‑1 (MQ=255)
                                                                   aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                       >  1:3433707/1‑91 (MQ=255)
                                                                   aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                       <  1:20336/91‑1 (MQ=255)
                                                                    aaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATctct                      <  2:1616330/91‑1 (MQ=255)
                                                                    aaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAAAGAAAAAAATTAATGGTTAACATctct                      <  2:2372992/91‑1 (MQ=255)
                                                                     aTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtg                     >  1:252784/1‑91 (MQ=255)
                                                                      tGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTGAAAGCTAAAATAGAAAAAAATTAATGGTTAACATCTCtgt                    <  2:3188135/91‑1 (MQ=255)
                                                                        atcacaaatcacaaATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttg                  <  1:2255277/91‑1 (MQ=255)
                                                                         taaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                 >  2:1963868/5‑91 (MQ=21)
                                                                         taaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                 >  2:2889556/5‑91 (MQ=21)
                                                                         taaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                 >  1:876616/5‑91 (MQ=21)
                                                                         taaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                 >  1:1843626/5‑91 (MQ=21)
                                                                          cacaaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                <  2:1145873/91‑1 (MQ=255)
                                                                          cacaaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                >  2:622880/1‑91 (MQ=255)
                                                                          aaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                >  1:3800803/4‑91 (MQ=21)
                                                                          aaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                >  2:3717809/4‑91 (MQ=21)
                                                                          aaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                >  2:132710/4‑91 (MQ=21)
                                                                           acaaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGttgttg               >  1:2600937/1‑91 (MQ=255)
                                                                           acaaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGttgttg               <  1:2509622/91‑1 (MQ=255)
                                                                           acaaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGttgttg               >  1:1728174/1‑91 (MQ=255)
                                                                            caaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTgg              >  2:2366776/1‑91 (MQ=255)
                                                                              aatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCa            >  2:52322/1‑91 (MQ=255)
                                                                               atcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCaa           <  2:1427724/91‑1 (MQ=255)
                                                                                tcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAat          >  1:1252872/1‑91 (MQ=37)
                                                                                tcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAat          >  2:3340353/1‑91 (MQ=37)
                                                                                   caaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAAtata       <  1:3497972/91‑1 (MQ=255)
                                                                                   caaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAAtata       <  2:1518161/91‑1 (MQ=255)
                                                                                   caaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATAGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAAtata       <  2:1302790/91‑1 (MQ=255)
                                                                                    aaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataa      >  2:1405498/1‑91 (MQ=255)
                                                                                     aatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaa     >  1:1750826/1‑91 (MQ=255)
                                                                                         |                                                                                          
CATTAAACATTAAAATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAA  >  USA300TCH1516_ALE/115666‑115844

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