Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,216,257 G→A 38.5% S178L (TCA→TTA)  sceD ← putative transglycosylase SceD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,216,2570GA38.5% 169.3 / 216.3 197S178L (TCA→TTA) sceDputative transglycosylase SceD
Reads supporting (aligned to +/- strand):  ref base G (61/60);  new base A (34/42);  total (95/102)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 203 total aligned reads displayed.

GGAGCATTTGCTGGTGATACACCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACTGATGAACCTTCTG  >  USA300TCH1516_ALE/2216167‑2216347
                                                                                          |                                                                                          
ggAGCATTTGCTGGTGATACACCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTg                                                                                            >  1:574009/1‑91 (MQ=255)
      tttGCTGGTGATACACCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTg                                                                                      <  1:2114791/91‑1 (MQ=255)
      tttGCTGGTGATACACCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTg                                                                                      <  2:3662290/91‑1 (MQ=255)
       ttGCTGGTGATACACCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTgg                                                                                     <  2:1668603/91‑1 (MQ=255)
        tGCTGGTGATACACCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGa                                                                                    <  2:3055015/91‑1 (MQ=255)
        tGCTGGTGATACACCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGa                                                                                    <  2:1435331/91‑1 (MQ=255)
         gCTGGTGATACACCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGAt                                                                                   <  1:625332/91‑1 (MQ=255)
         gCTGGTGATACACCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGAt                                                                                   >  2:2078606/1‑91 (MQ=255)
          cTGGTGATACACCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGAtt                                                                                  >  2:2341116/1‑91 (MQ=255)
              tGATACACCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACa                                                                              >  2:2451422/1‑91 (MQ=255)
               gATACACCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAg                                                                             <  1:2744496/91‑1 (MQ=255)
               gATACACCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAg                                                                             >  2:2398448/1‑91 (MQ=255)
                aTACACCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGc                                                                            <  2:3457603/91‑1 (MQ=255)
                  acacCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCAt                                                                          >  1:1989198/1‑91 (MQ=255)
                   cacCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATg                                                                         <  1:1211291/91‑1 (MQ=255)
                     ccTTTATATTTAGCAGGTGCTACTGAATCCCCAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGaa                                                                       >  2:3249566/1‑91 (MQ=255)
                     ccTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGaa                                                                       <  1:3336201/91‑1 (MQ=255)
                     ccTTTATAGTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGaa                                                                       <  2:1847202/91‑1 (MQ=255)
                      cTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAat                                                                      >  2:3604819/1‑91 (MQ=255)
                      cTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTAGATTTACAGCATGAat                                                                      <  1:1868347/91‑1 (MQ=255)
                         tataTTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATAtt                                                                   >  1:2626008/1‑91 (MQ=255)
                         tataTTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATAtt                                                                   >  2:2726913/1‑91 (MQ=255)
                          atatTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTg                                                                  >  2:3836595/1‑91 (MQ=255)
                           tgtttAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGc                                                                 <  2:1698511/89‑1 (MQ=255)
                           tatTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGc                                                                 >  1:2533171/1‑91 (MQ=255)
                           tatTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGc                                                                 <  1:1365466/91‑1 (MQ=255)
                           tatTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGc                                                                 <  1:3171999/91‑1 (MQ=255)
                            atTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGcc                                                                >  2:2867677/1‑91 (MQ=255)
                               tAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCAcc                                                             >  1:3915523/1‑91 (MQ=255)
                               tAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCAcc                                                             <  1:2957016/91‑1 (MQ=255)
                                 gCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTg                                                           <  1:512940/91‑1 (MQ=255)
                                 gCAGCTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTg                                                           <  1:3288346/91‑1 (MQ=255)
                                   aGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGAt                                                         >  1:998408/1‑91 (MQ=255)
                                   aGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGAt                                                         <  2:2093510/91‑1 (MQ=255)
                                   aGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGACACCTGAt                                                         <  2:2854421/91‑1 (MQ=255)
                                   aGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGAt                                                         <  1:2332937/91‑1 (MQ=255)
                                   aGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGAt                                                         <  1:158915/91‑1 (MQ=255)
                                    ggTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGAtt                                                        <  1:2695/91‑1 (MQ=255)
                                     gTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCCCCTGATTc                                                       >  2:3390333/1‑91 (MQ=255)
                                     gTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTc                                                       >  1:104210/1‑91 (MQ=255)
                                     gTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTc                                                       <  1:1275348/91‑1 (MQ=255)
                                     gTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTc                                                       >  1:2983573/1‑91 (MQ=255)
                                      tGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCa                                                      >  1:2513592/1‑91 (MQ=255)
                                      tGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCa                                                      >  1:2512325/1‑91 (MQ=255)
                                       gCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCAc                                                     >  2:3006463/1‑91 (MQ=255)
                                         tACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGt                                                   >  1:1077413/1‑91 (MQ=255)
                                          aCTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGtt                                                  >  1:692436/1‑91 (MQ=255)
                                          aCTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGtt                                                  >  2:991958/1‑91 (MQ=255)
                                           cTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTg                                                 <  1:2030346/91‑1 (MQ=255)
                                            tGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGa                                                <  2:1972249/91‑1 (MQ=255)
                                              aaTCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGc                                              <  2:228607/91‑1 (MQ=255)
                                              aaTCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGc                                              <  1:914775/91‑1 (MQ=255)
                                              aaTCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGc                                              <  1:2820919/91‑1 (MQ=255)
                                               aTCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCa                                             <  1:775834/91‑1 (MQ=255)
                                                tCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCaa                                            <  1:3867369/91‑1 (MQ=255)
                                                tCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCaa                                            <  2:827909/91‑1 (MQ=255)
                                                 cccAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAAt                                           <  2:1268151/91‑1 (MQ=255)
                                                 cccAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAAt                                           <  2:28656/91‑1 (MQ=255)
                                                   cAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAAttt                                         >  2:3835933/1‑91 (MQ=255)
                                                    aaGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATCTg                                        >  1:3593776/1‑91 (MQ=255)
                                                     aGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGt                                       >  1:898618/1‑91 (MQ=255)
                                                       ttGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGttt                                     <  2:1740200/91‑1 (MQ=255)
                                                       ttGATTGTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGttt                                     <  2:185579/91‑1 (MQ=255)
                                                         gATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTa                                   <  2:751565/91‑1 (MQ=255)
                                                           ttGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGa                                 <  1:1513954/91‑1 (MQ=255)
                                                           ttGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGa                                 <  1:1383349/91‑1 (MQ=255)
                                                           ttGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGa                                 >  1:428819/1‑91 (MQ=255)
                                                            tGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGAt                                >  1:2190268/1‑91 (MQ=255)
                                                             gTAAGAATTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCACTTAGCTTTAGATg                               >  2:3087908/1‑91 (MQ=255)
                                                              tAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGa                              <  2:1496953/91‑1 (MQ=255)
                                                                aGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGc                            >  1:1885097/1‑91 (MQ=255)
                                                                  aaTTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCAt                          >  2:3119374/1‑91 (MQ=255)
                                                                  aaTTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCAt                          <  1:1728721/91‑1 (MQ=255)
                                                                   aTTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCAtt                         <  1:735753/91‑1 (MQ=255)
                                                                   aTTGATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCAtt                         <  1:3717620/91‑1 (MQ=255)
                                                                    ttGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTc                        <  2:3544073/91‑1 (MQ=255)
                                                                      gATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGCTGAGCATTcac                      >  1:3584932/1‑91 (MQ=255)
                                                                      gATACTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTcac                      >  1:1012608/1‑91 (MQ=255)
                                                                         aCTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACAtt                   >  1:1117009/1‑91 (MQ=255)
                                                                          cTTACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACAttt                  <  2:2575314/91‑1 (MQ=255)
                                                                           ttACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTa                 >  1:1561455/1‑91 (MQ=255)
                                                                           ttACCAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTa                 >  1:2270165/1‑91 (MQ=255)
                                                                              ccAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACtg              >  1:872040/1‑91 (MQ=255)
                                                                              ccAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACtg              >  2:1084750/1‑91 (MQ=255)
                                                                              ccAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACtg              >  1:1750480/1‑91 (MQ=255)
                                                                              ccAGCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACtg              >  1:1670493/1‑91 (MQ=255)
                                                                               cAGCTGCACCTGATGTTGGATTTACAGCATGATTATTGCCACCTGATTCACGTTGGGCAATTTGTTTTAGATGAGCATTCACATTTACtga             <  2:1028428/91‑1 (MQ=255)
                                                                                aGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACtgat            <  1:2897135/91‑1 (MQ=255)
                                                                                 gCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACtgatg           >  2:976626/1‑91 (MQ=255)
                                                                                 gCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACtgatg           >  2:1803768/1‑91 (MQ=255)
                                                                                 gCTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACtgatg           <  1:927569/91‑1 (MQ=255)
                                                                                  cTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACtgatga          >  1:1887672/1‑91 (MQ=255)
                                                                                  cTGCACCTAATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACtgatga          >  1:3434266/1‑91 (MQ=255)
                                                                                    gCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACTGATGAAc        <  2:2965199/91‑1 (MQ=255)
                                                                                    gCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACTGATGAAc        >  1:1372422/1‑91 (MQ=255)
                                                                                    gCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACAATTACTGATGAAc        <  1:3052542/91‑1 (MQ=255)
                                                                                      aCCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACTGATGAACCt      >  2:2739336/1‑91 (MQ=255)
                                                                                      aCCTGACGTTGGATTTACAGCATGACTATTGCCACCTGCTTCACGTTGTGCAATTTGTTTTAGATGAGCATTCACATTTACTGATGATCCt      >  2:2652513/1‑91 (MQ=255)
                                                                                          gATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACTGATGAACCTTCtg  <  1:3083311/91‑1 (MQ=255)
                                                                                          gATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACTGATGAACCTTCtg  <  2:1211199/91‑1 (MQ=255)
                                                                                          |                                                                                          
GGAGCATTTGCTGGTGATACACCTTTATATTTAGCAGGTGCTACTGAATCCCAAGTTGATTGTAAGAATTGATACTTACCAGCTGCACCTGATGTTGGATTTACAGCATGAATATTGCCACCTGATTCACGTTGAGCAATTTGTTTTAGATGAGCATTCACATTTACTGATGAACCTTCTG  >  USA300TCH1516_ALE/2216167‑2216347

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: