Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 77,973 A→T 7.2% intergenic (+82/‑345) USA300HOU_RS00350 → / → USA300HOU_RS00355 Monoacylglycerol lipase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE77,9730AT 7.2% 799.9 / 36.8 247intergenic (+82/‑345)USA300HOU_RS00350/USA300HOU_RS00355Monoacylglycerol lipase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (120/108);  new base T (5/13);  total (126/121)
Fisher's exact test for biased strand distribution p-value = 5.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 248 total aligned reads displayed.

ATAACTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTC  >  USA300TCH1516_ALE/77884‑78064
                                                                                          |                                                                                           
ataACTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGAc                                                                                            >  1:2996764/1‑91 (MQ=255)
 taACTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACg                                                                                           <  1:3864884/91‑1 (MQ=255)
  aaCTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGa                                                                                          <  2:2938833/91‑1 (MQ=255)
   aCTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGAt                                                                                         <  2:251091/91‑1 (MQ=255)
     tAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAg                                                                                       <  2:217313/91‑1 (MQ=255)
     tAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAg                                                                                       >  2:3695386/1‑91 (MQ=255)
     tAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAg                                                                                       <  1:2416081/91‑1 (MQ=255)
     tAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAg                                                                                       >  1:1771196/1‑91 (MQ=255)
     tAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAg                                                                                       <  2:1584933/91‑1 (MQ=255)
       aaaCAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAt                                                                                     >  2:668411/1‑91 (MQ=255)
       aaaCAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAt                                                                                     <  1:76142/91‑1 (MQ=255)
        aaCAGTATAGTTCCGTGTATTTGACTATAAGAAATTATGAGGATATTAAAGATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAtt                                                                                    <  2:3603223/91‑1 (MQ=255)
         aCAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAttt                                                                                   >  1:1272250/1‑91 (MQ=255)
          cAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAtttt                                                                                  >  1:2107966/1‑91 (MQ=255)
            gTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAtttttt                                                                                >  1:2300456/1‑91 (MQ=255)
             tataGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTc                                                                               <  2:3180001/91‑1 (MQ=255)
                aGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATa                                                                            <  2:3406879/91‑1 (MQ=255)
                 gTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAg                                                                           <  1:3378065/91‑1 (MQ=21)
                  ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATGGATTTTTTCATAGc                                                                          <  2:3300652/91‑1 (MQ=255)
                  ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGc                                                                          <  2:3788390/91‑1 (MQ=255)
                  ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGc                                                                          <  1:588506/91‑1 (MQ=255)
                  ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGc                                                                          <  1:1752891/91‑1 (MQ=255)
                  ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGc                                                                          <  1:2965726/91‑1 (MQ=255)
                  ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATAATAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGc                                                                          <  2:3484656/91‑1 (MQ=255)
                   tCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCt                                                                         >  2:1372561/1‑91 (MQ=255)
                   tCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCt                                                                         >  1:3439076/1‑91 (MQ=255)
                   tCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCt                                                                         >  1:3312635/1‑91 (MQ=255)
                    ccGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTa                                                                        >  2:1390731/1‑91 (MQ=255)
                     cGTGTATTTGATTATAAGACATTATGAGGATATTAAAGATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAt                                                                       >  2:3520393/1‑91 (MQ=255)
                     cGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAt                                                                       <  2:122189/91‑1 (MQ=255)
                     cGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAt                                                                       <  1:1719054/91‑1 (MQ=255)
                     cGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTaa                                                                       <  2:3558103/91‑2 (MQ=18)
                      gtgtATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAg                                                                           >  1:439466/1‑86 (MQ=255)
                      gtgtATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTaat                                                                      <  2:2819344/91‑3 (MQ=18)
                         tATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAttttt                                                                   <  1:3686622/91‑1 (MQ=255)
                          aTTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAtttttt                                                                  >  1:1183007/1‑91 (MQ=255)
                           tttGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTtta                                                                 >  1:2280179/1‑91 (MQ=255)
                               aTTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagt                                                             <  1:3757961/91‑1 (MQ=255)
                               aTTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagt                                                             <  1:2642151/91‑1 (MQ=255)
                                 tataAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagtta                                                           <  2:100541/91‑1 (MQ=255)
                                   taAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCAt                                                                             >  2:1190687/1‑71 (MQ=255)
                                   taAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttata                                                         <  2:1832985/91‑1 (MQ=255)
                                   taAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttata                                                         <  2:1953272/91‑1 (MQ=255)
                                    aaGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttatag                                                        >  1:2612980/1‑91 (MQ=255)
                                      gAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATagag                                                      <  2:1024883/91‑1 (MQ=255)
                                       aaaTTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATagaga                                                     >  2:3874447/1‑91 (MQ=255)
                                       aaaTTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATagaga                                                     >  2:1955472/1‑91 (MQ=255)
                                          ttATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGTCGGTAGATTTTTTCATAGCTATTTTTTATAGTTATAGAgagga                                                  <  2:2518665/91‑1 (MQ=255)
                                          ttATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATTGCTATTTTTTATAGTTATAGAgagga                                                  <  2:2217044/91‑1 (MQ=255)
                                              gAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAg                                              >  1:2292555/1‑91 (MQ=255)
                                              gAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAg                                              >  1:3395310/1‑91 (MQ=255)
                                               aGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTTATAGCTATTTTTTATAGTTATAGAGAGGAGTAGa                                             >  2:2678982/1‑91 (MQ=255)
                                                ggATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGAc                                            >  2:833680/1‑91 (MQ=255)
                                                ggATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGAc                                            >  1:837430/1‑91 (MQ=255)
                                                    atTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTc                                        <  2:751199/91‑1 (MQ=255)
                                                     ttAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTcc                                       >  2:1291217/1‑91 (MQ=255)
                                                     ttAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTcc                                       >  1:624948/1‑91 (MQ=255)
                                                     ttAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTcc                                       >  1:2548717/1‑91 (MQ=255)
                                                       aaaCATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAg                                     >  2:3571538/1‑91 (MQ=255)
                                                          cATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGAct                                  >  2:1856627/1‑91 (MQ=255)
                                                           aTACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGActc                                 <  1:3297641/91‑1 (MQ=255)
                                                           aTACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGActc                                 <  2:1972868/91‑1 (MQ=255)
                                                            tACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATGGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGActct                                <  2:2078207/91‑1 (MQ=255)
                                                              cTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTg                              >  2:3045412/1‑91 (MQ=255)
                                                               tAAGATTAGCTTTGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTgg                             <  1:3519033/91‑1 (MQ=255)
                                                                aaGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGa                            <  2:931320/91‑1 (MQ=255)
                                                                aaGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGa                            >  1:3085584/1‑91 (MQ=255)
                                                                 aGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAt                           <  1:2126401/91‑1 (MQ=255)
                                                                 aGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAt                           <  1:3484368/91‑1 (MQ=255)
                                                                 aGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAt                           >  1:1486946/1‑91 (MQ=255)
                                                                 aGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAt                           >  1:1112896/1‑91 (MQ=255)
                                                                  gATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAtt                          >  2:437998/1‑91 (MQ=255)
                                                                  gATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAtt                          <  1:2421542/91‑1 (MQ=255)
                                                                  aaatAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAtt                          >  2:563238/4‑91 (MQ=17)
                                                                   aatAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAttt                         >  2:2673880/3‑91 (MQ=17)
                                                                   aTTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAttt                         <  1:2812970/91‑1 (MQ=255)
                                                                    atAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAtttt                        >  1:536936/2‑91 (MQ=17)
                                                                      aGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCCGACTCTTGGATTTTaa                      >  2:3708859/1‑91 (MQ=255)
                                                                        cTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAAt                    >  2:421622/1‑91 (MQ=255)
                                                                         tATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATc                   <  1:1072397/91‑1 (MQ=255)
                                                                          aTGAAGAAATCTATGACGATAGATTTTTTCATTGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATcc                  >  2:2386706/1‑91 (MQ=255)
                                                                          aTGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATcc                  >  2:907138/1‑91 (MQ=255)
                                                                           tGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCg                 <  1:838507/91‑1 (MQ=255)
                                                                           tGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCg                 <  1:838564/91‑1 (MQ=255)
                                                                           tGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCg                 <  1:2493558/91‑1 (MQ=255)
                                                                            gaagaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGt                >  1:95056/1‑91 (MQ=255)
                                                                            gaagaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGt                >  1:100186/1‑91 (MQ=255)
                                                                             aagaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATACGTa               <  2:1099938/91‑1 (MQ=255)
                                                                              agaaATCTATGACGATAGATTTTCTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaa              <  2:3613280/91‑1 (MQ=255)
                                                                               gaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaaa             <  1:3728070/91‑1 (MQ=255)
                                                                                 aaTCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaaaaa           >  2:2690268/1‑91 (MQ=255)
                                                                                  aTCTATGACGATAGATTTTTTCATAGCTATTTTTTATTGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaaaaaa          >  2:2739711/1‑91 (MQ=255)
                                                                                  aTCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaaaaaa          >  2:2333834/1‑91 (MQ=255)
                                                                                   tCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaaaaaaa         >  1:3485394/1‑91 (MQ=255)
                                                                                      aTGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACa      >  2:3050457/1‑91 (MQ=255)
                                                                                       tGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACaa     >  1:1866217/1‑91 (MQ=255)
                                                                                         gACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAg    >  1:3337180/1‑91 (MQ=255)
                                                                                         gACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAg    <  2:3273593/91‑1 (MQ=255)
                                                                                         gACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAg    <  2:1644069/91‑1 (MQ=255)
                                                                                         gACAATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAATAAACAAg    >  2:3441415/1‑91 (MQ=255)
                                                                                          |                                                                                           
ATAACTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTC  >  USA300TCH1516_ALE/77884‑78064

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