Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 265,636 C→A 7.6% intergenic (+36/+128) 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,6360CA 7.6% 686.5 / 31.4 198intergenic (+36/+128)ugpQ_2/scn_1Glycerophosphodiester phosphodiesterase, cytoplasmic/Staphylococcal complement inhibitor
Reads supporting (aligned to +/- strand):  ref base C (83/100);  new base A (11/4);  total (94/104)
Fisher's exact test for biased strand distribution p-value = 5.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 207 total aligned reads displayed.

AGAAAAGAAAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTATTTGTATGGCAATTT  >  USA300TCH1516_ALE/265546‑265726
                                                                                          |                                                                                          
 gaaaagaaaGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATAcc                                                                                           >  2:942749/1‑91 (MQ=255)
 gaaaagaaaGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATAcc                                                                                           >  1:1385700/1‑91 (MQ=255)
  aaaagaaaGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCt                                                                                          >  2:3425475/1‑91 (MQ=255)
   aaagaaagACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCtt                                                                                         >  2:1709454/1‑91 (MQ=255)
   aaagaaagACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCtt                                                                                         >  2:546589/1‑91 (MQ=255)
       aaagaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCaga                                                                                     >  1:432866/1‑91 (MQ=255)
       aaagaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCaga                                                                                     >  2:1454732/1‑91 (MQ=255)
       aaagaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCaga                                                                                     <  2:2086911/91‑1 (MQ=255)
        aagaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCagag                                                                                    <  1:1995411/91‑1 (MQ=255)
        aagaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCagag                                                                                    >  1:1496177/1‑91 (MQ=255)
        aagaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCagag                                                                                    >  1:798558/1‑91 (MQ=255)
         agaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAgg                                                                                   <  1:2333650/91‑1 (MQ=255)
         agaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAgg                                                                                   <  1:3703566/91‑1 (MQ=255)
         agaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAgg                                                                                   <  1:669085/91‑1 (MQ=255)
         agaCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAgg                                                                                   <  1:1571186/91‑1 (MQ=255)
          gaCGAAACATACGTCGATCGTTCTATAAGATTTTTGTTTGAATAATATAAACAAAGACCTCTGAAGGTATCATCTTGATAACTTTAGAGGt                                                                                  <  2:1766055/91‑1 (MQ=18)
           aCGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTc                                                                                 >  2:1938856/1‑91 (MQ=255)
            cGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCt                                                                                >  2:1909156/1‑91 (MQ=255)
            cGAAACATACTTCGATCGTTCTATAAGAATGTTGTGTGAATAATATAAACAAAGACCTCTAAAGGTATCATCTTGATAACTTTAGAGGTCt                                                                                <  2:2828277/91‑1 (MQ=18)
             gAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCtt                                                                               <  1:3558522/91‑1 (MQ=255)
               aaCATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCtttt                                                                             >  2:39476/1‑91 (MQ=255)
                 cATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGATCTCAAGATGATACCTTCAGAGGTCTTTTTa                                                                           >  1:3104184/1‑91 (MQ=255)
                  aTACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTaa                                                                          <  1:65290/91‑1 (MQ=255)
                    aCTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATg                                                                        >  1:1382630/1‑91 (MQ=255)
                    aCTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATg                                                                        >  2:931112/1‑91 (MQ=255)
                        cGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGttac                                                                    <  2:453932/91‑3 (MQ=255)
                        cGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGc                                                                    <  2:337256/91‑1 (MQ=255)
                        cGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGc                                                                    >  1:1608086/1‑91 (MQ=255)
                         gATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGcc                                                                   >  1:197323/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:3729964/91‑1 (MQ=255)
                             gTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCt                                                               >  2:1837544/1‑91 (MQ=255)
                              ttCTATAAGAATTTTGTTTGAATAATATAAACAAAGGCCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTa                                                              <  2:3058870/91‑1 (MQ=255)
                              ttCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTa                                                              <  1:1405227/91‑1 (MQ=255)
                              ttCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTa                                                              <  1:3358872/91‑1 (MQ=255)
                               tCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTAt                                                             >  1:1449476/1‑91 (MQ=255)
                                 tataAGAATTTTGTTTGAATAATATAAACAAAGACCTGTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATgg                                                           >  1:3262659/1‑91 (MQ=255)
                                 tataAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATgg                                                           >  2:3686307/1‑91 (MQ=255)
                                  ataAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATggg                                                          <  1:86469/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                                                     >  1:714187/1‑91 (MQ=255)
                                        aTTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGc                                                    <  1:3310268/91‑1 (MQ=255)
                                         ttttGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCa                                                   >  2:734445/1‑91 (MQ=255)
                                          tttgtttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCaa                                                  <  2:2808322/91‑1 (MQ=255)
                                           ttgtttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAAt                                                 <  1:916474/91‑1 (MQ=255)
                                           ttgtttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAAt                                                 >  1:630831/1‑91 (MQ=255)
                                             gtttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCg                                               >  2:1227119/1‑91 (MQ=255)
                                             gtttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCg                                               >  1:2374599/1‑91 (MQ=255)
                                             gtttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCg                                               <  1:2161866/91‑1 (MQ=255)
                                              tttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGt                                              >  2:2006688/1‑91 (MQ=255)
                                              tttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGt                                              <  2:1630632/91‑1 (MQ=255)
                                              tttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGt                                              >  2:1778335/1‑91 (MQ=255)
                                               ttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGtt                                             >  1:3537528/1‑91 (MQ=255)
                                               ttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGtt                                             <  2:798063/91‑1 (MQ=255)
                                               ttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGtt                                             >  1:1213376/1‑91 (MQ=255)
                                               ttgAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGACTTTTTAATGTTGCCATCTATGGGATAGGCAATCGtt                                             <  1:3259736/91‑1 (MQ=255)
                                               ttgAATAATATAAACAAAGACCTCTAAAGTGATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGtt                                             <  2:2175437/91‑1 (MQ=255)
                                                    taataTAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCAtt                                        <  1:2816926/91‑1 (MQ=255)
                                                    taataTAAACAAAGACCTCTAAAGTTATCAAGATGACACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCAtt                                        <  1:2085654/91‑1 (MQ=255)
                                                     aataTAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCAtac                                       <  2:2447456/91‑3 (MQ=255)
                                                       tataAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGt                                     >  1:1087140/1‑91 (MQ=255)
                                                        ataAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGtt                                    <  1:266602/91‑1 (MQ=255)
                                                          aaaCAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTta                                  <  2:1177834/91‑1 (MQ=255)
                                                           aaCAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTtat                                 >  2:3831826/1‑91 (MQ=255)
                                                            aCAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTtata                                >  1:2240729/1‑91 (MQ=255)
                                                             aaaaGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTatat                               <  1:3412226/90‑1 (MQ=255)
                                                              aaaGACCTCTGAAGGTATCATCTTGATAACTTTAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATAtt                              >  2:914308/1‑91 (MQ=17)
                                                               aaGACCTCTGAAGGTATCATCTTGATAACTTTAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTc                             >  1:2597628/1‑91 (MQ=18)
                                                                 gACCTCTGAAGGTATCATCTTGATAACTTTAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCat                           >  1:1701812/1‑91 (MQ=21)
                                                                  aCCTCTGAAGGTATCATCTTGATAACTTTAGAGGTCTTTTTAATGTTGCCATATATGGGATAGGCAATCGTTTCATTAGTTTATATTAata                          >  1:3832413/1‑91 (MQ=18)
                                                                  aCCTCTGAAGGTATCATCTTGATAACTTTAGAGGACTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCata                          >  2:358108/1‑91 (MQ=21)
                                                                  aCCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCata                          >  2:2147749/1‑91 (MQ=255)
                                                                   ccTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCatat                         <  2:1824615/91‑1 (MQ=255)
                                                                   ccTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCatat                         <  1:182556/91‑1 (MQ=255)
                                                                   ccTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCatat                         <  1:182047/91‑1 (MQ=255)
                                                                     tctAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGa                       >  1:2722426/1‑91 (MQ=255)
                                                                     tctAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGa                       >  1:1910755/1‑91 (MQ=255)
                                                                     tctAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGa                       >  2:3462280/1‑91 (MQ=255)
                                                                       tAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACa                     >  2:1161423/1‑91 (MQ=255)
                                                                         aaGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAg                   <  1:286998/91‑1 (MQ=255)
                                                                          aGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTACGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGt                  <  2:2999249/91‑1 (MQ=255)
                                                                            ttATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTAt                <  1:2641975/91‑1 (MQ=255)
                                                                              aTCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTAttt              <  2:267657/91‑1 (MQ=255)
                                                                                  aGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTATTTGTAt          <  1:1333184/91‑1 (MQ=255)
                                                                                    ttgatAACTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTATTTGTATgg        >  2:2733802/2‑91 (MQ=255)
                                                                                     tgatACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTATTTGTATGGc       >  2:1084009/1‑91 (MQ=255)
                                                                                      gatACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTATTTGTATGGCa      >  2:439246/1‑91 (MQ=255)
                                                                                       atACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTATTTGTATGGCaa     <  2:3054229/91‑1 (MQ=255)
                                                                                       atACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTATTTGTATGGCaa     >  2:1383168/1‑91 (MQ=255)
                                                                                        tACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTATTTGTATGGCAAt    >  2:3466281/1‑91 (MQ=255)
                                                                                        tACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTATTTGTATGGCAAt    <  2:1382630/91‑1 (MQ=255)
                                                                                         aCCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTATTTGTATGGCAAtt   >  2:2600748/1‑91 (MQ=255)
                                                                                         aCCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTATTTGTATGGCAAtt   >  2:2603595/1‑91 (MQ=255)
                                                                                          ccTTCAGAGGTCTTTTTAATGCTGCCATCTATGGGATAGGCAATCGTTTCATACGTTTATATTCATATGACAAGAATTTGTATGGCAAttt  >  2:501242/1‑91 (MQ=255)
                                                                                          |                                                                                          
AGAAAAGAAAGACGAAACATACTTCGATCGTTCTATAAGAATTTTGTTTGAATAATATAAACAAAGACCTCTAAAGTTATCAAGATGATACCTTCAGAGGTCTTTTTAATGTTGCCATCTATGGGATAGGCAATCGTTTCATTCGTTTATATTCATATGACAAGTATTTGTATGGCAATTT  >  USA300TCH1516_ALE/265546‑265726

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: