Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,218,734 A→C 11.5% intergenic (+225/+38) USA300HOU_RS06065 → / ← USA300HOU_RS06070 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,218,7340AC11.5% 332.8 / 30.8 104intergenic (+225/+38)USA300HOU_RS06065/USA300HOU_RS06070hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (44/48);  new base C (5/7);  total (49/55)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 109 total aligned reads displayed.

TTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAA  >  USA300TCH1516_ALE/1218644‑1218822
                                                                                          |                                                                                         
ttAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGacta                                                                                           >  1:2207049/1‑91 (MQ=255)
   aaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGTATTAAATAATAAATAAACACCTCACCATACGAAAGGACTActt                                                                                        >  1:4147222/1‑91 (MQ=255)
   aaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActt                                                                                        >  1:2717950/1‑91 (MQ=255)
   aaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActt                                                                                        >  1:3826636/1‑91 (MQ=255)
     tgatgaTACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttc                                                                                      >  1:2047458/1‑91 (MQ=255)
     tgatgaTACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttc                                                                                      >  1:1210362/1‑91 (MQ=255)
      gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                                    <  1:1601879/91‑1 (MQ=25)
      gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                                    <  1:1585471/91‑1 (MQ=25)
      gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                                    <  1:3947360/91‑1 (MQ=25)
      gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                                    <  1:1585444/91‑1 (MQ=25)
      gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATATAAATAAACACCTCACAAAAGAAAGTAGTCctttctta                                                                                   <  1:663742/91‑2 (MQ=21)
        tgatACGTTATTTAAGAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttcttt                                                                                  <  1:3426939/91‑1 (MQ=255)
         gatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCtttt                                                                                 >  1:1824839/1‑91 (MQ=255)
         gatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCtttt                                                                                 >  1:4019707/1‑91 (MQ=255)
         gatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttcttatg                                                                                 <  1:1884895/91‑4 (MQ=21)
          atACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTtg                                                                                <  1:694064/91‑1 (MQ=255)
           tACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTtgt                                                                               <  1:985161/91‑1 (MQ=255)
           tACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTtgt                                                                               >  1:3665857/1‑91 (MQ=255)
               ttATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAgg                                                                           <  1:3816301/91‑1 (MQ=255)
                tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgt                                                                          <  1:2599431/91‑1 (MQ=255)
                tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTCTTGTGAGgt                                                                          >  1:3962244/1‑91 (MQ=255)
                 aTTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtg                                                                         <  1:3770732/91‑1 (MQ=255)
                  tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                                        <  1:648349/91‑1 (MQ=255)
                  tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                                        <  1:2254227/91‑1 (MQ=255)
                  tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                                        <  1:2254066/91‑1 (MQ=255)
                       aaaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAtt                                                                   >  1:201982/1‑91 (MQ=255)
                        aaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAttt                                                                  >  1:1532662/1‑91 (MQ=255)
                        aaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAttt                                                                  <  1:3499544/91‑1 (MQ=255)
                        aaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAttt                                                                  >  1:1127949/1‑91 (MQ=255)
                         acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTg                                                                 <  1:2737147/91‑1 (MQ=255)
                         acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCGCACAAAAGAAAGTAGGCCTTTCTTATGGTGAGGTGTTTAttta                                                                 <  1:1734314/91‑2 (MQ=9)
                           acGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTc                                                               >  1:1977915/1‑91 (MQ=255)
                            cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                              >  1:373790/1‑91 (MQ=255)
                            cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                              <  1:2098546/91‑1 (MQ=255)
                             gATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGa                                                             <  1:1231524/91‑1 (MQ=255)
                                tcTTCGTTTCGATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACaa                                                          <  1:2457997/90‑1 (MQ=255)
                                   tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                                       >  1:2174579/1‑91 (MQ=255)
                                   tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                                       >  1:1343173/1‑91 (MQ=255)
                                   tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                                       >  1:593848/1‑91 (MQ=255)
                                   tCGATTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                                       <  1:3732165/91‑1 (MQ=255)
                                    cGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACaatcat                                                      <  1:1713389/91‑5 (MQ=255)
                                           aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                               >  1:1879615/1‑91 (MQ=255)
                                           aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                               >  1:3203564/1‑91 (MQ=255)
                                              aaCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTg                                            <  1:2064164/91‑1 (MQ=255)
                                              aaCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTg                                            <  1:3470896/91‑1 (MQ=255)
                                               aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                           <  1:3898920/91‑1 (MQ=255)
                                               aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                           <  1:1320388/91‑1 (MQ=255)
                                                cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                          <  1:3542953/91‑1 (MQ=255)
                                                cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                          <  1:1011621/91‑1 (MQ=255)
                                                 gCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGttt                                         <  1:2694856/91‑1 (MQ=255)
                                                  cATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTc                                        <  1:123855/91‑1 (MQ=255)
                                                  cATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTc                                        >  1:4032216/1‑91 (MQ=255)
                                                   aTTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCa                                       >  1:3787701/1‑91 (MQ=255)
                                                    ttAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAt                                      <  1:3947531/91‑1 (MQ=255)
                                                     taaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtt                                     >  1:3094013/1‑91 (MQ=255)
                                                      aaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttt                                    >  1:1290563/1‑91 (MQ=255)
                                                      aaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttt                                    <  1:3310153/91‑1 (MQ=255)
                                                      aaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttt                                    >  1:3359560/1‑91 (MQ=255)
                                                          aataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTg                                <  1:4101985/91‑1 (MQ=255)
                                                          aataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTg                                <  1:2646079/91‑1 (MQ=255)
                                                           ataaatagACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                               >  1:3616726/1‑91 (MQ=255)
                                                           ataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                               <  1:1028794/91‑1 (MQ=255)
                                                           ataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                               <  1:2529880/91‑1 (MQ=255)
                                                            taaataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGaa                              >  1:4050174/1‑91 (MQ=255)
                                                               ataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATaa                           <  1:2315104/91‑1 (MQ=255)
                                                               ataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATaa                           <  1:2315416/91‑1 (MQ=255)
                                                                taaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAg                          >  1:3971644/1‑91 (MQ=255)
                                                                taaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAg                          <  1:97389/91‑1 (MQ=255)
                                                                taaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAg                          <  1:2447652/91‑1 (MQ=255)
                                                                taaaCACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAg                          >  1:3973087/1‑91 (MQ=14)
                                                                  aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt                        <  1:3588558/91‑1 (MQ=255)
                                                                  aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt                        <  1:1039404/91‑1 (MQ=255)
                                                                  aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt                        <  1:1544599/91‑1 (MQ=255)
                                                                  aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt                        <  1:3022825/91‑1 (MQ=255)
                                                                   acacCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTc                       <  1:444594/91‑1 (MQ=255)
                                                                   acacCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTc                       <  1:444593/91‑1 (MQ=255)
                                                                    cacCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCa                      >  1:2754694/1‑91 (MQ=255)
                                                                    cacCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCa                      >  1:2754869/1‑91 (MQ=255)
                                                                    cacCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCa                      >  1:1492270/1‑91 (MQ=17)
                                                                     acCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCAt                     >  1:1643654/1‑91 (MQ=255)
                                                                       cTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATaa                   >  1:34598/1‑91 (MQ=255)
                                                                        tCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAAt                  >  1:102284/1‑91 (MQ=255)
                                                                         cACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATGAGTTCATAATc                 <  1:659202/91‑1 (MQ=255)
                                                                         cACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATc                 <  1:2890507/91‑1 (MQ=255)
                                                                            caaaaGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTCTGATTGTTTCATTTTTTGAATAAGTTCATAATCAg               >  1:2335172/4‑91 (MQ=18)
                                                                            cATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAgg              >  1:1513701/1‑91 (MQ=255)
                                                                            cATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAgg              >  1:1166397/1‑91 (MQ=255)
                                                                             aTAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgt             >  1:2021047/1‑91 (MQ=255)
                                                                              tAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtg            >  1:3148655/1‑91 (MQ=255)
                                                                              aaaGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgt             >  1:925577/2‑91 (MQ=21)
                                                                               aaGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtg            >  1:1405002/1‑91 (MQ=21)
                                                                                 gAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAg         >  1:3000395/1‑91 (MQ=255)
                                                                                 gAAAGGACTACTTTCTTTTGTGAGGTGTTTATATGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAg         <  1:3535507/91‑1 (MQ=255)
                                                                                  aaaGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGc        >  1:1093531/1‑91 (MQ=255)
                                                                                     ggACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATa     >  1:1590488/1‑91 (MQ=255)
                                                                                     ggACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATa     >  1:3161892/1‑91 (MQ=255)
                                                                                     ggACTACTTTCTGTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATa     >  1:3185388/1‑91 (MQ=255)
                                                                                       actactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACa   >  1:925149/1‑91 (MQ=255)
                                                                                       actactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACa   <  1:1564050/91‑1 (MQ=255)
                                                                                        ctactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACaa  <  1:2377352/91‑1 (MQ=255)
                                                                                          |                                                                                         
TTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAA  >  USA300TCH1516_ALE/1218644‑1218822

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: