Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,733,446 T→A 5.7% intergenic (‑27/+73) 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,4460TA 5.7% 565.4 / 14.7 169intergenic (‑27/+73)csbD_2/cymRStress response protein CsbD/HTH‑type transcriptional regulator CymR
Reads supporting (aligned to +/- strand):  ref base T (82/76);  new base A (3/7);  total (85/84)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

TAACGTTACCTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAATATAAAACATGTA  >  USA300TCH1516_ALE/1733356‑1733536
                                                                                          |                                                                                          
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:656217/1‑91 (MQ=255)
       aCCTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATAc                                                                                     <  2:1353662/91‑1 (MQ=255)
       aCCTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATAc                                                                                     >  2:733176/1‑91 (MQ=255)
       aCCTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATAc                                                                                     >  1:654569/1‑91 (MQ=255)
        ccTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATAcc                                                                                    >  1:1089687/1‑91 (MQ=255)
        ccTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATAAATATATAca                                                                                    <  2:3884984/91‑2 (MQ=255)
        ccTACTGTTTCTTTAACATTTCCTTTAAATTCATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATATATATATTTATATAcc                                                                                    >  2:2971092/1‑91 (MQ=255)
         cTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCa                                                                                   >  2:2178409/1‑91 (MQ=255)
          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)
            cTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTa                                                                                <  2:1321972/91‑1 (MQ=255)
               tttCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAagg                                                                             >  1:2121423/1‑91 (MQ=255)
               tttCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAagg                                                                             <  2:2019972/91‑1 (MQ=255)
                ttCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAagga                                                                            >  1:3045847/1‑91 (MQ=255)
                ttCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAagga                                                                            >  1:1155506/1‑91 (MQ=255)
                  cTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAaggagg                                                                          <  1:817820/91‑1 (MQ=255)
                  cTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAaggagg                                                                          <  1:105768/91‑1 (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)
                      aaCATTTCGTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCg                                                                      <  2:3375593/91‑1 (MQ=255)
                           ttCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTaaa                                                                 <  2:1255259/91‑1 (MQ=255)
                           ttCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTaaa                                                                 <  2:3818825/91‑1 (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                                                           <  2:2847324/91‑1 (MQ=255)
                                 tAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATtat                                                           <  1:319222/91‑1 (MQ=255)
                                  aaaTTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATtata                                                          >  2:1043640/1‑91 (MQ=255)
                                  aaaTTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATtata                                                          >  2:2465265/1‑91 (MQ=255)
                                    aTTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTataat                                                        >  1:2012514/1‑91 (MQ=255)
                                      tGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAatat                                                      <  2:3335905/91‑1 (MQ=255)
                                      tGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAatat                                                      >  1:133204/1‑91 (MQ=255)
                                      tGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAatat                                                      <  2:2656095/91‑1 (MQ=255)
                                      tGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAatat                                                      <  1:629386/91‑1 (MQ=255)
                                        aTCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGa                                                    <  2:2165209/91‑1 (MQ=255)
                                              cttacttTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGtttt                                              <  1:852633/91‑1 (MQ=255)
                                              cttacttTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGtttt                                              <  1:2291826/91‑1 (MQ=255)
                                                tacttTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTa                                            <  1:1907129/91‑1 (MQ=255)
                                                   tttCGTCTGCCACAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGt                                         <  2:274679/91‑1 (MQ=255)
                                                    ttCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGtt                                        >  2:3549763/1‑91 (MQ=255)
                                                     tCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTta                                       >  2:1667834/1‑91 (MQ=255)
                                                     tCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTta                                       >  1:560849/1‑91 (MQ=255)
                                                      cGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTtat                                      >  1:588567/1‑91 (MQ=255)
                                                      cGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTtat                                      >  1:2027508/1‑91 (MQ=255)
                                                       gTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTtata                                     <  1:1395459/91‑1 (MQ=255)
                                                        tCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTtatag                                    <  1:2541674/91‑1 (MQ=255)
                                                        tCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTtatag                                    >  1:785657/1‑91 (MQ=255)
                                                         cTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATAAAGTTTTTATGTtatagt                                   <  2:1201560/91‑1 (MQ=255)
                                                          tGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAGGTTTTTATGTtatagta                                  <  2:906689/91‑1 (MQ=255)
                                                             cATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTtatagtatag                               >  2:337540/1‑91 (MQ=255)
                                                             cATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTtatagtatag                               <  1:691486/91‑1 (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:2551210/1‑91 (MQ=255)
                                                                     aaCCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTtata                       >  2:1790730/1‑91 (MQ=255)
                                                                     aaCCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTtata                       >  2:1790414/1‑91 (MQ=255)
                                                                     aaCCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTtata                       >  2:689798/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)
                                                                           ctTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTa                 <  2:551381/91‑1 (MQ=255)
                                                                           ctTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTa                 >  2:453030/1‑91 (MQ=255)
                                                                           ctTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTa                 <  2:551780/91‑1 (MQ=255)
                                                                           ctTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTa                 >  2:2941892/1‑91 (MQ=255)
                                                                           ctTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTa                 <  2:2404277/91‑1 (MQ=255)
                                                                           ctTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTa                 >  1:1255787/1‑91 (MQ=255)
                                                                            ttGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTaa                >  2:261011/1‑91 (MQ=255)
                                                                             tGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTaaa               >  2:3269800/1‑91 (MQ=255)
                                                                             tGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTaaa               >  1:2634053/1‑91 (MQ=255)
                                                                              ggATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAat              <  2:3272188/91‑1 (MQ=255)
                                                                              ggATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAat              <  2:3341890/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)
                                                                               gATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAata             >  2:3097655/1‑91 (MQ=255)
                                                                                  gTATATAAATATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAATATaa          >  1:37036/1‑91 (MQ=255)
                                                                                    atatatTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAATATaaaa        <  2:588567/91‑1 (MQ=255)
                                                                                    atatatTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAATATaaaa        <  2:1427562/91‑1 (MQ=255)
                                                                                     tatatTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGCTTAAATATAAAAc       >  1:1669895/1‑91 (MQ=255)
                                                                                      atatTTATATACCACTAAGGAGGTCCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAATATAAAACa      <  2:3613386/91‑1 (MQ=255)
                                                                                       tatTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAATATAAAACAt     >  2:1142848/1‑91 (MQ=255)
                                                                                          |                                                                                          
TAACGTTACCTACTGTTTCTTTAACATTTCCTTTAAATTGATCGAACTTACTTTCGTCTGCCATAGTGAAACCTCCTTGGATGTATATATTTATATACCACTAAGGAGGTTCGCTAAACATTATAATATGAAGTTTTTATGTTATAGTATAGTATTTATACGATTAAATATAAAACATGTA  >  USA300TCH1516_ALE/1733356‑1733536

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: