Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,078,152 A→T 11.3% Q510L (CAG→CTG)  purL → Phosphoribosylformylglycinamidine synthase subunit PurL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,078,1520AT11.3% 306.4 / 26.7 97Q510L (CAG→CTG) purLPhosphoribosylformylglycinamidine synthase subunit PurL
Reads supporting (aligned to +/- strand):  ref base A (42/44);  new base T (6/5);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

ATCGAAATTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATAAC  >  USA300TCH1516_ALE/1078062‑1078237
                                                                                          |                                                                                     
aTCGAAATTTAAGTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACa                                                                                       >  1:3196124/1‑91 (MQ=255)
  cGAAATTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGt                                                                                     >  1:1359610/1‑91 (MQ=255)
     aaTTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGa                                                                                  <  1:951025/91‑1 (MQ=255)
     aaTTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGa                                                                                  <  1:1034782/91‑1 (MQ=255)
     aaTTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGa                                                                                  >  1:529608/1‑91 (MQ=255)
     aaTTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGa                                                                                  <  1:1089309/91‑1 (MQ=255)
      aTTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGAt                                                                                 <  1:687600/91‑1 (MQ=255)
      aTTTAATAGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGAt                                                                                 <  1:2069813/91‑1 (MQ=255)
         tAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGATAGa                                                                              <  1:3352845/91‑1 (MQ=255)
          aaTTGCCGTGGGTGCAACACCATTAGCAATGATAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGAt                                                                             >  1:2183357/1‑91 (MQ=255)
          aaTTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGAt                                                                             <  1:1655483/91‑1 (MQ=255)
          aaTTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGAt                                                                             <  1:2663381/91‑1 (MQ=255)
          aaTTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGAt                                                                             <  1:1422222/91‑1 (MQ=255)
            ttGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTc                                                                           <  1:3650288/91‑1 (MQ=255)
            ttGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTc                                                                           <  1:3650383/91‑1 (MQ=255)
            ttGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTc                                                                           <  1:3650346/91‑1 (MQ=255)
             tGCCGTGGGTGCAACACCATTAGCATTGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCa                                                                          >  1:729860/1‑91 (MQ=255)
             tGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTNAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGATAGAtttc                                                                          <  1:3519076/91‑3 (MQ=255)
             tGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCa                                                                          >  1:729720/1‑91 (MQ=255)
                cGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACg                                                                       >  1:1301887/1‑91 (MQ=255)
                cGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACg                                                                       <  1:766168/91‑1 (MQ=255)
                cGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACg                                                                       >  1:2559643/1‑91 (MQ=255)
                 gTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGa                                                                      <  1:2883037/91‑1 (MQ=255)
                 gTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGa                                                                      >  1:3061599/1‑91 (MQ=255)
                 gTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGa                                                                      >  1:3061661/1‑91 (MQ=255)
                     gTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAgg                                                                  >  1:2718107/1‑91 (MQ=255)
                      tGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGt                                                                 <  1:357555/91‑1 (MQ=255)
                      tGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGt                                                                 <  1:2768007/91‑1 (MQ=255)
                      tGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGt                                                                 <  1:357680/91‑1 (MQ=255)
                          acacCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATgg                                                             >  1:1705226/1‑91 (MQ=255)
                          acacCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATgg                                                             >  1:3065055/1‑91 (MQ=255)
                          acacCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATgg                                                             <  1:147274/91‑1 (MQ=255)
                             ccATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAg                                                          >  1:1944718/1‑91 (MQ=255)
                              cATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGa                                                         <  1:3409451/91‑1 (MQ=255)
                              cATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGa                                                         <  1:3408971/91‑1 (MQ=255)
                              aaTTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGa                                                         >  1:2614349/2‑91 (MQ=255)
                               aTTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGaa                                                        >  1:3279711/1‑91 (MQ=255)
                                  aGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCa                                                     >  1:106453/1‑91 (MQ=255)
                                  aGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCa                                                     <  1:124174/91‑1 (MQ=255)
                                   gCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCAt                                                    <  1:4027992/91‑1 (MQ=255)
                                   gCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCAt                                                    >  1:1395933/1‑91 (MQ=255)
                                     aaTGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGc                                                  <  1:3839510/91‑1 (MQ=255)
                                       tGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGa                                                <  1:2961627/91‑1 (MQ=255)
                                       tGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGa                                                <  1:893822/91‑1 (MQ=255)
                                       tGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGa                                                <  1:2051682/91‑1 (MQ=255)
                                        gACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGAc                                               >  1:254843/1‑91 (MQ=255)
                                         aCAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACa                                              <  1:2369300/91‑1 (MQ=255)
                                         aCAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACa                                              <  1:3050179/91‑1 (MQ=255)
                                         aCAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACa                                              <  1:3050693/91‑1 (MQ=255)
                                           aGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACAtt                                            >  1:1124083/1‑91 (MQ=255)
                                           aGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACAtt                                            >  1:1124098/1‑91 (MQ=255)
                                              ttGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCtt                                         >  1:1326329/1‑91 (MQ=255)
                                              ttGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCtt                                         >  1:418078/1‑91 (MQ=255)
                                              ttGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCtt                                         <  1:3084200/91‑1 (MQ=255)
                                                gTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTaa                                       >  1:1104455/1‑91 (MQ=255)
                                                gTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTaa                                       >  1:2818862/1‑91 (MQ=255)
                                                 tttAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAg                                      <  1:3881928/91‑1 (MQ=255)
                                                  ttAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGa                                     >  1:2433662/1‑91 (MQ=255)
                                                  ttAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGa                                     <  1:649452/91‑1 (MQ=255)
                                                    aaaTTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGaca                                   >  1:608493/1‑91 (MQ=255)
                                                    aaaTTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGaca                                   >  1:1393013/1‑91 (MQ=255)
                                                    aaaTTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGaca                                   >  1:71772/1‑91 (MQ=255)
                                                     aaTTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGacac                                  >  1:1172156/1‑91 (MQ=255)
                                                     aaTTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGacac                                  >  1:1171859/1‑91 (MQ=255)
                                                     aaTTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGacac                                  >  1:1171427/1‑91 (MQ=255)
                                                          tGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCagta                             <  1:2984007/91‑1 (MQ=255)
                                                             ttCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGtt                          >  1:2989378/1‑91 (MQ=255)
                                                             ttCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGtt                          >  1:3420205/1‑91 (MQ=255)
                                                              tctcCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGttt                         >  1:3089105/1‑91 (MQ=255)
                                                              tctcCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGttt                         >  1:1606881/1‑91 (MQ=255)
                                                              tctcCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGttt                         >  1:3950937/1‑91 (MQ=255)
                                                                   tGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGt                    <  1:1009415/91‑1 (MQ=255)
                                                                   tGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGt                    <  1:1009272/91‑1 (MQ=255)
                                                                   tGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGt                    <  1:3472663/91‑1 (MQ=255)
                                                                   tGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGt                    <  1:331414/91‑1 (MQ=255)
                                                                   tGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGt                    <  1:1009597/91‑1 (MQ=255)
                                                                    gAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTa                   <  1:1222613/91‑1 (MQ=255)
                                                                     aaaaGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTaa                  <  1:1682604/91‑1 (MQ=255)
                                                                     aaaaGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTaa                  <  1:1682366/91‑1 (MQ=255)
                                                                      aaagaaagAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAAt                 >  1:2721532/1‑91 (MQ=255)
                                                                      aaagaaagAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAAt                 <  1:72571/91‑1 (MQ=255)
                                                                       aagaaagaAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATg                <  1:2728676/91‑1 (MQ=255)
                                                                       aagaaagaAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATg                <  1:2410295/91‑1 (MQ=255)
                                                                       aagaaagaAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATg                <  1:2043295/91‑1 (MQ=255)
                                                                          aaagaaaTCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTAt             <  1:2596330/91‑1 (MQ=255)
                                                                           aagaaaTCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATc            <  1:3406543/91‑1 (MQ=255)
                                                                            agaaaTCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCt           >  1:47542/1‑91 (MQ=255)
                                                                                aTCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTtat       >  1:1531383/1‑91 (MQ=255)
                                                                                aTCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTtat       >  1:1861754/1‑91 (MQ=255)
                                                                                 tCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTtata      <  1:97427/91‑1 (MQ=255)
                                                                                    aTCAACTGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATaa   >  1:3667080/1‑91 (MQ=255)
                                                                                    aTCAACTGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATaa   >  1:3569707/1‑91 (MQ=255)
                                                                                     tCAACTGTTGATAGATTTAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATaac  >  1:1777358/1‑91 (MQ=255)
                                                                                     tCAACTGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATaac  >  1:2682605/1‑91 (MQ=255)
                                                                                     tCAACTGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATaac  >  1:3915314/1‑91 (MQ=255)
                                                                                     tCAACTGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATaac  >  1:3172516/1‑91 (MQ=255)
                                                                                     tCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATaac  >  1:158360/1‑91 (MQ=255)
                                                                                     tCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATaac  >  1:35590/1‑91 (MQ=255)
                                                                                          |                                                                                     
ATCGAAATTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATAAC  >  USA300TCH1516_ALE/1078062‑1078237

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