Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,929,573 C→A 5.3% R31L (CGC→CTC)  USA300TCH1516_01793 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,929,5730CA 5.3% 989.8 / 28.3 282R31L (CGC→CTC) USA300TCH1516_01793hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (138/129);  new base A (7/8);  total (145/137)
Fisher's exact test for biased strand distribution p-value = 7.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 293 total aligned reads displayed.

GAATGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTCTCTTGTCA  >  USA300TCH1516_ALE/1929483‑1929663
                                                                                          |                                                                                          
 aaTGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTcgcg                                                                                           >  2:2368654/1‑91 (MQ=255)
 aaTGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTcgcg                                                                                           <  2:2980621/91‑1 (MQ=255)
  aTGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGt                                                                                          <  1:301563/91‑1 (MQ=255)
   tGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGta                                                                                         <  1:960968/91‑1 (MQ=255)
   tGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGta                                                                                         >  2:2659932/1‑91 (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:572086/91‑1 (MQ=255)
        ccGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAAtt                                                                                    <  1:1052478/91‑1 (MQ=255)
        ccGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAAtt                                                                                    <  1:856517/91‑1 (MQ=255)
         cGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAAttt                                                                                   <  1:1803755/91‑1 (MQ=255)
         cGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAAttt                                                                                   >  2:3296122/1‑91 (MQ=255)
          gtgtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATACCCGCCATTTACAAACGTCGAAGGGGTGAAATCATACTCGCGTATAATTTc                                                                                  <  2:1011517/91‑1 (MQ=255)
          gtgtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTc                                                                                  >  2:2239210/1‑91 (MQ=255)
          gtgtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTc                                                                                  <  2:1734129/91‑1 (MQ=255)
          gtgtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTc                                                                                  <  1:2762201/91‑1 (MQ=255)
          gtgtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTc                                                                                  >  1:3699583/1‑91 (MQ=255)
            gtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATTATACGCGAGTATGATTTCac                                                                                <  2:2903041/91‑2 (MQ=38)
              tttGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATTATACGCGAGTATGATTTCacac                                                                              <  1:954253/91‑4 (MQ=37)
                  ttGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAg                                                                          <  2:576260/90‑1 (MQ=255)
                  aTGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAg                                                                          >  1:1301338/1‑91 (MQ=255)
                   tGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAgg                                                                         >  2:3715806/1‑91 (MQ=255)
                    gTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGc                                                                        <  1:1685780/91‑1 (MQ=255)
                     ttGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCt                                                                       <  1:3290304/91‑1 (MQ=255)
                     ttGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCt                                                                       >  2:2716136/1‑91 (MQ=255)
                      tGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCtt                                                                      <  2:1362486/91‑1 (MQ=255)
                      tGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCtt                                                                      <  2:2469584/91‑1 (MQ=255)
                      tGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCtt                                                                      >  2:2692626/1‑91 (MQ=255)
                      tGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCtt                                                                      <  1:1883552/91‑1 (MQ=255)
                       gCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTa                                                                     >  2:1472799/1‑91 (MQ=255)
                        cTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTAc                                                                    <  2:3125412/91‑1 (MQ=255)
                        cTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTAc                                                                    <  1:3773839/91‑1 (MQ=255)
                         tttATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTAcc                                                                   >  1:3207333/1‑91 (MQ=255)
                         tttATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTAcc                                                                   <  1:1004179/91‑1 (MQ=255)
                          ttATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCa                                                                  >  1:1051673/1‑91 (MQ=255)
                            aTCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAtt                                                                <  2:791756/91‑1 (MQ=255)
                             tCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAttt                                                               >  2:3032788/1‑91 (MQ=255)
                             tCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAttt                                                               >  2:3015868/1‑91 (MQ=255)
                              ccATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAtttt                                                              <  1:1409420/91‑1 (MQ=255)
                                 ttttCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTtat                                                           >  1:2195575/1‑91 (MQ=255)
                                 ttttCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTtat                                                           >  2:3293795/1‑91 (MQ=255)
                                 ttttCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTtat                                                           >  1:2902388/1‑91 (MQ=255)
                                  tttCCTACATTTTATAACCGCCATTTACAAACGGCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTtata                                                          <  1:3230906/91‑1 (MQ=255)
                                    tCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTtatata                                                        <  1:2163262/91‑1 (MQ=255)
                                     ccTACATTTTATAACCGCCATTTACAAACTTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATaa                                                       <  2:2548907/91‑1 (MQ=255)
                                     ccTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATaa                                                       >  1:1602930/1‑91 (MQ=255)
                                     ccTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATaa                                                       <  2:1997600/91‑1 (MQ=255)
                                     ccTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATaa                                                       <  1:3515184/91‑1 (MQ=255)
                                      cTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAAt                                                      <  1:1442593/91‑1 (MQ=255)
                                        aCATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCt                                                    >  1:1917667/1‑91 (MQ=255)
                                          aTTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTaa                                                  <  1:1371795/91‑1 (MQ=255)
                                           ttttATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAAc                                                 <  1:2692626/91‑1 (MQ=255)
                                            tttATAACCGCCATTTACAAACGTCGAAGGTATGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAAcc                                                >  1:315205/1‑91 (MQ=255)
                                              tataACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAt                                              >  1:2363582/1‑91 (MQ=255)
                                              tataACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAt                                              <  2:61383/91‑1 (MQ=255)
                                                 aaCCGCCATTTACAAACGTCGAAGGTGTGACATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAGTTg                                           >  1:2370878/1‑91 (MQ=255)
                                                 aaCCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTg                                           <  1:2934840/91‑1 (MQ=255)
                                                 aaCCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTg                                           >  1:1658687/1‑91 (MQ=255)
                                                  aCCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGt                                          <  1:3848303/91‑1 (MQ=255)
                                                     gCCATTTACAAACGTCGCAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAAc                                       <  1:3329109/91‑1 (MQ=255)
                                                        aTTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTa                                    <  1:1258096/91‑1 (MQ=255)
                                                           tACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAAc                                 >  1:2275037/1‑91 (MQ=255)
                                                             cAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAActc                               <  2:3548829/91‑1 (MQ=255)
                                                                aCGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGa                            >  2:2293538/1‑91 (MQ=255)
                                                                 cGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGaa                           >  1:833717/1‑91 (MQ=255)
                                                                 cGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGaa                           <  2:1388283/91‑1 (MQ=255)
                                                                  gTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAAc                          <  2:961944/91‑1 (MQ=255)
                                                                    ctaaGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTa                        <  2:1655640/89‑1 (MQ=255)
                                                                    cGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTa                        <  2:2043409/91‑1 (MQ=255)
                                                                     gAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTaa                       >  1:2850180/1‑91 (MQ=255)
                                                                       aGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAAt                     >  1:3652312/1‑91 (MQ=255)
                                                                        ggTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTTAACTAAATg                    <  1:1954801/91‑1 (MQ=255)
                                                                        ggTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATg                    >  2:1851959/1‑91 (MQ=255)
                                                                        ggTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATg                    <  2:3153009/91‑1 (MQ=255)
                                                                           gtgAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtc                 >  1:1897958/1‑91 (MQ=255)
                                                                            tgAAATTATACGCGAGTATGATTTCACACCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtct                >  1:2391131/1‑91 (MQ=14)
                                                                            tgAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTaaa                      <  2:1296293/85‑1 (MQ=255)
                                                                             gAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAAc                                       <  2:955859/67‑1 (MQ=255)
                                                                             gAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtctt               <  2:2421245/91‑1 (MQ=255)
                                                                             gAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtctt               <  1:3627729/91‑1 (MQ=255)
                                                                              aaaTTATACGCGAGTATGATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttc              >  2:397414/1‑91 (MQ=14)
                                                                              aaaTTATACGCGAGTATGATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttc              >  2:3496069/1‑91 (MQ=14)
                                                                              aaaTTATACGCGAGTATGATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttc              >  2:2839773/1‑91 (MQ=14)
                                                                              aaaTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttc              <  1:1479528/91‑1 (MQ=255)
                                                                              aaaTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttc              >  1:2037020/1‑91 (MQ=255)
                                                                              aaaTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttc              <  2:902940/91‑1 (MQ=255)
                                                                              aaaTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttc              <  2:3806755/91‑1 (MQ=255)
                                                                               aaTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttct             <  1:3534114/91‑1 (MQ=255)
                                                                               aaTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttct             <  1:678825/91‑1 (MQ=255)
                                                                                aTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATcttctt            >  2:2488497/1‑91 (MQ=255)
                                                                                 tCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttcttt           <  2:528281/91‑1 (MQ=255)
                                                                                  cATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttctttc          >  2:2119723/1‑91 (MQ=255)
                                                                                  cATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttctttc          >  1:3075477/1‑91 (MQ=255)
                                                                                    tACTCGCGTATGATTTCATTCCTCGGCTTACCATTTTTATATAATCTAAGCATTTGTAAGTTAAACTCTGAACTAAATGATCTTCTTtctc        <  2:2955984/91‑1 (MQ=17)
                                                                                    tACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTtctc        >  2:1084803/1‑91 (MQ=255)
                                                                                     aCTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTctct       >  1:2118436/1‑91 (MQ=255)
                                                                                     aCTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTctct       >  1:435201/1‑91 (MQ=255)
                                                                                       tCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTCTCTTg     >  2:1750886/1‑91 (MQ=255)
                                                                                         gcgTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTCTCTTGTc   <  1:3346579/91‑1 (MQ=255)
                                                                                          |                                                                                          
GAATGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTCTCTTGTCA  >  USA300TCH1516_ALE/1929483‑1929663

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