Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 77,975 Δ1 bp 7.0% intergenic (+84/‑343) 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,9750G. 7.0% 830.2 / 52.0 244intergenic (+84/‑343)USA300HOU_RS00350/USA300HOU_RS00355Monoacylglycerol lipase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (118/107);  new base . (5/12);  total (125/119)
Fisher's exact test for biased strand distribution p-value = 8.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

Only 100 of 247 total aligned reads displayed.

TAACTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCA  >  USA300TCH1516_ALE/77885‑78065
                                                                                           |                                                                                          
taACTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACg                                                                                            <  1:3864884/91‑1 (MQ=255)
 aaCTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGa                                                                                           <  2:2938833/91‑1 (MQ=255)
  aCTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGAt                                                                                          <  2:251091/91‑1 (MQ=255)
    tAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAg                                                                                        >  1:1771196/1‑91 (MQ=255)
    tAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAg                                                                                        <  2:217313/91‑1 (MQ=255)
    tAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAg                                                                                        >  2:3695386/1‑91 (MQ=255)
    tAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAg                                                                                        <  1:2416081/91‑1 (MQ=255)
      aaaCAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAt                                                                                      >  2:668411/1‑91 (MQ=255)
       aaCAGTATAGTTCCGTGTATTTGACTATAAGAAATTATGAGGATATTAAAGATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAtt                                                                                     <  2:3603223/91‑1 (MQ=255)
        aCAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAttt                                                                                    <  2:1489254/91‑1 (MQ=255)
        aCAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAttt                                                                                    >  1:1272250/1‑91 (MQ=255)
         cAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAtttt                                                                                   >  1:2107966/1‑91 (MQ=255)
          aGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAttttt                                                                                  >  1:759578/1‑91 (MQ=255)
           gTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAtttttt                                                                                 >  1:2300456/1‑91 (MQ=255)
             ataGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCa                                                                               >  2:2281843/1‑91 (MQ=255)
               aGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATa                                                                             <  2:3406879/91‑1 (MQ=255)
                gTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAg                                                                            <  1:3378065/91‑1 (MQ=21)
                 ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGc                                                                           <  2:3788390/91‑1 (MQ=255)
                 ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGc                                                                           <  1:1752891/91‑1 (MQ=255)
                 ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGc                                                                           <  1:2966255/91‑1 (MQ=255)
                 ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATAATAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGc                                                                           <  2:3484656/91‑1 (MQ=255)
                  tCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCt                                                                          >  1:3312635/1‑91 (MQ=255)
                  tCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCt                                                                          >  1:3439076/1‑91 (MQ=255)
                  tCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCt                                                                          <  2:3186018/91‑1 (MQ=21)
                   ccGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTa                                                                         >  2:1390731/1‑91 (MQ=255)
                    cGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAt                                                                        <  1:1719054/91‑1 (MQ=255)
                    cGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTaa                                                                        <  2:3558103/91‑2 (MQ=18)
                     gtgtATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTaat                                                                       <  2:2819344/91‑3 (MQ=18)
                      tgtATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAttt                                                                      <  1:3657874/91‑1 (MQ=255)
                         aTTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAtttttt                                                                   >  1:1183007/1‑91 (MQ=255)
                          tttGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTtta                                                                  >  1:2280179/1‑91 (MQ=255)
                              aTTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagt                                                              <  1:2642151/91‑1 (MQ=255)
                              aTTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagt                                                              <  2:1149553/91‑1 (MQ=255)
                              aTTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagt                                                              <  1:3757625/91‑1 (MQ=255)
                                tataAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagtta                                                            <  2:100541/91‑1 (MQ=255)
                                 ataAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGCAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttat                                                           >  2:3475970/1‑91 (MQ=255)
                                 ataAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttat                                                           <  2:1944184/91‑1 (MQ=255)
                                  taAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttata                                                          <  1:421007/91‑1 (MQ=255)
                                      aaaTTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATagaga                                                      >  2:3874447/1‑91 (MQ=255)
                                      aaaTTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATagaga                                                      >  2:1955472/1‑91 (MQ=255)
                                       aaTTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATagagag                                                     <  2:1964577/91‑1 (MQ=255)
                                         ttATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGTCGGTAGATTTTTTCATAGCTATTTTTTATAGTTATAGAgagga                                                   <  2:2518665/91‑1 (MQ=255)
                                           aTGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGt                                                 >  2:3320720/1‑91 (MQ=255)
                                             gAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAg                                               >  1:3395310/1‑91 (MQ=255)
                                             gAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAg                                               >  2:3526954/1‑91 (MQ=255)
                                             gAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAg                                               >  1:2292555/1‑91 (MQ=255)
                                              aGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTTATAGCTATTTTTTATAGTTATAGAGAGGAGTAGa                                              >  2:2678982/1‑91 (MQ=255)
                                               ggATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGAc                                             >  1:837430/1‑91 (MQ=255)
                                               ggATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGAc                                             >  2:120217/1‑91 (MQ=255)
                                               ggATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGAc                                             >  2:833680/1‑91 (MQ=255)
                                                   atTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTc                                         <  2:751199/91‑1 (MQ=255)
                                                    ttAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTcc                                        >  1:624948/1‑91 (MQ=255)
                                                    ttAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTcc                                        >  1:2548717/1‑91 (MQ=255)
                                                    ttAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTcc                                        >  2:1291217/1‑91 (MQ=255)
                                                      aaaCATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAg                                      >  2:3571538/1‑91 (MQ=255)
                                                         cATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGAct                                   >  2:1856627/1‑91 (MQ=255)
                                                         cATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGAct                                   >  1:2607304/1‑91 (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)
                                                               aaGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGa                             >  2:1367004/1‑91 (MQ=255)
                                                                aGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAt                            <  2:2138056/91‑1 (MQ=255)
                                                                aGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAt                            <  1:3484368/91‑1 (MQ=255)
                                                                aGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAt                            >  1:1112896/1‑91 (MQ=255)
                                                                aGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAt                            <  1:2126401/91‑1 (MQ=255)
                                                                aGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAt                            >  1:1486946/1‑91 (MQ=255)
                                                                 gATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAtt                           >  2:1668862/1‑91 (MQ=255)
                                                                 gATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAtt                           <  1:2421542/91‑1 (MQ=255)
                                                                 gATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAtt                           >  2:437998/1‑91 (MQ=255)
                                                                 aaatAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAtt                           >  2:563238/4‑91 (MQ=17)
                                                                  aatAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAttt                          >  2:2673880/3‑91 (MQ=17)
                                                                   atAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAtttt                         >  1:536936/2‑91 (MQ=17)
                                                                     aGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCCGACTCTTGGATTTTaa                       >  2:3708859/1‑91 (MQ=255)
                                                                       cTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAAt                     >  2:421622/1‑91 (MQ=255)
                                                                       cTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAAt                     >  1:2470617/1‑91 (MQ=255)
                                                                        tATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATc                    <  1:1072397/91‑1 (MQ=255)
                                                                         aTGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATcc                   >  2:907138/1‑91 (MQ=255)
                                                                          tGAAGAAATCTATGACGNTAGATNTTTTCATAGCTATTTTTTANAGTTATAGAGAGGAGNAGACTGTCCAGACTCTTGGATTTTAAATNCg                  <  2:1020852/91‑1 (MQ=17)
                                                                          tGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCg                  <  1:838564/91‑1 (MQ=255)
                                                                          tGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCg                  <  1:838466/91‑1 (MQ=255)
                                                                           gaagaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGt                 >  1:100186/1‑91 (MQ=255)
                                                                            aagaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTa                >  1:542583/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)
                                                                                    tATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaaaaaaaaa        <  1:2859455/91‑2 (MQ=255)
                                                                                    tATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaaaaaaaaa        <  1:932701/91‑2 (MQ=255)
                                                                                     aTGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACa       >  1:3201488/1‑91 (MQ=255)
                                                                                     aTGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACa       >  2:3050457/1‑91 (MQ=255)
                                                                                      tGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACaa      >  1:1866217/1‑91 (MQ=255)
                                                                                        gACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAg     <  2:1644069/91‑1 (MQ=255)
                                                                                        gACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAg     <  2:3273593/91‑1 (MQ=255)
                                                                                        gACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAg     >  1:3337180/1‑91 (MQ=255)
                                                                                         aCGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAACAAGTc   <  2:2442998/91‑1 (MQ=255)
                                                                                           gATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCa  <  1:267964/91‑1 (MQ=255)
                                                                                           |                                                                                          
TAACTAAAACAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCA  >  USA300TCH1516_ALE/77885‑78065

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: