Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 265,622 T→G 6.4% intergenic (+22/+142) ugpQ_2 → / ← scn_1 Glycerophosphodiester phosphodiesterase, cytoplasmic/Staphylococcal complement inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE265,6220TG 6.4% 654.4 / 25.9 196intergenic (+22/+142)ugpQ_2/scn_1Glycerophosphodiester phosphodiesterase, cytoplasmic/Staphylococcal complement inhibitor
Reads supporting (aligned to +/- strand):  ref base T (82/100);  new base G (9/4);  total (92/104)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 204 total aligned reads displayed.

GATCAGATTAAAGAAGAAAAGAAAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTAT  >  USA300TCH1516_ALE/265532‑265712
                                                                                          |                                                                                          
ggtCAGATTAAAGAAGAAAAGAAAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGt                                                                                            >  1:354729/3‑91 (MQ=255)
gATCAGATTAAAGAAGAAAAGAAAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGt                                                                                            <  2:3319334/91‑1 (MQ=255)
 aTCAGATTAAAGAAGAAAAGAAAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGtt                                                                                           >  2:2530737/1‑91 (MQ=255)
          aagaagAAAAGAAAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAgatg                                                                                  >  2:2641975/1‑91 (MQ=255)
           agaagaAAAGAAAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAgatga                                                                                 <  1:1610753/91‑1 (MQ=255)
             aagaaaagaaAGACGAAACATAGCTCGATCGTTCTATAAGAATTTTGTTGGAATAATATTAACAAAGACCTCTAAAGTTATCAAGATGATa                                                                               <  2:2954137/91‑1 (MQ=255)
             aagaaaagaaAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATa                                                                               <  1:2628814/91‑1 (MQ=255)
             aagaaaagaaAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATa                                                                               <  2:1628735/91‑1 (MQ=255)
              agaaaagaaaGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATAc                                                                              <  1:2691/91‑1 (MQ=255)
                aaaagaaaGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCt                                                                            >  2:3425475/1‑91 (MQ=255)
                 aaagaaagACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCtt                                                                           >  2:1709454/1‑91 (MQ=255)
                 aaagaaagACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCtt                                                                           <  2:2158864/91‑1 (MQ=255)
                 aaagaaagACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCtt                                                                           >  2:546589/1‑91 (MQ=255)
                    gaaagaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCag                                                                        <  2:3722809/91‑1 (MQ=255)
                    gaaagaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCag                                                                        <  1:567828/91‑1 (MQ=255)
                    gaaagaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCag                                                                        <  2:758424/91‑1 (MQ=255)
                     aaagaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCaga                                                                       >  1:432866/1‑91 (MQ=255)
                      aagaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCagag                                                                      >  1:798558/1‑91 (MQ=255)
                       agaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAgg                                                                     <  1:3703566/91‑1 (MQ=255)
                       agaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAgg                                                                     <  1:2333650/91‑1 (MQ=255)
                        gaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGt                                                                    >  1:3262299/1‑91 (MQ=255)
                        gaCGAAACATACGTCGATCGTTCTATAAGATTTTTGTTTGAATAATATAAACAAAGACCTCTGAAGGTATCATCTTGATAACTTTAGAGGt                                                                    <  2:1766055/91‑1 (MQ=18)
                         aCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTc                                                                   >  2:1938856/1‑91 (MQ=255)
                          cGAAACATACTTCGATCGTTCTATAAGAATGTTGTGTGAATAATATAAACAAAGACCTCTAAAGGTATCATCTTGATAACTTTAGAGGTCt                                                                  <  2:2828277/91‑1 (MQ=18)
                           gAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCtt                                                                 <  2:201724/91‑1 (MQ=255)
                           gAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCtt                                                                 <  1:3558522/91‑1 (MQ=255)
                                aTACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTGAAGGTATCATCTTGATAACTTTAGAGGTCtttgttt                                                            <  2:2473489/91‑5 (MQ=14)
                                aTACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTGAAGGTATCATCTTGATAACTTTAGAGGTCtttgttt                                                            <  2:3868254/91‑5 (MQ=17)
                                aTACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTaa                                                            <  2:1477615/91‑1 (MQ=255)
                                aTACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTaa                                                            <  1:65441/91‑1 (MQ=255)
                                  aCTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATg                                                          >  1:1382845/1‑91 (MQ=255)
                                  aCTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATg                                                          >  2:931112/1‑91 (MQ=255)
                                     tCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTg                                                       <  2:3763880/91‑1 (MQ=255)
                                     tCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTg                                                       <  2:93503/91‑1 (MQ=255)
                                      cGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGttac                                                      <  2:453932/91‑3 (MQ=255)
                                      cGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGc                                                      <  2:337256/91‑1 (MQ=255)
                                       gATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGcc                                                     >  1:198232/1‑91 (MQ=255)
                                       gATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAACACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGcc                                                     >  2:3193163/1‑91 (MQ=255)
                                         tCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCAt                                                   >  2:1351898/1‑91 (MQ=255)
                                           gTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGcctatc                                                 <  1:1073024/91‑5 (MQ=255)
                                           gTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCt                                                 >  2:1837544/1‑91 (MQ=255)
                                           gTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCt                                                 >  2:1333184/1‑91 (MQ=255)
                                            ttCTATAAGAATTTTGTTTGAATAATATAAACAAAGGCCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTa                                                <  2:3058870/91‑1 (MQ=255)
                                                ataAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATggg                                            >  2:1226507/1‑91 (MQ=255)
                                                  aaGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGAt                                          <  1:1578490/91‑1 (MQ=255)
                                                   aGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATa                                         <  2:2278186/91‑1 (MQ=255)
                                                     gaTTTTGTTTGAATAATATAAACAAAGTCCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAgg                                       <  1:2258785/90‑1 (MQ=255)
                                                     aaTTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCACGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAgg                                       <  1:3868024/91‑1 (MQ=255)
                                                     aaTTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAgg                                       <  1:1031212/91‑1 (MQ=255)
                                                     aaTTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAgg                                       <  1:2235735/91‑1 (MQ=255)
                                                     aaTTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAgg                                       >  2:2422668/1‑91 (MQ=255)
                                                     aaTTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAgg                                       <  2:1919181/91‑1 (MQ=255)
                                                     aaTTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAgg                                       >  1:714187/1‑91 (MQ=255)
                                                       ttttGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCa                                     <  2:3541704/91‑1 (MQ=255)
                                                       ttttGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCa                                     >  2:734445/1‑91 (MQ=255)
                                                        tttgtttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCaa                                    <  2:2808322/91‑1 (MQ=255)
                                                         ttgtttgAATAATATAGACAAAGACCTCGAAAGTTATCAGGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAAt                                   <  1:915408/91‑1 (MQ=255)
                                                           gtttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCg                                 >  1:2374599/1‑91 (MQ=255)
                                                            tttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTCATGTTGCCATCTATGGGATAGGCAATCGt                                >  1:3106673/1‑91 (MQ=255)
                                                            tttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGt                                >  2:1778335/1‑91 (MQ=255)
                                                            tttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGt                                >  2:2005971/1‑91 (MQ=255)
                                                             ttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGtt                               <  2:798063/91‑1 (MQ=255)
                                                             ttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGtt                               >  1:1213376/1‑91 (MQ=255)
                                                                 ataataTAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCAt                           <  2:1017139/91‑1 (MQ=255)
                                                                  taataTAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCAtt                          <  1:1564843/91‑1 (MQ=255)
                                                                  taataTAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCAtt                          <  1:2816926/91‑1 (MQ=255)
                                                                  taataTAAACAAAGACCTCTAAAGTTATCAAGATGACACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCAtt                          <  1:2085695/91‑1 (MQ=255)
                                                                     tataAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGt                       >  1:639002/1‑91 (MQ=255)
                                                                     tataAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGt                       >  1:1087140/1‑91 (MQ=255)
                                                                      ataAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGtt                      <  2:432866/91‑1 (MQ=255)
                                                                      ataAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGtt                      <  1:266536/91‑1 (MQ=255)
                                                                       taAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGttt                     >  1:1836860/1‑91 (MQ=255)
                                                                         aaCAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCACTCGTTTCATTCGTTtat                   >  1:3320769/1‑91 (MQ=255)
                                                                         aaCAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTtat                   >  2:1405647/1‑91 (MQ=255)
                                                                         aaCAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTtat                   >  1:1395546/1‑91 (MQ=255)
                                                                         aaCAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTtat                   >  1:2589778/1‑91 (MQ=255)
                                                                          aCAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTtata                  >  1:2240729/1‑91 (MQ=255)
                                                                           cAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTatat                 <  2:3023634/91‑1 (MQ=255)
                                                                            aaaGACCTCTGAAGGTATCATCTTGATAACTTTAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATAtt                >  2:914308/1‑91 (MQ=17)
                                                                             aaGACCTCTGAAGGTATCATCTTGATAACTTTAGAGGACTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTc               >  1:3270200/1‑91 (MQ=18)
                                                                              aGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCa              <  1:1930585/91‑1 (MQ=255)
                                                                              aGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCa              <  1:1448989/91‑1 (MQ=255)
                                                                               gACCTCTGAAGGTATCATCTTGATAACTTTAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCat             >  1:1701812/1‑91 (MQ=21)
                                                                               gACCTCTGAAGGTATCATCTTGATAACTTTAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCat             >  2:488524/1‑91 (MQ=21)
                                                                                aCCTCTGAAGGTATCATCTTGATAACTTTAGAGGTCTTTTTAATGTTGCCATATATGGGATAGGCAATCGTTTCATTAGTTTATATTAata            >  1:3832413/1‑91 (MQ=18)
                                                                                aCCTCTGAAGGTATCATCTTGATAACTTTAGAGGACTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCata            >  2:358108/1‑91 (MQ=21)
                                                                                aCCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATAttattc            >  2:2141065/1‑87 (MQ=255)
                                                                                 ccTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCatat           <  1:182490/91‑1 (MQ=255)
                                                                                 ccTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCatat           <  2:1824615/91‑1 (MQ=255)
                                                                                  ctctAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATg          <  1:3425475/91‑1 (MQ=255)
                                                                                   tctAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGa         >  2:3462280/1‑91 (MQ=255)
                                                                                   tctAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGa         >  1:1910755/1‑91 (MQ=255)
                                                                                   tctAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGa         >  1:2622619/1‑91 (MQ=255)
                                                                                     tAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACa       >  2:1161423/1‑91 (MQ=255)
                                                                                      aaaGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACaa      >  1:1458088/1‑91 (MQ=255)
                                                                                       aaGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAg     <  1:2048213/91‑1 (MQ=255)
                                                                                       aaGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAg     <  1:286998/91‑1 (MQ=255)
                                                                                       aaGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAg     >  1:2917740/1‑91 (MQ=255)
                                                                                        aGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGt    <  1:1381510/91‑1 (MQ=255)
                                                                                         gTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTa   >  1:3380829/1‑91 (MQ=255)
                                                                                          |                                                                                          
GATCAGATTAAAGAAGAAAAGAAAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTAT  >  USA300TCH1516_ALE/265532‑265712

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: