Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 77,971:1 +C 7.1% intergenic (+80/‑347) 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,9711.C 7.1% 1030.8 / 71.4 238intergenic (+80/‑347)USA300HOU_RS00350/USA300HOU_RS00355Monoacylglycerol lipase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (116/105);  new base C (5/12);  total (121/117)
Fisher's exact test for biased strand distribution p-value = 8.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.13e-01

Only 100 of 244 total aligned reads displayed.

CTTATAACTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAA  >  USA300TCH1516_ALE/77881‑78061
                                                                                           |                                                                                          
cTTATAACTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTAt                                                                                             <  1:1660616/91‑1 (MQ=255)
 ttATAACTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATg                                                                                           >  2:580569/1‑91 (MQ=255)
    taACTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTCTGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACg                                                                                        <  1:3869715/91‑1 (MQ=255)
     aaCTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGa                                                                                       <  2:1483364/91‑1 (MQ=255)
     aaCTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGa                                                                                       <  2:2938833/91‑1 (MQ=255)
      aCTAAAACAGTATAGTTCCGTGTATTTTATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGAt                                                                                      <  2:421185/91‑1 (MQ=255)
      aCTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGAt                                                                                      <  2:251091/91‑1 (MQ=255)
       cTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATa                                                                                     >  1:731158/1‑91 (MQ=255)
        tAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAg                                                                                    <  2:1584933/91‑1 (MQ=255)
        tAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAg                                                                                    <  1:856913/91‑1 (MQ=255)
        tAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAg                                                                                    >  1:541582/1‑91 (MQ=255)
         aaaaCAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTTAGATTAGCTATGAAGAAATCTATGACGATAGa                                                                                   >  1:3297932/1‑91 (MQ=255)
          aaaCAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAt                                                                                  >  2:668411/1‑91 (MQ=255)
          aaaCAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAt                                                                                  <  1:76142/91‑1 (MQ=255)
           aaCAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAtt                                                                                 >  1:596746/1‑91 (MQ=255)
            aCAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAttt                                                                                >  1:1272250/1‑91 (MQ=255)
             cAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAtttt                                                                               >  1:2108949/1‑91 (MQ=255)
               gTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAtttttt                                                                             >  1:2300456/1‑91 (MQ=255)
                tataGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTc                                                                            <  2:3180001/91‑1 (MQ=255)
                 ataGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAtttttttc                                                                           >  2:2156943/1‑89 (MQ=255)
                 ataGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCa                                                                           >  2:3787068/1‑91 (MQ=255)
                 ataGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCa                                                                           >  2:32792/1‑91 (MQ=255)
                   aGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATa                                                                         <  1:3036191/91‑1 (MQ=21)
                    gTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAg                                                                        <  1:1991842/91‑1 (MQ=21)
                    gTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAg                                                                        <  1:3378065/91‑1 (MQ=21)
                     ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGc                                                                       <  1:1752749/91‑1 (MQ=255)
                      tCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCt                                                                      >  1:3312864/1‑91 (MQ=255)
                      tCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCt                                                                      >  2:1372561/1‑91 (MQ=255)
                      tCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCt                                                                      >  1:3439076/1‑91 (MQ=255)
                      tCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCt                                                                      <  1:2292167/91‑1 (MQ=21)
                        cGTGTATTTGATTATAAGACATTATGAGGATATTAAAGATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAt                                                                    >  2:3520393/1‑91 (MQ=255)
                        cGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAt                                                                    <  2:122189/91‑1 (MQ=255)
                        cGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAt                                                                    >  2:3476304/1‑91 (MQ=255)
                         gtgtATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAg                                                                        <  2:439466/86‑1 (MQ=255)
                         gtgtATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAtt                                                                   >  1:238495/1‑91 (MQ=255)
                         gtgtATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTaat                                                                   <  2:2819344/91‑3 (MQ=18)
                          tgtATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAttt                                                                  <  1:3657874/91‑1 (MQ=255)
                             aTTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAtttttt                                                               >  1:1183007/1‑91 (MQ=255)
                              tttGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTtta                                                              >  2:3588362/1‑91 (MQ=255)
                                tGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatt                                                            >  2:426604/1‑90 (MQ=255)
                                  aTTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagt                                                          <  1:3757625/91‑1 (MQ=255)
                                  aTTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagt                                                          <  2:3435126/91‑1 (MQ=255)
                                    tataAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAGTCTATGACGATAGATTTTTTAATAGCTATTTTttatagtta                                                        <  2:101412/91‑1 (MQ=255)
                                     ataAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttat                                                       <  2:1944184/91‑1 (MQ=255)
                                      taAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCAt                                                                          >  2:1190687/1‑71 (MQ=255)
                                      taAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttata                                                      >  1:1753012/1‑91 (MQ=255)
                                      taAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttata                                                      <  2:1832985/91‑1 (MQ=255)
                                         gAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATagag                                                   <  2:1024883/91‑1 (MQ=255)
                                          aaaTTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATagaga                                                  >  2:1955472/1‑91 (MQ=255)
                                           aaTTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATagagag                                                 <  2:1964577/91‑1 (MQ=255)
                                             ttATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGTCGGTAGATTTTTTCATAGCTATTTTTTATAGTTATAGAgagga                                               <  2:2518665/91‑1 (MQ=255)
                                             ttATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATTGCTATTTTTTATAGTTATAGAgagga                                               <  2:2217044/91‑1 (MQ=255)
                                             ttATGAGGATATAAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAgaggt                                               >  1:2980298/1‑90 (MQ=255)
                                              tATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAgaggag                                              <  1:2229015/91‑1 (MQ=255)
                                               aTGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGt                                             >  2:3320720/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)
                                                       atTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTc                                     <  2:751199/91‑1 (MQ=255)
                                                        ttAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTcc                                    <  1:1830812/91‑1 (MQ=255)
                                                             cATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGAct                               >  2:1856627/1‑91 (MQ=255)
                                                             cATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGAct                               >  1:2607304/1‑91 (MQ=255)
                                                               tACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATGGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGActct                             <  2:2078207/91‑1 (MQ=255)
                                                               tACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGActct                             >  1:2623302/1‑91 (MQ=255)
                                                                 cTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTg                           >  2:3045412/1‑91 (MQ=255)
                                                                  tAAGATTAGCTTTGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTgg                          <  1:3519033/91‑1 (MQ=255)
                                                                  tAAGAGTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTgg                          <  2:3569031/91‑1 (MQ=255)
                                                                   aaGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGa                         <  2:931320/91‑1 (MQ=255)
                                                                    aGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAt                        <  2:3543353/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                       <  1:2421542/91‑1 (MQ=255)
                                                                      aatAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAttt                      >  2:2673880/3‑91 (MQ=17)
                                                                      aTTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAttt                      <  1:493453/91‑1 (MQ=255)
                                                                       ttAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAtttt                     <  2:2300456/91‑1 (MQ=255)
                                                                           cTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAAt                 >  1:2470617/1‑91 (MQ=255)
                                                                           cTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTGAAAt                 >  1:2512827/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)
                                                                             aTGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATcc               >  2:3697193/1‑91 (MQ=255)
                                                                              tGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCg              <  1:838269/91‑1 (MQ=255)
                                                                              tGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCg              <  1:838507/91‑1 (MQ=255)
                                                                               gaagaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGt             >  1:100186/1‑91 (MQ=255)
                                                                               gaagaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGt             >  1:95056/1‑91 (MQ=255)
                                                                                aagaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTa            >  1:2702281/1‑91 (MQ=255)
                                                                                aagaaATCTATGACGAGAGATTTTTGCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTa            <  2:1101084/91‑1 (MQ=255)
                                                                                 agaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaa           <  2:3629105/91‑1 (MQ=255)
                                                                                 agaaATCTATGACGATAGATTTTTCCATAGCTATTTTTTATAGCTATAGAAAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaa           <  2:3611959/91‑1 (MQ=255)
                                                                                 agaaATCTAAGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCATaa           <  2:1346070/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)
                                                                                      tCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaaaaaaa      >  1:3485394/1‑91 (MQ=255)
                                                                                        tATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaaaaaaaaa    <  1:2859455/91‑2 (MQ=255)
                                                                                         aTGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGAct                               <  2:3857081/63‑1 (MQ=255)
                                                                                         aTGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGAct                               >  1:3857081/1‑63 (MQ=255)
                                                                                         aTGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACa   >  1:3201488/1‑91 (MQ=255)
                                                                                         aTGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACa   >  2:3050457/1‑91 (MQ=255)
                                                                                           |                                                                                          
CTTATAACTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAA  >  USA300TCH1516_ALE/77881‑78061

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