Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 659,371 T→A 5.8% intergenic (+177/‑735) USA300TCH1516_00590 → / → USA300HOU_RS03165 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE659,3710TA 5.8% 900.3 / 28.4 267intergenic (+177/‑735)USA300TCH1516_00590/USA300HOU_RS03165hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (136/114);  new base A (13/4);  total (149/118)
Fisher's exact test for biased strand distribution p-value = 8.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 318 total aligned reads displayed.

GTACGGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAGGA  >  USA300TCH1516_ALE/659281‑659461
                                                                                          |                                                                                          
gTACGGTAACTTCTGTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATtttttatttt                                                                                            <  1:2433561/91‑1 (MQ=255)
 tACGGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATtttttattttt                                                                                           >  2:3405928/1‑91 (MQ=255)
 tACGGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATtttttattttt                                                                                           >  1:1233403/1‑91 (MQ=255)
 tACGGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATtttttattttt                                                                                           >  1:1480390/1‑91 (MQ=255)
    ggTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATc                                                                                        >  1:2110038/1‑91 (MQ=255)
      tAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCta                                                                                      >  1:2567277/1‑91 (MQ=255)
         cttcttTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaaat                                                                                   <  1:2885626/91‑1 (MQ=255)
          ttcttTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATTTTAATAGGGATGTTAATTTTTTATTTTTATCtaaata                                                                                  <  1:3332370/91‑1 (MQ=255)
          ttcttTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaaata                                                                                  <  1:716559/91‑1 (MQ=255)
            cgttAAGGTTCAGATTGAAATACAACTAGCGGTCACTGTTGTTTTTATCATGTTAATAGGGATGTTAATTTGTTATTTTTATCgaaataaa                                                                                <  2:180894/89‑1 (MQ=255)
             tttAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATaaaa                                                                               >  2:1962179/1‑91 (MQ=255)
             tttAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATaaaa                                                                               <  1:1930179/91‑1 (MQ=255)
               tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATTTAAATtataaa                                                                             >  2:1308333/1‑91 (MQ=255)
               tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATaaaaaa                                                                             <  1:1296866/91‑1 (MQ=255)
               tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATaaaaaa                                                                             <  2:3706152/91‑1 (MQ=255)
                 ttGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGTGATGTTAATTTTTTATTTTTATCTAAATAAAAAAtt                                                                           <  2:556749/90‑1 (MQ=255)
                 attttCAGTTTGAATAACAACTAGCATTCATTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAAtt                                                                           <  2:2029808/88‑1 (MQ=255)
                  tGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTa                                                                          >  1:3496957/1‑91 (MQ=255)
                  tGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTa                                                                          >  1:1782119/1‑91 (MQ=255)
                    ttCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAAc                                                                        <  1:1551671/91‑1 (MQ=255)
                     tCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACa                                                                       >  2:3289306/1‑91 (MQ=255)
                     tCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACa                                                                       >  2:1915670/1‑91 (MQ=255)
                     tCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACa                                                                       >  2:3591936/1‑91 (MQ=255)
                      cAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAt                                                                      >  1:3634289/1‑91 (MQ=255)
                      cAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAt                                                                      >  1:2179338/1‑91 (MQ=255)
                       aGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                     <  1:41790/91‑1 (MQ=255)
                       aGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                     <  1:2611665/91‑1 (MQ=255)
                       aGTTTGAAAAACAACTAGCAATCACTGTATTTATTATCAGGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                     <  2:21693/91‑1 (MQ=255)
                        gTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTa                                                                    >  1:2347831/1‑91 (MQ=255)
                          ttGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaa                                                                  <  1:2219051/91‑1 (MQ=255)
                           tGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaaa                                                                 <  2:1572075/91‑1 (MQ=255)
                            gAAAAACAACTAGCAATCACTGTATTTTTTATCATGTCAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaaaa                                                                >  1:3344169/1‑91 (MQ=255)
                             aaaaaCAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAAt                                                               >  2:302575/1‑91 (MQ=255)
                             aaaaaCAACTAGCAATCACTGTATTTTTTATCATGTTAATAAGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAAt                                                               <  2:1037663/91‑1 (MQ=255)
                              aaaaCAACTAGCAATCACTGTATTTTTTATCATCTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAAtt                                                              <  1:3245744/91‑1 (MQ=255)
                                aacaacTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAAtttt                                                            <  2:2967965/91‑1 (MQ=255)
                                aacaacTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAAtttt                                                            <  1:1852174/91‑1 (MQ=255)
                                  caacTAGCAATCACTGTTTTTTTTATCATGTTTATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTa                                                          >  2:2445298/1‑91 (MQ=255)
                                  caacTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTa                                                          >  2:2325889/1‑91 (MQ=255)
                                    acTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAAt                                                        >  1:1586907/1‑91 (MQ=255)
                                     cTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATa                                                       <  1:1954418/91‑1 (MQ=255)
                                       aGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAcc                                                     <  1:2970647/91‑1 (MQ=255)
                                       aGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAcc                                                     <  1:667664/91‑1 (MQ=255)
                                           aTCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACg                                                 <  1:1116793/91‑1 (MQ=255)
                                            tCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGt                                                >  1:1042239/1‑91 (MQ=255)
                                            tCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGt                                                >  1:27586/1‑91 (MQ=255)
                                               cTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGt                                             >  2:3948518/1‑91 (MQ=255)
                                                  tATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTtaa                                          <  1:1181811/91‑1 (MQ=255)
                                                   aTTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTtaat                                         >  2:62608/1‑91 (MQ=255)
                                                   aTTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTtaat                                         <  2:3879795/91‑1 (MQ=255)
                                                   aTTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTtaat                                         <  2:1501968/91‑1 (MQ=255)
                                                   aTTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTtaat                                         >  2:2657902/1‑91 (MQ=255)
                                                   aTTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTtaat                                         >  2:3077305/1‑91 (MQ=255)
                                                    ttttttATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTCATACCAACGCGGTtaata                                        >  2:1223275/1‑91 (MQ=255)
                                                    ttttttATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTtaata                                        >  1:1568859/1‑91 (MQ=255)
                                                      ttttATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAg                                      <  2:2589003/91‑1 (MQ=255)
                                                        ttATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGaaa                                   >  2:3324697/1‑90 (MQ=255)
                                                         tATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAAt                                   >  2:817991/1‑91 (MQ=255)
                                                         tATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAAt                                   <  1:3858236/91‑1 (MQ=255)
                                                           tCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATag                                 <  1:1247875/91‑1 (MQ=255)
                                                           tCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATag                                 >  2:1666944/1‑91 (MQ=255)
                                                            cATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATaga                                >  2:3104969/1‑91 (MQ=255)
                                                            cATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATaga                                >  2:1581744/1‑91 (MQ=255)
                                                            cATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATaga                                >  2:3072980/1‑91 (MQ=255)
                                                                ttAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTg                            >  2:1664095/1‑91 (MQ=255)
                                                                ttAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTg                            >  1:2458566/1‑91 (MQ=255)
                                                                  aaTAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAGATTTTTAATACCACCGTGGTTAATAAGAATAGAGTTGTa                          <  2:3752360/91‑1 (MQ=255)
                                                                    tAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTAtt                        >  2:2412161/1‑91 (MQ=255)
                                                                    aaGGGATGTAAAGTTTTTATTTATATCTAAATAAAAAATTAAGATTAAAAATTTTTAATACCAGCGTGGTTAATAAGAATAGTGTTGTAtt                        <  2:498840/90‑1 (MQ=255)
                                                                       ggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAAt                     >  1:643503/1‑91 (MQ=255)
                                                                        gATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAGTTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATa                    >  1:3148970/1‑91 (MQ=255)
                                                                        aaTGTTAATTTTTTATTTAGATAAAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATa                    >  1:1671670/2‑91 (MQ=21)
                                                                         aTGTTAATTTTTTATTTTTTTTTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATAc                   >  2:1657581/1‑91 (MQ=255)
                                                                         aTGTTAATTTTTTATTTAGATAAAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATAc                   >  2:3911902/1‑91 (MQ=25)
                                                                         aTGTTAATTTTTTATTTAGATAAAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATAc                   >  1:3485075/1‑91 (MQ=25)
                                                                          tGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATAcc                  <  1:2038820/91‑1 (MQ=255)
                                                                           gTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCg                 >  1:2620216/1‑91 (MQ=255)
                                                                            ttAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGa                >  2:721482/1‑91 (MQ=255)
                                                                              aaTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACt              <  1:3023608/91‑1 (MQ=255)
                                                                               aTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACtt             <  2:2802205/91‑1 (MQ=255)
                                                                               aTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACtt             >  1:1425492/1‑91 (MQ=255)
                                                                                ttttttATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACttt            >  1:2287558/1‑91 (MQ=255)
                                                                                ttttttATTTAGATAAAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACttt            >  1:845616/1‑91 (MQ=37)
                                                                                   tttatttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCaa         <  1:2950523/91‑1 (MQ=255)
                                                                                    ttatttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCaaa        >  2:2543527/1‑91 (MQ=255)
                                                                                    ttatttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCaaa        >  1:3797602/1‑91 (MQ=255)
                                                                                    ttatttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCaaa        <  2:3169505/91‑1 (MQ=255)
                                                                                     tatttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAAc       >  2:1120022/1‑91 (MQ=255)
                                                                                     tatttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCAACTTTCAAAc       <  1:3770364/91‑1 (MQ=255)
                                                                                      atttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACa      >  2:2288390/1‑91 (MQ=255)
                                                                                      atttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACa      >  2:780506/1‑91 (MQ=255)
                                                                                       tttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACaa     >  2:894709/1‑91 (MQ=255)
                                                                                       tttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACaa     <  2:2183397/91‑1 (MQ=255)
                                                                                       tttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACaa     <  2:1686947/91‑1 (MQ=255)
                                                                                       tttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACaa     >  1:3351974/1‑91 (MQ=255)
                                                                                       tttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACaa     >  1:1525168/1‑91 (MQ=255)
                                                                                        ttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAg    >  2:2220368/1‑91 (MQ=255)
                                                                                        ttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAg    <  1:766128/91‑1 (MQ=255)
                                                                                        ttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAg    >  1:579192/1‑91 (MQ=255)
                                                                                          ttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAGGa  <  1:2833488/91‑1 (MQ=255)
                                                                                          |                                                                                          
GTACGGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAGGA  >  USA300TCH1516_ALE/659281‑659461

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: