Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,313,118 A→T 5.6% intergenic (+22/+275) rny_1 → / ← USA300HOU_RS06485 Ribonuclease Y/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,313,1180AT 5.6% 431.8 / 13.4 124intergenic (+22/+275)rny_1/USA300HOU_RS06485Ribonuclease Y/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (60/57);  new base T (5/2);  total (65/59)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 132 total aligned reads displayed.

TTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGAAAGTA  >  USA300TCH1516_ALE/1313028‑1313208
                                                                                          |                                                                                          
ttACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTa                                                                                            <  1:3192168/91‑1 (MQ=255)
ttACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTa                                                                                            <  1:212539/91‑1 (MQ=255)
 tACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAg                                                                                           <  1:182824/91‑1 (MQ=255)
 tACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAg                                                                                           <  1:185188/91‑1 (MQ=255)
     atatCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAg                                                                                       >  1:3839831/1‑91 (MQ=255)
     atatCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAg                                                                                       >  1:3841279/1‑91 (MQ=255)
     atatCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAg                                                                                       >  1:1776098/1‑91 (MQ=255)
     atatCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAg                                                                                       <  1:2048087/91‑1 (MQ=255)
      tatCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAgg                                                                                      >  1:1444372/1‑91 (MQ=255)
         ccTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAg                                                                                   <  1:3153267/91‑1 (MQ=255)
         ccTGGTCATATCAAGGTCACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAg                                                                                   <  1:3551049/91‑1 (MQ=255)
          cTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGc                                                                                  <  1:179279/91‑1 (MQ=255)
           tGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCt                                                                                 >  1:2910140/1‑91 (MQ=255)
            ggTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACTAATTTGTGAGGGAgaca                                                                                <  1:3951028/91‑4 (MQ=39)
             gTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCttt                                                                               >  1:2017538/1‑91 (MQ=255)
               cATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCttttt                                                                             <  1:3777419/91‑1 (MQ=255)
                atatCAAGGTGACAGTTGTCCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCtttttt                                                                            >  1:519586/1‑91 (MQ=255)
                 tatgaaGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTa                                                                           <  1:3734007/87‑1 (MQ=255)
                 tatCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTa                                                                           <  1:266885/91‑1 (MQ=255)
                 tatCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTa                                                                           <  1:1446104/91‑1 (MQ=255)
                   tCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAg                                                                         <  1:3870980/91‑1 (MQ=255)
                    cAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGt                                                                        >  1:3813589/1‑91 (MQ=255)
                    cAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAAAAATTTTTGTCTCCCTCACCAATTAGTGAGGGAGCGTTTTTAAGt                                                                        >  1:894728/1‑91 (MQ=255)
                     aaGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTGTGTCTCCCTCACAGATTAGTGAGGGAGCTTTTTCAAGtt                                                                       >  1:3865672/1‑91 (MQ=255)
                      aGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTg                                                                      <  1:1899412/91‑1 (MQ=255)
                       ggTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGt                                                                     <  1:3838483/91‑1 (MQ=255)
                       ggTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGt                                                                     >  1:1810944/1‑91 (MQ=255)
                       ggTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGt                                                                     >  1:3893292/1‑91 (MQ=255)
                        gTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTa                                                                    <  1:3998692/91‑1 (MQ=255)
                         tGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAg                                                                   >  1:2589789/1‑91 (MQ=255)
                          gACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGt                                                                  <  1:1159103/91‑1 (MQ=255)
                          gACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGt                                                                  <  1:1159440/91‑1 (MQ=255)
                           aCAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTc                                                                 >  1:1697130/1‑91 (MQ=255)
                           aCAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTc                                                                 >  1:2815755/1‑91 (MQ=255)
                           aCAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTc                                                                 >  1:1696278/1‑91 (MQ=255)
                           aCAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTc                                                                 <  1:1901417/91‑1 (MQ=255)
                            cAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCt                                                                <  1:1494068/91‑1 (MQ=255)
                            cAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCt                                                                >  1:4132924/1‑91 (MQ=255)
                            cAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCt                                                                >  1:2694803/1‑91 (MQ=255)
                              gttgttCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTa                                                              >  1:2262429/1‑91 (MQ=255)
                               ttgttCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTaa                                                             <  1:3349987/91‑1 (MQ=255)
                               ttgttCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTaa                                                             >  1:1700107/1‑91 (MQ=255)
                                 gttCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATc                                                           <  1:2500779/91‑1 (MQ=255)
                                     gagagaCTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGt                                                       <  1:3815003/91‑1 (MQ=255)
                                      agaggCAAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGtt                                                      <  1:1118116/91‑1 (MQ=255)
                                      agagaCTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGtt                                                      >  1:1092491/1‑91 (MQ=255)
                                       gagaCTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTa                                                     >  1:1205890/1‑91 (MQ=255)
                                           cTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACa                                                 >  1:1029218/1‑91 (MQ=255)
                                            tAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAg                                                <  1:2667820/91‑1 (MQ=255)
                                            tAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAg                                                <  1:2015403/91‑1 (MQ=255)
                                             agagCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGc                                               <  1:1611415/91‑1 (MQ=255)
                                              gagCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCa                                              >  1:3020246/1‑91 (MQ=255)
                                               agCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCAc                                             >  1:2778924/1‑91 (MQ=255)
                                               agCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCAc                                             >  1:727524/1‑91 (MQ=255)
                                               agCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCAc                                             <  1:30869/91‑1 (MQ=255)
                                                gCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACt                                            >  1:1896473/1‑91 (MQ=255)
                                                 cAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACtt                                           >  1:2401230/1‑91 (MQ=255)
                                                   gTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTa                                         <  1:1041038/91‑1 (MQ=255)
                                                   gTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTa                                         <  1:1037916/91‑1 (MQ=255)
                                                    tAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTAt                                        <  1:1773338/91‑1 (MQ=255)
                                                    tAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTAt                                        <  1:812199/91‑1 (MQ=255)
                                                    tAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTAt                                        <  1:812345/91‑1 (MQ=255)
                                                    tAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTAt                                        <  1:2831519/91‑1 (MQ=255)
                                                       aaTATGCGAAATCATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTAACgg                                     >  1:3209623/1‑91 (MQ=255)
                                                        atatGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGt                                    <  1:3703185/91‑1 (MQ=255)
                                                          atGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGtaa                                  >  1:4047042/1‑91 (MQ=255)
                                                          atGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAAGGAGCTTTTTTAAGTTGTAGTCTTAATTTAGTTAGACAGCACTTTATAGGtaa                                  >  1:3865950/1‑91 (MQ=255)
                                                           tGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGtaat                                 <  1:3472032/91‑1 (MQ=255)
                                                            gCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGtaata                                <  1:1757036/91‑1 (MQ=255)
                                                             cGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGtaataa                               >  1:760198/1‑91 (MQ=255)
                                                              gAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTATTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAAc                              <  1:3497040/91‑1 (MQ=255)
                                                              gAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAAc                              <  1:1604337/91‑1 (MQ=255)
                                                              gAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAAc                              <  1:1317382/91‑1 (MQ=255)
                                                                   aatTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTatat                         <  1:3380785/91‑1 (MQ=255)
                                                                   aatTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTatat                         <  1:3379494/91‑1 (MQ=255)
                                                                      ttttGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTaa                      >  1:460041/1‑91 (MQ=255)
                                                                      ttttGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTaa                      >  1:24361/1‑91 (MQ=255)
                                                                      ttttGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACCATATTaa                      >  1:2066740/1‑91 (MQ=255)
                                                                       tttGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTaaa                     >  1:280270/1‑91 (MQ=255)
                                                                           gctcCCTCACTAATTTGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATGTTAAACagt                 >  1:2641694/2‑91 (MQ=39)
                                                                            ctcCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagta                <  1:3814026/91‑1 (MQ=255)
                                                                             tcCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagtag               >  1:1589877/1‑91 (MQ=255)
                                                                             tcCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagtag               >  1:4018261/1‑91 (MQ=255)
                                                                             tcCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagtag               >  1:109314/1‑91 (MQ=255)
                                                                              cccTCACTAATTTGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagtagt              >  1:2897435/1‑91 (MQ=255)
                                                                              cccTCACAAATTAGTGGGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagtagt              >  1:1209444/1‑91 (MQ=255)
                                                                              cccTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagtagt              >  1:3239860/1‑91 (MQ=255)
                                                                               ccTCACTAATTTGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGtt             >  1:1001339/1‑91 (MQ=255)
                                                                                cTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTCAATCTAATTAGACAGCACTTTATCGGTAATAACTATAGTAAACAGTGGTTa            >  1:3203182/1‑91 (MQ=255)
                                                                                  cactAATTTGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTAtt          >  1:760715/1‑91 (MQ=255)
                                                                                    caAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTg        >  1:2889343/1‑91 (MQ=255)
                                                                                      aaTTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGaa      >  1:1191871/1‑91 (MQ=255)
                                                                                      aaTTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGaa      >  1:3242481/1‑91 (MQ=255)
                                                                                       aTTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGaaa     >  1:2318630/1‑91 (MQ=255)
                                                                                         tAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGAAAGt   >  1:3491404/1‑91 (MQ=255)
                                                                                         tAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGAAAGt   >  1:2589826/1‑91 (MQ=255)
                                                                                          aGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGAAAGTa  >  1:4023252/1‑91 (MQ=255)
                                                                                          aGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGAAAGTa  >  1:3738942/1‑91 (MQ=255)
                                                                                          aGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGAAAGTa  >  1:2469740/1‑91 (MQ=255)
                                                                                          aGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGAAAGTa  >  1:2560361/1‑91 (MQ=255)
                                                                                          |                                                                                          
TTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTGAAAGTA  >  USA300TCH1516_ALE/1313028‑1313208

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: ‑