Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,833,387 G→C 10.3% intergenic (‑723/+367) lipA_2 ← / ← hisI Lipase 1/Phosphoribosyl‑AMP cyclohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,833,3870GC10.3% 702.3 / 56.2 214intergenic (‑723/+367)lipA_2/hisILipase 1/Phosphoribosyl‑AMP cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base G (94/98);  new base C (7/15);  total (101/113)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 217 total aligned reads displayed.

TTCGGGATGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCATTTTTTCCCTT  >  USA300TCH1516_ALE/2833298‑2833477
                                                                                         |                                                                                          
 tCGGGATGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGta                                                                                          <  2:2309513/91‑1 (MQ=255)
  cGGGATGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGtat                                                                                         >  1:154034/1‑91 (MQ=255)
  cGGGATGCATATCTGAGTTCTAAGGGAGAGTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATAGATTTGTATTTATACtat                                                                                         <  2:1279766/91‑1 (MQ=255)
  cGGGATACATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGtat                                                                                         >  1:3018850/1‑91 (MQ=255)
   gggATGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGtata                                                                                        >  1:512825/1‑91 (MQ=255)
    ggATGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATaa                                                                                       >  1:372049/1‑91 (MQ=255)
    ggATGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATaa                                                                                       >  1:94845/1‑91 (MQ=255)
    ggATGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATaa                                                                                       >  2:1814387/1‑91 (MQ=255)
     gATGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATaaa                                                                                      <  1:409554/91‑1 (MQ=255)
     gATGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATaaa                                                                                      <  1:2710396/91‑1 (MQ=255)
      aTGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTCCATATATTTGCATTTATAGTATAAAt                                                                                     >  2:3364933/1‑91 (MQ=255)
      aTGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAAt                                                                                     <  2:813352/91‑1 (MQ=255)
      aTGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAAt                                                                                     >  2:3810880/1‑91 (MQ=255)
      aTGCAGATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAAt                                                                                     <  2:3574677/91‑1 (MQ=255)
        gCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATAc                                                                                   <  1:3585223/91‑1 (MQ=255)
        gCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATACTATAAATAc                                                                                   <  1:1105097/91‑1 (MQ=255)
         cATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATACTATAAATACa                                                                                  <  1:2012216/91‑1 (MQ=255)
          atatCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATACTATAAATACaa                                                                                 <  1:3181123/91‑1 (MQ=255)
           tatCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACaaa                                                                                <  2:1598809/91‑1 (MQ=255)
           tatCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACaaa                                                                                >  1:2426532/1‑91 (MQ=255)
           tatCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACaaa                                                                                <  2:536005/91‑1 (MQ=255)
           tatCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATACTATAAATACaaa                                                                                <  2:752220/91‑1 (MQ=255)
            atCTTGGTTCTGATGGAGGTTGTCTTTTAAATAGGCTATGACATTACTATTTCAGGCTAGATATATTTGTAGTTATACTATAAATACaaat                                                                               <  2:3857729/91‑2 (MQ=255)
                tAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATg                                                                           <  1:1425758/91‑1 (MQ=255)
                  gTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTa                                                                         <  1:755162/91‑1 (MQ=255)
                  gTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTa                                                                         >  2:651014/1‑91 (MQ=255)
                   ttCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTaa                                                                        <  2:875597/91‑1 (MQ=255)
                     cTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTaaaa                                                                      >  2:3848000/1‑90 (MQ=255)
                     cTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAg                                                                      <  2:2969796/91‑1 (MQ=255)
                     cTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAg                                                                      <  1:764952/91‑1 (MQ=255)
                        aTGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACa                                                                   >  2:986556/1‑91 (MQ=255)
                           gagaTTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACt                                                                >  1:1222247/1‑91 (MQ=255)
                           gagaTTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACt                                                                >  1:993526/1‑91 (MQ=255)
                           gagaTTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACt                                                                <  1:1264871/91‑1 (MQ=255)
                              aTTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATg                                                             >  1:461418/1‑91 (MQ=255)
                              aTTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATg                                                             <  2:2701928/91‑1 (MQ=255)
                              aTTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATg                                                             >  1:238890/1‑91 (MQ=255)
                               ttGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGa                                                            >  1:1924251/1‑91 (MQ=255)
                                tGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGat                                                           <  2:1319745/91‑1 (MQ=255)
                                 gTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGata                                                          >  1:2019292/1‑91 (MQ=255)
                                  tCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGatat                                                         >  2:60427/1‑91 (MQ=255)
                                  tCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGatat                                                         >  2:392569/1‑91 (MQ=255)
                                   cTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATAtt                                                        <  2:238890/91‑1 (MQ=255)
                                    ttttAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTa                                                       >  2:3508210/1‑91 (MQ=255)
                                      ttAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTGTAAATACAAAGAATGTAAAGACAACTATGATATTAAt                                                     <  2:1479754/91‑1 (MQ=255)
                                         aaTAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAgg                                                  <  2:438936/91‑1 (MQ=255)
                                         aaTAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAgg                                                  <  2:2284501/91‑1 (MQ=255)
                                              gCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGgtagta                                             <  1:3473241/91‑1 (MQ=255)
                                              gCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGgtagta                                             <  1:3771246/91‑1 (MQ=255)
                                               cTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTaa                                            >  1:1710213/1‑91 (MQ=255)
                                                  tGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGa                                         <  1:1577416/91‑1 (MQ=255)
                                                  tGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGa                                         >  1:194685/1‑91 (MQ=255)
                                                   gACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGAc                                        >  2:91448/1‑91 (MQ=255)
                                                    aCATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACa                                       <  2:1797146/91‑1 (MQ=255)
                                                    aCATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACa                                       <  2:3747689/91‑1 (MQ=255)
                                                     cATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACaa                                      >  1:2239284/1‑91 (MQ=255)
                                                         tCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACAt                                  <  2:1231885/91‑1 (MQ=255)
                                                           tATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATaa                                >  1:2440039/1‑91 (MQ=255)
                                                             tttCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAAct                              >  1:3042307/1‑91 (MQ=255)
                                                              ttCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAActc                             <  2:1935993/91‑1 (MQ=255)
                                                              ttCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAActc                             <  2:1936178/91‑1 (MQ=255)
                                                              ttCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAActc                             <  2:3771105/91‑1 (MQ=255)
                                                                cAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCtt                           <  1:3145488/91‑1 (MQ=255)
                                                                 aGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCttt                          <  2:1691246/91‑1 (MQ=255)
                                                                  ggCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTa                         >  2:1989472/1‑91 (MQ=255)
                                                                  ggCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTa                         >  2:3234901/1‑91 (MQ=255)
                                                                    cTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGc                       >  2:3482802/1‑91 (MQ=255)
                                                                    cTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGc                       <  1:3399448/91‑1 (MQ=255)
                                                                     tAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCt                      <  2:1924251/91‑1 (MQ=255)
                                                                       gATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCttt                    >  1:1079899/1‑91 (MQ=255)
                                                                       gATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCttt                    <  2:3165403/91‑1 (MQ=255)
                                                                       gATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCttt                    >  2:2806859/1‑91 (MQ=255)
                                                                        atatatTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCtttt                   >  1:2835492/1‑91 (MQ=255)
                                                                        atatatTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCtttt                   <  1:131779/91‑1 (MQ=255)
                                                                        atatatTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCtttt                   <  1:23339/91‑1 (MQ=255)
                                                                         tatatTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCttttt                  <  2:2838924/91‑1 (MQ=255)
                                                                         tatatTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCttttt                  <  2:3284563/91‑1 (MQ=255)
                                                                         tatatTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCtattt                  >  1:2344605/1‑91 (MQ=255)
                                                                           tatTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCg                <  1:3535189/91‑1 (MQ=255)
                                                                            atTTGTATTTATAGTATAAATACAAAGAATGTAGAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGc               <  1:120606/91‑1 (MQ=255)
                                                                            atTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGc               >  1:2150046/1‑91 (MQ=255)
                                                                            atTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGc               >  1:3850229/1‑91 (MQ=255)
                                                                             tttGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGcc              <  2:2298398/91‑1 (MQ=255)
                                                                             tttGTATTTATACTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGcc              >  2:2353498/1‑91 (MQ=255)
                                                                              ttGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCa             <  1:3633526/91‑1 (MQ=255)
                                                                               tGTATTTATACTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCAt            >  2:1411219/1‑91 (MQ=255)
                                                                                gTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCAtt           <  2:223248/91‑1 (MQ=255)
                                                                                gTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCAtt           >  2:1653999/1‑91 (MQ=255)
                                                                                   tttATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCAttttt        >  1:686702/1‑91 (MQ=255)
                                                                                   tttATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCAttttt        <  2:920282/91‑1 (MQ=255)
                                                                                   tttATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCAttttt        >  1:805388/1‑91 (MQ=255)
                                                                                   tttATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCAttttt        >  2:2010453/1‑91 (MQ=255)
                                                                                   tttATACTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCAttttt        >  1:3342536/1‑91 (MQ=255)
                                                                                    ttATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCAtttttt       >  2:2977384/1‑91 (MQ=255)
                                                                                      ataGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCATTTTTTcc     >  2:563331/1‑91 (MQ=255)
                                                                                      ataGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCATTTTTTcc     >  1:1653215/1‑91 (MQ=255)
                                                                                      ataGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCATTTTTTcc     >  1:185036/1‑91 (MQ=255)
                                                                                        aGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCATTTTTTCCCt   <  1:2352105/91‑1 (MQ=255)
                                                                                         gTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCATTTTTTCCCtt  <  2:1542832/91‑1 (MQ=255)
                                                                                         gTAAAAACACAAAGAAAGTAAAGACAACTTTGATCTTAATAAGTCATAACGACAACCATGACTCTTTAGCTTTTTCGCCATTTTTGCCCtt  <  2:2617070/91‑1 (MQ=255)
                                                                                         |                                                                                          
TTCGGGATGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCATTTTTTCCCTT  >  USA300TCH1516_ALE/2833298‑2833477

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: