Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,218,730 G→T 11.4% intergenic (+221/+42) 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,7300GT11.4% 347.4 / 27.6 106intergenic (+221/+42)USA300HOU_RS06065/USA300HOU_RS06070hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (43/50);  new base T (5/7);  total (48/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

TAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATA  >  USA300TCH1516_ALE/1218640‑1218819
                                                                                          |                                                                                          
tAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAgg                                                                                            <  1:2377804/91‑1 (MQ=255)
 attttAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGa                                                                                           <  1:543140/88‑1 (MQ=255)
 aaCTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGa                                                                                           <  1:152962/91‑1 (MQ=255)
  aCTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGac                                                                                          <  1:3892062/91‑1 (MQ=255)
    ttAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGacta                                                                                        >  1:2207049/1‑91 (MQ=255)
     tAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGactac                                                                                       >  1:2093287/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:1210362/1‑91 (MQ=255)
         tgatgaTACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttc                                                                                   >  1:2047458/1‑91 (MQ=255)
         tgatgaTACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttc                                                                                   <  1:421807/91‑1 (MQ=255)
          gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                                 <  1:1585444/91‑1 (MQ=25)
          gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                                 <  1:1601879/91‑1 (MQ=25)
          gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATATAAATAAACACCTCACAAAAGAAAGTAGTCctttctta                                                                                <  1:663742/91‑2 (MQ=21)
            tgatACGTTATTTAAGAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttcttt                                                                               <  1:3426939/91‑1 (MQ=255)
             gatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCtttt                                                                              >  1:4019707/1‑91 (MQ=255)
             gatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCtttt                                                                              >  1:1824839/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)
                    tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgt                                                                       >  1:3077971/1‑91 (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)
                            aaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAttt                                                               <  1:3499544/91‑1 (MQ=255)
                            aaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAttt                                                               >  1:1127949/1‑91 (MQ=255)
                            aaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAttt                                                               >  1:3505261/1‑91 (MQ=255)
                            aaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAttt                                                               >  1:1532662/1‑91 (MQ=255)
                             acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCGCACAAAAGAAAGTAGGCCTTTCTTATGGTGAGGTGTTTAttta                                                              <  1:1734314/91‑2 (MQ=9)
                               acGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTc                                                            >  1:1977915/1‑91 (MQ=255)
                                cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                           <  1:2098546/91‑1 (MQ=255)
                                cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                           >  1:373790/1‑91 (MQ=255)
                                 gATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGa                                                          <  1:1231524/91‑1 (MQ=255)
                                    tcTTCGTTTCGATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACaa                                                       <  1:2457997/90‑1 (MQ=255)
                                    aCTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACaa                                                       <  1:2458255/91‑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)
                                                  aaCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTg                                         <  1:3470896/91‑1 (MQ=255)
                                                  aaCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTg                                         <  1:2064164/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:4032216/1‑91 (MQ=255)
                                                      cATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTc                                     <  1:123855/91‑1 (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:3310153/91‑1 (MQ=255)
                                                          aaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttt                                 >  1:3359560/1‑91 (MQ=255)
                                                              aataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTg                             <  1:2646079/91‑1 (MQ=255)
                                                              aataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTg                             <  1:4101985/91‑1 (MQ=255)
                                                               ataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                            <  1:2529880/91‑1 (MQ=255)
                                                               ataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                            <  1:1028794/91‑1 (MQ=255)
                                                                taaataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGaa                           >  1:4050174/1‑91 (MQ=255)
                                                                  aataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATa                         <  1:1340590/91‑1 (MQ=255)
                                                                   ataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATaa                        <  1:2315104/91‑1 (MQ=255)
                                                                    taaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAg                       <  1:97389/91‑1 (MQ=255)
                                                                    taaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAg                       >  1:3971644/1‑91 (MQ=255)
                                                                    taaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAg                       <  1:2447652/91‑1 (MQ=255)
                                                                    taaaCACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAg                       >  1:3973087/1‑91 (MQ=14)
                                                                      aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt                     <  1:1544599/91‑1 (MQ=255)
                                                                      aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt                     <  1:1039404/91‑1 (MQ=255)
                                                                      aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt                     <  1:3588558/91‑1 (MQ=255)
                                                                      aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt                     <  1:3022825/91‑1 (MQ=255)
                                                                       acacCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTc                    <  1:444594/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:1643653/1‑91 (MQ=255)
                                                                         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)
                                                                                     gAAAGCACTACTTTCTTTTGTGAGGAGTTTATTTGTCGACAACTTGATTATGAGTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAg      <  1:3857051/91‑1 (MQ=255)
                                                                                      aaaGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGc     >  1:1093531/1‑91 (MQ=255)
                                                                                         ggACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATa  >  1:3161892/1‑91 (MQ=255)
                                                                                         ggACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATa  >  1:1590488/1‑91 (MQ=255)
                                                                                         ggACTACTTTCTGTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATa  >  1:3185388/1‑91 (MQ=255)
                                                                                          |                                                                                          
TAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATA  >  USA300TCH1516_ALE/1218640‑1218819

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: