Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 265,629 G→C 6.9% intergenic (+29/+135) 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,6290GC 6.9% 655.5 / 25.7 189intergenic (+29/+135)ugpQ_2/scn_1Glycerophosphodiester phosphodiesterase, cytoplasmic/Staphylococcal complement inhibitor
Reads supporting (aligned to +/- strand):  ref base G (78/98);  new base C (9/4);  total (87/102)
Fisher's exact test for biased strand distribution p-value = 9.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 198 total aligned reads displayed.

TAAAGAAGAAAAGAAAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTATTTGTAT  >  USA300TCH1516_ALE/265540‑265718
                                                                                         |                                                                                         
tAAAGAAGAAAAGAAAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAga                                                                                          <  2:764310/91‑1 (MQ=255)
  aagaagAAAAGAAAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAgatg                                                                                        >  2:2641975/1‑91 (MQ=255)
   agaagaAAAGAAAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAgatga                                                                                       <  1:1610753/91‑1 (MQ=255)
    gaagaaAAGAAAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAgatgat                                                                                      >  2:931554/1‑91 (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)
       gaaaagaaaGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATAcc                                                                                   >  2:942749/1‑91 (MQ=255)
        aaaagaaaGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCt                                                                                  >  2:3425475/1‑91 (MQ=255)
         aaagaaagACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCtt                                                                                 <  2:2158864/91‑1 (MQ=255)
         aaagaaagACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCtt                                                                                 >  2:546589/1‑91 (MQ=255)
            gaaagaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCag                                                                              <  1:567828/91‑1 (MQ=255)
            gaaagaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCag                                                                              <  2:3722809/91‑1 (MQ=255)
             aaagaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCaga                                                                             >  1:432866/1‑91 (MQ=255)
              aagaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCagag                                                                            >  1:798558/1‑91 (MQ=255)
              aagaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCagag                                                                            <  1:1078631/91‑1 (MQ=255)
               agaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAgg                                                                           <  1:3703566/91‑1 (MQ=255)
               agaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAgg                                                                           <  1:2333650/91‑1 (MQ=255)
                 aCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTc                                                                         >  2:1938856/1‑91 (MQ=255)
                  cGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCt                                                                        >  2:1909156/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)
                       cATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGATCTCAAGATGATACCTTCAGAGGTCTTTTTa                                                                   >  1:3104184/1‑91 (MQ=255)
                        aTACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTGAAGGTATCATCTTGATAACTTTAGAGGTCtttgttt                                                                  <  2:3868254/91‑5 (MQ=17)
                        aTACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTaa                                                                  <  2:1477615/91‑1 (MQ=255)
                        aTACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTaa                                                                  <  1:65441/91‑1 (MQ=255)
                         tACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAAt                                                                 <  2:1813828/91‑1 (MQ=255)
                          aCTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATg                                                                >  1:1382845/1‑91 (MQ=255)
                             tCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTg                                                             <  2:3763880/91‑1 (MQ=255)
                              cGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGttac                                                            <  2:453932/91‑3 (MQ=255)
                              cGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGc                                                            <  2:337256/91‑1 (MQ=255)
                               gATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGcc                                                           >  1:1055905/1‑91 (MQ=255)
                               gATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGcc                                                           >  1:198232/1‑91 (MQ=255)
                               gATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAACACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGcc                                                           >  2:3193163/1‑91 (MQ=255)
                                 tCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCAt                                                         >  2:1351898/1‑91 (MQ=255)
                                   gTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCt                                                       >  2:1333184/1‑91 (MQ=255)
                                   gTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCt                                                       >  2:1836901/1‑91 (MQ=255)
                                    ttCTATAAGAATTTTGTTTGAATAATATAAACAAAGGCCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTa                                                      <  2:3058870/91‑1 (MQ=255)
                                    ttCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTa                                                      <  2:91923/91‑1 (MQ=255)
                                    ttCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTa                                                      <  1:3358872/91‑1 (MQ=255)
                                       tataAGAATTTTGTTTGAATAATATAAACAAAGACCTGTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATgg                                                   >  1:3262659/1‑91 (MQ=255)
                                        ataAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATggg                                                  >  2:1226507/1‑91 (MQ=255)
                                        ataAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATggg                                                  <  1:86469/91‑1 (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:2235735/91‑1 (MQ=255)
                                             aaTTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAgg                                             <  1:1031212/91‑1 (MQ=255)
                                             aaTTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAgg                                             >  1:714187/1‑91 (MQ=255)
                                              aTTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGc                                            <  1:3310268/91‑1 (MQ=255)
                                               ttttGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCa                                           >  2:734445/1‑91 (MQ=255)
                                               ttttGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCa                                           <  2:3541704/91‑1 (MQ=255)
                                                tttgtttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCaa                                          <  2:2808322/91‑1 (MQ=255)
                                                 ttgtttgAATAATATAGACAAAGACCTCGAAAGTTATCAGGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAAt                                         <  1:915408/91‑1 (MQ=255)
                                                  tgtttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATc                                        <  2:2713881/91‑1 (MQ=255)
                                                   gtttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCg                                       >  1:2374599/1‑91 (MQ=255)
                                                    tttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGt                                      >  2:2005971/1‑91 (MQ=255)
                                                     ttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGtt                                     >  1:3537528/1‑91 (MQ=255)
                                                     ttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGACTTTTTAATGTTGCCATCTATGGGATAGGCAATCGtt                                     <  1:3259736/91‑1 (MQ=255)
                                                      tgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGttt                                    <  1:1195554/91‑1 (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)
                                                           aataTAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCAtac                               <  2:2447456/91‑3 (MQ=255)
                                                             tataAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGt                             >  1:639002/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)
                                                                aaaCAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTta                          >  2:1768182/1‑91 (MQ=255)
                                                                 aaCAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTtat                         >  1:1395546/1‑91 (MQ=255)
                                                                 aaCAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTtat                         >  2:1405647/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)
                                                                   aaaaGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTatat                       <  1:3412226/90‑1 (MQ=255)
                                                                    aaaGACCTCTGAAGGTATCATCTTGATAACTTTAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATAtt                      >  1:1269514/1‑91 (MQ=17)
                                                                    aaaGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATAtt                      >  2:1745772/1‑91 (MQ=255)
                                                                      aGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCa                    <  1:1448989/91‑1 (MQ=255)
                                                                      aGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCa                    <  1:1930585/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)
                                                                           tctAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGa               >  1:2622619/1‑91 (MQ=255)
                                                                           tctAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGa               >  2:3462280/1‑91 (MQ=255)
                                                                           tctAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGa               >  1:1910755/1‑91 (MQ=255)
                                                                             tAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACa             >  2:806904/1‑91 (MQ=255)
                                                                             tAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACa             >  2:1161423/1‑91 (MQ=255)
                                                                              aaaGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACaa            >  1:1458088/1‑91 (MQ=255)
                                                                              aaaGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACaa            <  2:2491648/91‑1 (MQ=255)
                                                                               aaGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAg           <  1:2048213/91‑1 (MQ=255)
                                                                               aaGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAg           <  1:286998/91‑1 (MQ=255)
                                                                                aGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGt          <  1:1381510/91‑1 (MQ=255)
                                                                                        aGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTATTTGTAt  <  2:1661943/91‑1 (MQ=255)
                                                                                         |                                                                                         
TAAAGAAGAAAAGAAAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTATTTGTAT  >  USA300TCH1516_ALE/265540‑265718

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: