Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,733,445 T→A 5.6% intergenic (‑26/+74) csbD_2 ← / ← cymR Stress response protein CsbD/HTH‑type transcriptional regulator CymR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,733,4450TA 5.6% 549.1 / 15.9 162intergenic (‑26/+74)csbD_2/cymRStress response protein CsbD/HTH‑type transcriptional regulator CymR
Reads supporting (aligned to +/- strand):  ref base T (78/74);  new base A (3/6);  total (82/80)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 175 total aligned reads displayed.

GTAACGTTACCTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAATATAAAACAT  >  USA300TCH1516_ALE/1733355‑1733533
                                                                                          |                                                                                        
 tAACGTTACCTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATAtt                                                                                         <  2:1843739/91‑1 (MQ=255)
  aaCGTTACCTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATAttt                                                                                        >  2:1620507/1‑91 (MQ=255)
      ttACCTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTCCTTTCGTCTGCCATCGTGAAACCTCCTTGGACGTATATACTtctat                                                                                    >  2:2351201/1‑91 (MQ=37)
       tACCTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTtatata                                                                                   >  1:2140417/1‑91 (MQ=255)
       tACCTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTtatata                                                                                   <  1:3299279/91‑1 (MQ=255)
        aCCTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATAc                                                                                  >  1:654569/1‑91 (MQ=255)
        aCCTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATAc                                                                                  <  2:1353662/91‑1 (MQ=255)
        aCCTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATAc                                                                                  >  1:656217/1‑91 (MQ=255)
        aCCTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATAc                                                                                  >  2:733176/1‑91 (MQ=255)
         ccTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATAcc                                                                                 >  1:1089687/1‑91 (MQ=255)
          cTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCa                                                                                >  2:2178409/1‑91 (MQ=255)
           tACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTAGTGGTATATAAATATATAcatc                                                                               <  2:3508315/91‑4 (MQ=25)
           tACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTAGTGGTATATAAATATATAcatc                                                                               <  2:1565792/91‑4 (MQ=25)
            aCTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATTGTGAAACCTCCTTGGATGTATATATTTATATACCACt                                                                              <  1:3280168/91‑1 (MQ=255)
            aCTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATCGTGAAACCTCCTTGGATGTATATATTTATATACCACt                                                                              <  1:3280249/91‑1 (MQ=255)
            aCTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACt                                                                              <  1:2744785/91‑1 (MQ=255)
            aCTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACt                                                                              <  2:3302258/91‑1 (MQ=255)
            aCTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACt                                                                              <  2:408060/91‑1 (MQ=255)
            aCTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACt                                                                              >  2:557059/1‑91 (MQ=255)
             cTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTa                                                                             <  2:1321972/91‑1 (MQ=255)
                tttCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAagg                                                                          >  1:2121423/1‑91 (MQ=255)
                tttCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAagg                                                                          <  2:2019972/91‑1 (MQ=255)
                 ttCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAagga                                                                         >  1:1155506/1‑91 (MQ=255)
                 ttCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAagga                                                                         >  1:3045847/1‑91 (MQ=255)
                   cTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAaggagg                                                                       <  1:105768/91‑1 (MQ=255)
                   cTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAaggagg                                                                       <  1:817820/91‑1 (MQ=255)
                    tttAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGt                                                                      <  2:33546/91‑1 (MQ=255)
                    tttAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGt                                                                      >  2:3563791/1‑91 (MQ=255)
                    tttAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTAGTGGTATATAAATATATACATCCAAGGAGGt                                                                      <  1:1940607/91‑1 (MQ=17)
                     ttAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTAAATATTCATATACCACTAaggacgat                                                                     >  1:171424/1‑87 (MQ=255)
                     ttAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTAGTGGTATATAAATATATACATCCAAGGAGGtt                                                                     <  1:1695190/91‑1 (MQ=17)
                      tAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTc                                                                    <  2:2452800/91‑1 (MQ=255)
                        aCATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGc                                                                  <  2:2725765/91‑1 (MQ=255)
                            ttCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTaaa                                                              <  2:1255259/91‑1 (MQ=255)
                             tCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAAc                                                             >  2:2956465/1‑91 (MQ=255)
                               cTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACAt                                                           >  2:95575/1‑91 (MQ=255)
                                tttAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACAtt                                                          >  1:3779404/1‑91 (MQ=255)
                                 ttAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATta                                                         >  2:2270722/1‑91 (MQ=255)
                                 ttAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATta                                                         >  2:789364/1‑91 (MQ=255)
                                 ttAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATGTATATACCACTAAGGAGGTTCGCTAAACATta                                                         <  2:3809812/91‑1 (MQ=255)
                                  tAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATtat                                                        <  1:319222/91‑1 (MQ=255)
                                  tAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATtat                                                        <  2:2847324/91‑1 (MQ=255)
                                   aaaTTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATtata                                                       >  2:2465265/1‑91 (MQ=255)
                                   aaaTTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATtata                                                       >  2:1043640/1‑91 (MQ=255)
                                     aTTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTataat                                                     >  1:2012514/1‑91 (MQ=255)
                                       tGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAatat                                                   <  1:629386/91‑1 (MQ=255)
                                       tGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAatat                                                   <  2:2656095/91‑1 (MQ=255)
                                       tGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAatat                                                   >  1:133204/1‑91 (MQ=255)
                                         aTCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGa                                                 <  2:2165209/91‑1 (MQ=255)
                                              acttacttTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGttt                                            >  2:3793047/1‑91 (MQ=255)
                                               cttacttTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGtttt                                           <  1:852633/91‑1 (MQ=255)
                                               cttacttTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGtttt                                           <  1:2291826/91‑1 (MQ=255)
                                                 tacttTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTa                                         <  1:1907129/91‑1 (MQ=255)
                                                      tCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTta                                    >  1:560849/1‑91 (MQ=255)
                                                      tCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTta                                    >  2:1667834/1‑91 (MQ=255)
                                                       cGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTtat                                   >  1:2027508/1‑91 (MQ=255)
                                                       cGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTtat                                   >  1:588567/1‑91 (MQ=255)
                                                        gTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTtata                                  <  1:1395459/91‑1 (MQ=255)
                                                         tCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTtatag                                 >  1:785657/1‑91 (MQ=255)
                                                         tCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTtatag                                 <  1:2541674/91‑1 (MQ=255)
                                                          cTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATAAAGTTTTTATGTtatagt                                <  2:1201560/91‑1 (MQ=255)
                                                           tGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAGGTTTTTATGTtatagta                               <  2:906689/91‑1 (MQ=255)
                                                              cATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTtatagtatag                            <  1:691486/91‑1 (MQ=255)
                                                              cATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTtatagtatag                            >  2:337540/1‑91 (MQ=255)
                                                                  gTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTAtt                        >  1:3772885/1‑91 (MQ=255)
                                                                   tGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTAttt                       <  2:1482801/91‑1 (MQ=255)
                                                                   tGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTAttt                       <  1:2889213/91‑1 (MQ=255)
                                                                    gAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTta                      <  1:1281247/91‑1 (MQ=255)
                                                                    gAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTta                      <  1:845685/91‑1 (MQ=255)
                                                                     aaaCCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTtat                     >  1:2421762/1‑91 (MQ=255)
                                                                      aaCCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAGACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTtata                    >  1:3871965/1‑91 (MQ=255)
                                                                      aaCCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTtata                    >  2:689798/1‑91 (MQ=255)
                                                                      aaCCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTtata                    >  2:2551210/1‑91 (MQ=255)
                                                                      aaCCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTtata                    >  2:1790730/1‑91 (MQ=255)
                                                                      aaCCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTtata                    >  2:1790414/1‑91 (MQ=255)
                                                                       aCCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATAc                   >  1:2160886/1‑91 (MQ=255)
                                                                        cctcctTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACg                  >  1:959402/1‑91 (MQ=255)
                                                                          tcctTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGAt                >  1:3115776/1‑91 (MQ=255)
                                                                          tcctTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGAt                >  1:640857/1‑91 (MQ=255)
                                                                          tcctTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGAt                >  2:3001034/1‑91 (MQ=255)
                                                                            ctTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTa              <  2:470668/91‑1 (MQ=255)
                                                                            ctTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTa              >  1:1255787/1‑91 (MQ=255)
                                                                            ctTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTa              <  2:551381/91‑1 (MQ=255)
                                                                            ctTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTa              <  2:2404277/91‑1 (MQ=255)
                                                                            ctTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATCAAGTTTTTATGTTATAGTATAGTATTTATACGATTa              >  2:3831992/1‑91 (MQ=255)
                                                                             ttGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTaa             >  2:261011/1‑91 (MQ=255)
                                                                              tGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTaaa            >  2:3269800/1‑91 (MQ=255)
                                                                              tGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTaaa            >  2:3170127/1‑91 (MQ=255)
                                                                              tGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTaaa            >  1:2634053/1‑91 (MQ=255)
                                                                               ggATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAat           <  2:3341778/91‑1 (MQ=255)
                                                                               ggATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAat           <  2:2566704/91‑1 (MQ=255)
                                                                               ggATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAat           <  2:2108833/91‑1 (MQ=255)
                                                                               agtggTATATAAATATATACATCCAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAat           >  1:487276/5‑91 (MQ=25)
                                                                                gtggTATATAAATATATACATCCAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAata          >  1:1009232/4‑91 (MQ=25)
                                                                                gATGTATATCTTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAata          >  1:1951374/1‑91 (MQ=255)
                                                                                   gTATATAAATATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAATATaa       >  1:37036/1‑91 (MQ=255)
                                                                                    tatataTTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAATATaaa      >  2:2778483/1‑91 (MQ=255)
                                                                                     atatatTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAATATaaaa     <  2:1427562/91‑1 (MQ=255)
                                                                                      tatatTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGCTTAAATATAAAAc    >  1:1669895/1‑91 (MQ=255)
                                                                                        tatTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAATATAAAACAt  >  2:1142848/1‑91 (MQ=255)
                                                                                          |                                                                                        
GTAACGTTACCTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAATATAAAACAT  >  USA300TCH1516_ALE/1733355‑1733533

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: