Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,929,570 T→G 5.2% E32A (GAG→GCG)  USA300TCH1516_01793 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,929,5700TG 5.2% 968.7 / 33.8 287E32A (GAG→GCG) USA300TCH1516_01793hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (137/135);  new base G (7/8);  total (144/143)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 292 total aligned reads displayed.

ATTGAATGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTCTCTTG  >  USA300TCH1516_ALE/1929480‑1929660
                                                                                          |                                                                                          
aTTGAATGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACt                                                                                            <  1:115647/91‑1 (MQ=255)
  tGAATGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTcg                                                                                          <  2:1831995/91‑1 (MQ=255)
     aTGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGt                                                                                       <  1:1950791/91‑1 (MQ=255)
      tGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGta                                                                                      <  1:960968/91‑1 (MQ=255)
      tGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGta                                                                                      >  2:1077014/1‑91 (MQ=255)
      tGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGta                                                                                      >  2:2659932/1‑91 (MQ=255)
        tACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATTATACGCGAGtatg                                                                                    <  2:2366199/91‑2 (MQ=255)
         aCCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATaa                                                                                   <  1:49397/91‑1 (MQ=255)
          cccGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATTATACGCGAGtatgat                                                                                  <  1:917595/91‑4 (MQ=255)
          cccGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATTATACGCGAGtatgat                                                                                  <  2:3582527/91‑4 (MQ=255)
          cccGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAAt                                                                                  >  1:3652106/1‑91 (MQ=255)
          cccGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAAt                                                                                  <  1:572086/91‑1 (MQ=255)
           ccGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAAtt                                                                                 <  1:1052478/91‑1 (MQ=255)
            cGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAAttt                                                                                <  1:1551952/91‑1 (MQ=255)
            cGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAAttt                                                                                >  2:3296121/1‑91 (MQ=255)
             gtgtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTc                                                                               <  2:1010303/91‑1 (MQ=255)
             gtgtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTc                                                                               <  1:838699/91‑1 (MQ=255)
             gtgtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTc                                                                               >  1:1340113/1‑91 (MQ=255)
               gtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATTATACGCGAGTATGATTTCac                                                                             <  2:2903288/91‑2 (MQ=38)
                ttttGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCAtt                                                                            >  1:1704617/1‑91 (MQ=255)
                 tttGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATTATACGCGAGTATGATTTCacac                                                                           <  1:954253/91‑4 (MQ=37)
                 tttGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTc                                                                           >  1:3343462/1‑91 (MQ=255)
                  ttGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTcc                                                                          >  2:2089197/1‑91 (MQ=255)
                    gATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTa                                                                        >  1:369824/1‑91 (MQ=255)
                     ttGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAg                                                                       <  2:576260/90‑1 (MQ=255)
                     aTGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAg                                                                       >  1:1299975/1‑91 (MQ=255)
                      tGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAgg                                                                      >  2:3715806/1‑91 (MQ=255)
                      tGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAgg                                                                      >  1:2162342/1‑91 (MQ=255)
                        ttGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCt                                                                    <  1:3290751/91‑1 (MQ=255)
                         tGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCtt                                                                   >  2:2692626/1‑91 (MQ=255)
                          gCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTa                                                                  >  2:2030768/1‑91 (MQ=255)
                          gCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTa                                                                  <  2:3143773/91‑1 (MQ=255)
                          gCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTa                                                                  >  2:1472799/1‑91 (MQ=255)
                           cTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTAc                                                                 <  1:3773839/91‑1 (MQ=255)
                            tttATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTAcc                                                                <  1:1004179/91‑1 (MQ=255)
                             ttATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCa                                                               >  1:1051673/1‑91 (MQ=255)
                              tATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAt                                                              >  1:3196896/1‑91 (MQ=255)
                               aTCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAtt                                                             <  2:791756/91‑1 (MQ=255)
                                 ccATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAtttt                                                           <  1:1675842/91‑1 (MQ=255)
                                    ttttCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTtat                                                        >  1:2195272/1‑91 (MQ=255)
                                        ccTACATTTTATAACCGCCATTTACAAACTTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATaa                                                    <  2:2548907/91‑1 (MQ=255)
                                         cTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAAt                                                   >  1:2316668/1‑91 (MQ=255)
                                         cTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAAt                                                   >  1:1897103/1‑91 (MQ=255)
                                          tACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATc                                                  <  1:3402550/91‑1 (MQ=255)
                                           aCATTTTGTAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCt                                                 <  2:3453818/91‑1 (MQ=255)
                                           aCATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCt                                                 >  1:1917633/1‑91 (MQ=255)
                                           aCATTTTATAACCGCCATTTACAAACGTCGAAGCTGTGAAATCATACTCGCGTATAATTACCTTCCTAGGCTTACCATTTTTATATAATCt                                                 >  1:3105140/1‑91 (MQ=255)
                                                 tataACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAt                                           <  2:61383/91‑1 (MQ=255)
                                                  ataACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAtt                                          >  1:3220346/1‑91 (MQ=255)
                                                   taACCGCCATTTGCAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAttt                                         >  1:2989342/1‑91 (MQ=255)
                                                   taACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAttt                                         >  2:2841049/1‑91 (MQ=255)
                                                    aaCCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTg                                        <  1:909415/91‑1 (MQ=255)
                                                    aaCCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTg                                        <  1:3037180/91‑1 (MQ=255)
                                                     aCCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGt                                       <  1:3848303/91‑1 (MQ=255)
                                                      ccGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTa                                      >  2:1891552/1‑91 (MQ=255)
                                                       cGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTaa                                     <  2:3507543/91‑1 (MQ=255)
                                                         ccATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTCGTAACt                                   >  2:2816120/1‑91 (MQ=255)
                                                          cATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACtt                                  >  1:1383709/1‑91 (MQ=255)
                                                           aTTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTa                                 <  1:1635907/91‑1 (MQ=255)
                                                             ttACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCGTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTaaa                               <  2:2997391/91‑1 (MQ=255)
                                                             ttACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTaaa                               >  1:2889620/1‑91 (MQ=255)
                                                              tACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAAc                              >  1:2377504/1‑91 (MQ=255)
                                                                cAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAActc                            <  2:3548829/91‑1 (MQ=255)
                                                                 aaaCGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAActct                           <  1:2716136/91‑1 (MQ=255)
                                                                  aaCGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTg                          <  1:2798694/91‑1 (MQ=255)
                                                                  aaCGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTATTATATAATCTAACCATTTATACCTTAAACTCTg                          >  1:3518680/1‑91 (MQ=255)
                                                                   atggcGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGa                         >  2:2199890/5‑91 (MQ=255)
                                                                    cGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGaa                        <  2:1388283/91‑1 (MQ=255)
                                                                     gTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAAc                       <  2:961944/91‑1 (MQ=255)
                                                                       ctaaGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTa                     <  2:1655640/89‑1 (MQ=255)
                                                                             tgtgAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAt               <  1:2368654/91‑1 (MQ=255)
                                                                              gtgAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtc              <  1:1985163/91‑1 (MQ=255)
                                                                                gAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAAc                                    <  2:955859/67‑1 (MQ=255)
                                                                                gAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtctt            <  2:2421191/91‑1 (MQ=255)
                                                                                 aaaTTATACGCGAGTATGATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttc           >  2:397414/1‑91 (MQ=14)
                                                                                 aaaTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttc           <  2:902940/91‑1 (MQ=255)
                                                                                 aaaTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttc           <  1:1479528/91‑1 (MQ=255)
                                                                                 aaaTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttc           <  2:3806755/91‑1 (MQ=255)
                                                                                  aaTCATACTCGCGTATAATTTCATTTATTGGCTTACCATTTTCATATAATCTAACCATTTATAACTTAAACTCTGAACGAAATGTtcttct          >  2:3220157/1‑91 (MQ=14)
                                                                                  aaTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttct          <  1:678825/91‑1 (MQ=255)
                                                                                   aTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATcttctt         >  2:1736438/1‑91 (MQ=255)
                                                                                   aTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATcttctt         >  2:2488497/1‑91 (MQ=255)
                                                                                    tCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttcttt        <  2:528281/91‑1 (MQ=255)
                                                                                     cATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttctttc       >  1:2266324/1‑91 (MQ=255)
                                                                                     cATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttctttc       <  2:2283986/91‑1 (MQ=255)
                                                                                     cATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttctttc       <  2:2461336/91‑1 (MQ=255)
                                                                                     cATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttctttc       <  2:1351947/91‑1 (MQ=255)
                                                                                     cATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttctttc       >  2:1284832/1‑91 (MQ=255)
                                                                                      aTACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttctttct      <  2:1792952/91‑1 (MQ=255)
                                                                                       tACTCGCGTATGATTTCATTCCTGGGCTTACCATTTTTATATAATCTAAGCATTTGTAAGTTAAACTCTGAACTAAATGATCTTCTTtctc     <  2:2956056/91‑1 (MQ=17)
                                                                                       tACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTtctc     <  1:2470621/91‑1 (MQ=255)
                                                                                       tACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTtctc     >  1:1840416/1‑91 (MQ=255)
                                                                                        aCTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATTTAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTctct    <  1:3240398/91‑1 (MQ=255)
                                                                                        aCTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTctct    <  1:2919518/91‑1 (MQ=255)
                                                                                        aCTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTctct    <  1:1435140/91‑1 (MQ=255)
                                                                                        aCTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTctct    >  1:435201/1‑91 (MQ=255)
                                                                                        aCTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTctct    >  2:3032507/1‑91 (MQ=255)
                                                                                        aCTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTctct    >  1:2117420/1‑91 (MQ=255)
                                                                                        aCTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTctct    >  2:3033231/1‑91 (MQ=255)
                                                                                        aCTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACCCTGAACTAAATGATCTTCTTTctct    >  1:3600224/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATTGAATGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTCTCTTG  >  USA300TCH1516_ALE/1929480‑1929660

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