Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,931,599 A→T 16.4% L1217I (TTA→ATA)  USA300HOU_RS09605 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,931,5990AT16.4% 1365.9 / 225.6 511L1217I (TTA→ATA) USA300HOU_RS09605hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (244/183);  new base T (58/26);  total (302/209)
Fisher's exact test for biased strand distribution p-value = 5.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 566 total aligned reads displayed.

TGTACCACCATTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAGTTTTCCA  >  USA300TCH1516_ALE/1931509‑1931689
                                                                                          |                                                                                           
tGTACCACCATTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTtaa                                                                                             <  2:1046117/91‑1 (MQ=255)
  tACCACCATTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTtaata                                                                                           <  2:545471/91‑1 (MQ=255)
  tACCACCATTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTtaata                                                                                           <  2:225260/91‑1 (MQ=255)
   accaccATTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTtattaa                                                                                          <  1:2514243/91‑1 (MQ=255)
   accaccATTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTtaataa                                                                                          >  1:2010059/1‑91 (MQ=255)
     caccaTTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTATTAAAc                                                                                        <  1:2218548/91‑1 (MQ=255)
     caccaTTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAAc                                                                                        <  1:1021066/91‑1 (MQ=255)
     caccaTTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAAc                                                                                        <  2:3633924/91‑1 (MQ=255)
      accaTTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAAct                                                                                       <  1:54196/91‑1 (MQ=255)
       ccaTTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAActt                                                                                      >  2:3270267/1‑91 (MQ=255)
        caTTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAGTAGACAATTTTTTAGTATTAAGAAAGTTTATTAAActtt                                                                                     <  2:1367835/91‑1 (MQ=255)
        caTTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAActtt                                                                                     <  1:155306/91‑1 (MQ=255)
        caTTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAActtt                                                                                     <  1:2924668/91‑1 (MQ=255)
         aTTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAActttc                                                                                    <  1:302077/91‑1 (MQ=255)
         aTTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAActttc                                                                                    <  2:616776/91‑1 (MQ=255)
         aTTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAActttc                                                                                    >  2:89177/1‑91 (MQ=255)
           tCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAActttc‑tt                                                                                 <  2:1438988/91‑1 (MQ=255)
           tCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAActttc‑tt                                                                                 <  1:1882413/91‑1 (MQ=255)
            cTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTGATAAActttc‑ttt                                                                                >  2:2684611/1‑91 (MQ=255)
            cTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAActttc‑ttt                                                                                >  1:2635921/1‑91 (MQ=255)
             ttGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑tttt                                                                               <  2:2478668/91‑1 (MQ=255)
              tGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑ttttt                                                                              >  2:2378397/1‑91 (MQ=255)
              tGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTGTAATAAACTTTC‑ttttt                                                                              >  2:1823/1‑91 (MQ=255)
               gTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑tttttt                                                                             <  1:466800/91‑1 (MQ=255)
               gTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑tttttt                                                                             >  1:3739877/1‑91 (MQ=255)
                      tcGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGt                                                                      <  1:1704803/91‑1 (MQ=255)
                         aaGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTaaa                                                                   >  2:3483878/1‑91 (MQ=255)
                         aaGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTaaa                                                                   >  2:1187624/1‑91 (MQ=255)
                          agagAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTCTAATACACTTTC‑TTTCTTAATAAGTAAAc                                                                  >  2:58187/1‑91 (MQ=255)
                           gagaAATTTTTTCTAGATATATTTTCCATCTAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAAca                                                                 <  2:2181486/91‑1 (MQ=255)
                           gagaAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAAca                                                                 <  2:222678/91‑1 (MQ=255)
                              aaaTTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAAcattc                                                              <  2:308094/91‑1 (MQ=255)
                                aTTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAAcattcat                                                            >  2:2108533/1‑91 (MQ=255)
                                aTTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAAcattcat                                                            <  2:1707716/91‑1 (MQ=255)
                                aTTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAAcattcat                                                            >  2:1677037/1‑91 (MQ=255)
                                aTTTTTTATAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAAcattcat                                                            <  1:1963612/91‑1 (MQ=255)
                                 ttttttCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAAcattcatt                                                           <  2:3021659/91‑1 (MQ=255)
                                 ttttttCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAAcattcatt                                                           <  1:69982/91‑1 (MQ=255)
                                 ttttttCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAAcattcatt                                                           <  1:3356555/91‑1 (MQ=255)
                                 ttttttCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAAcattcatt                                                           <  2:416825/91‑1 (MQ=255)
                                   ttttCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTAt                                                         <  2:2751250/91‑1 (MQ=255)
                                   ttttCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTAt                                                         <  2:1280314/91‑1 (MQ=255)
                                    tttGTAGATCTATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATc                                                        <  1:3557033/91‑1 (MQ=255)
                                      tCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCaa                                                      >  1:3488194/1‑91 (MQ=255)
                                       cTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCaaa                                                     >  2:3297300/1‑91 (MQ=255)
                                       cTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCaaa                                                     >  1:1775490/1‑91 (MQ=255)
                                       cTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCaaa                                                     >  1:2769400/1‑91 (MQ=255)
                                       cTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCaaa                                                     >  1:3202831/1‑91 (MQ=255)
                                           atatatTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCg                                                 <  1:2364113/91‑1 (MQ=255)
                                             atatTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAg                                               <  1:926295/91‑1 (MQ=255)
                                              tatTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGt                                              <  1:10719/91‑1 (MQ=255)
                                              tatTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGt                                              <  1:2603988/91‑1 (MQ=255)
                                               atTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTa                                             >  2:2945933/1‑91 (MQ=255)
                                                  ttCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAAtt                                          <  2:1594732/91‑1 (MQ=255)
                                                  ttCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAAtt                                          >  1:94195/1‑91 (MQ=255)
                                                   tccatccaATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAAttt                                         <  1:1388177/91‑1 (MQ=255)
                                                   tatatccaATAGACAATTTTTTTGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAAttt                                         >  2:3365978/4‑91 (MQ=255)
                                                     catccaATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAg                                       >  2:2020551/1‑91 (MQ=255)
                                                      atccaATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGa                                      >  2:3766744/1‑91 (MQ=255)
                                                       tccaATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTCTTTTTTTAATAATTAAACATTCATTATCAACTTCGAGTCATTTAGa                                      >  2:1062507/1‑91 (MQ=255)
                                                         caATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAAtt                                   >  1:77455/1‑91 (MQ=255)
                                                            tAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAAt                                >  2:2637964/1‑91 (MQ=255)
                                                             aGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGCGTAATTTAGAATTAAtt                               >  2:3930740/1‑90 (MQ=255)
                                                              gACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATAc                              <  1:2455085/91‑1 (MQ=255)
                                                              gACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATAc                              <  1:1338043/91‑1 (MQ=255)
                                                                cAATTTTTTTTTATTTAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACat                            >  2:765117/1‑91 (MQ=255)
                                                                cAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACat                            <  1:2082412/91‑1 (MQ=255)
                                                                   ttttttAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATAtt                         <  2:1951807/91‑1 (MQ=255)
                                                                    tttttAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATAttt                        >  2:3842781/1‑90 (MQ=255)
                                                                     ttttAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGCGTAATTTAGAATTAATACATATTGt                       >  1:3668244/1‑91 (MQ=255)
                                                                     ttttAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGt                       <  1:2836243/91‑1 (MQ=255)
                                                                          gTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGCGTAATTTAGAATTAATACATATTGTTTCAc                  >  2:250786/1‑91 (MQ=255)
                                                                          gTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCAc                  >  1:395676/1‑91 (MQ=255)
                                                                          gTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCAc                  >  2:2871531/1‑91 (MQ=255)
                                                                          gTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCAc                  <  1:1184974/91‑1 (MQ=255)
                                                                          gTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCAc                  >  1:3240764/1‑91 (MQ=255)
                                                                           tATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCAcc                 >  1:3072648/1‑91 (MQ=255)
                                                                             ttATGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTACTTTAGAATTAATACATATTGTTTCACCaa               >  2:721639/1‑91 (MQ=255)
                                                                                aGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAAtt            <  2:3579222/91‑1 (MQ=255)
                                                                                 gAAAGTTTATTAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAAttt           >  1:1307279/1‑91 (MQ=255)
                                                                                 gAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAAttt           <  1:3842181/91‑1 (MQ=255)
                                                                                   aaGTTTATTAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAg         >  2:2821781/1‑91 (MQ=255)
                                                                                     gTTTATTAAACTTTC‑TTTTTTAATATGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAGtt       >  2:2917586/1‑91 (MQ=255)
                                                                                     gTTTATTAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCACATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTCAGtt       >  2:3735566/1‑91 (MQ=255)
                                                                                     gTTTATTAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGCATTAATACATATTGTTTCACCAAATTTAGtt       >  2:174038/1‑91 (MQ=255)
                                                                                     gTTTATTAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCATATTTAGtt       >  2:3431548/1‑91 (MQ=255)
                                                                                     gTTTATTAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAGtt       >  1:843325/1‑91 (MQ=255)
                                                                                     gTTTATTAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAGtt       >  2:3167875/1‑91 (MQ=255)
                                                                                     gTTTATTAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAGtt       >  2:1006490/1‑91 (MQ=255)
                                                                                     gTTTATTAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAGtt       >  2:2326375/1‑91 (MQ=255)
                                                                                     gTTTATTAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAGtt       >  1:597046/1‑91 (MQ=255)
                                                                                      tttATTAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAGttt      >  2:1081993/1‑91 (MQ=255)
                                                                                      tttATTAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAGttt      >  1:1262123/1‑91 (MQ=255)
                                                                                       ttattaAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAACTTCGCGTAATTTAGAATTAATACATATTGTTTCACCAGATTTAGtttt     >  2:206299/1‑91 (MQ=255)
                                                                                       ttAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAGtttt     >  2:3808347/1‑91 (MQ=255)
                                                                                        taataaACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAGTTTTc    <  2:1528033/91‑1 (MQ=255)
                                                                                        taataaACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAGTTTTc    >  1:3428117/1‑91 (MQ=255)
                                                                                        taataaACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAGTTTTc    >  2:2140043/1‑91 (MQ=255)
                                                                                        taataaACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAGTTTTc    >  2:260679/1‑91 (MQ=255)
                                                                                         aataaACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAGTTTTcc   <  2:3348061/91‑1 (MQ=255)
                                                                                          |                                                                                           
TGTACCACCATTCTTGTCTATCTCGAAGAGAAATTTTTTCTAGATATATTTTCCATCCAATAGACAATTTTTTAGTATTAAGAAAGTTTAATAAACTTTC‑TTTTTTAATAAGTAAACATTCATTATCAAATTCGAGTAATTTAGAATTAATACATATTGTTTCACCAAATTTAGTTTTCCA  >  USA300TCH1516_ALE/1931509‑1931689

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