Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,448,484 A→G 9.7% intergenic (‑243/+34) yxeP_4 ← / ← hutI putative hydrolase YxeP/Imidazolonepropionase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,448,4840AG 9.7% 646.2 / 45.1 225intergenic (‑243/+34)yxeP_4/hutIputative hydrolase YxeP/Imidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base A (118/82);  new base G (11/11);  total (132/93)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 230 total aligned reads displayed.

AACTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACTTT  >  USA300TCH1516_ALE/2448394‑2448574
                                                                                          |                                                                                          
aaCTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTa                                                                                            >  2:133920/1‑91 (MQ=255)
     ttCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACACAAAGCTCTAAACATTCGTTTa                                                                                       >  2:2935504/1‑91 (MQ=255)
        tAACTTTTATAATTGAGACTCTATACAATAACGTGTTCTCGAGTATACTTGTTTTTTCAAACCACAACAAGCTCTACACGTTAGTTTtaat                                                                                    >  2:907297/1‑90 (MQ=255)
        tAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAAc                                                                                    <  2:1978632/91‑1 (MQ=255)
        tAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAAc                                                                                    >  2:3236044/1‑91 (MQ=255)
        tAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAAc                                                                                    >  1:1804360/1‑91 (MQ=255)
        tAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAAc                                                                                    >  2:2933348/1‑91 (MQ=255)
           cTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCaa                                                                                 <  2:857499/91‑1 (MQ=255)
               tataATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCt                                                                             >  1:1018053/1‑91 (MQ=255)
               tataATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCt                                                                             >  1:3742845/1‑91 (MQ=255)
               tataATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCt                                                                             <  1:638432/91‑1 (MQ=255)
               tataATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCt                                                                             <  2:554825/91‑1 (MQ=255)
                 taATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCCCAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTa                                                                           >  2:3112482/1‑91 (MQ=255)
                 taATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTa                                                                           <  2:1179715/91‑1 (MQ=37)
                 taATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTa                                                                           <  2:2644406/91‑1 (MQ=255)
                 taATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTa                                                                           >  1:2800822/1‑91 (MQ=255)
                 taATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTag                                                                          <  1:752815/91‑1 (MQ=25)
                  aaTTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTag                                                                          <  1:2014597/91‑1 (MQ=37)
                  aaTTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTag                                                                          <  1:3466787/91‑1 (MQ=37)
                  aaTTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTag                                                                          <  2:3275688/91‑1 (MQ=255)
                   aTTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTaga                                                                         >  2:1426301/1‑91 (MQ=255)
                    ttGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTagag                                                                        <  1:2698465/91‑1 (MQ=25)
                    ttGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTagag                                                                        >  2:1609453/1‑91 (MQ=255)
                     tGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGc                                                                       <  2:3154728/91‑1 (MQ=255)
                     tGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGc                                                                       >  1:1331425/1‑91 (MQ=255)
                        gaCTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACTACAAAAAGCTCTAAACATTAGTTTACACCAATGCTTAGAGCttt                                                                    >  2:900610/1‑91 (MQ=255)
                        gaCTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTAGAGCttt                                                                    <  2:2993198/91‑1 (MQ=25)
                        gaCTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCttt                                                                    <  1:744297/91‑1 (MQ=255)
                        aaCTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCttt                                                                    <  2:669887/90‑1 (MQ=255)
                         aCTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTAGAGCtttt                                                                   <  2:2659660/91‑2 (MQ=25)
                         aCTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTc                                                                   >  2:278537/1‑91 (MQ=255)
                         aCTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTc                                                                   <  1:2238843/91‑1 (MQ=255)
                          ctctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCt                                                                  >  1:1021814/1‑91 (MQ=255)
                           tctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACTACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTa                                                                 >  2:1773385/1‑91 (MQ=255)
                           tctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTa                                                                 <  2:2536149/91‑1 (MQ=255)
                           tctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTa                                                                 <  1:1386217/91‑1 (MQ=255)
                           tctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTa                                                                 <  1:425166/91‑1 (MQ=255)
                             tataCAAAAACGTGTTCTCGAATATACTTGTTTTTACAACCCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAat                                                               >  1:1360743/1‑91 (MQ=255)
                              ataCAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAatt                                                              >  2:2089075/1‑91 (MQ=255)
                                aCAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAattat                                                            <  2:1029000/91‑1 (MQ=255)
                                 cAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAGATCGCTCTAAACATTAGTTTAAACCAATGCTTAGGGCTTTCTAattatt                                                           <  2:3290270/91‑1 (MQ=255)
                                 cAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAattatt                                                           >  2:77514/1‑91 (MQ=255)
                                   aaaaCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTAtttt                                                         <  2:25326/91‑1 (MQ=255)
                                   aaaaCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTAtttt                                                         >  1:2583615/1‑91 (MQ=255)
                                    aaaCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTa                                                        <  2:3855738/91‑1 (MQ=255)
                                     aaCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTAt                                                       <  1:1733738/91‑1 (MQ=255)
                                     aaCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTAt                                                       <  2:3358596/91‑1 (MQ=255)
                                      aCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATg                                                      >  1:1353423/1‑91 (MQ=255)
                                      aCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATg                                                      <  1:2886083/91‑1 (MQ=255)
                                      aCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATg                                                      >  1:3843105/1‑91 (MQ=255)
                                       cGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGc                                                     >  2:485154/1‑91 (MQ=255)
                                         tgtTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCtt                                                   <  2:267536/91‑1 (MQ=255)
                                           ttCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTa                                                 >  1:2773260/1‑91 (MQ=255)
                                           ttCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTa                                                 >  2:2014201/1‑91 (MQ=255)
                                           ttCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTa                                                 >  2:2343036/1‑91 (MQ=255)
                                              tcGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTaaaa                                              <  2:792176/91‑1 (MQ=255)
                                                  atatACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATa                                          >  1:1855725/1‑91 (MQ=255)
                                                   tataCTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATAc                                         <  1:1457723/91‑1 (MQ=255)
                                                   tataCTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATAc                                         >  1:3006181/1‑91 (MQ=255)
                                                      aCTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACtgt                                      >  1:332727/1‑91 (MQ=255)
                                                        ttGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTgtgt                                    >  1:2078942/1‑91 (MQ=255)
                                                         tGTTTTTACAAACCACAAAAAGCTCTAAAGATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGtt                                   <  1:3167457/91‑1 (MQ=255)
                                                         tGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGtt                                   >  1:2554663/1‑91 (MQ=255)
                                                          gTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTa                                  >  1:91616/1‑91 (MQ=255)
                                                             tttACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCt                               <  1:2544805/91‑1 (MQ=255)
                                                             tttACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCt                               <  1:1685764/91‑1 (MQ=255)
                                                               tACAAACCACAAAAAGCTCTAAACATTCGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTAc                             >  1:2559928/1‑91 (MQ=255)
                                                                aCAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACg                            >  1:2823350/1‑91 (MQ=255)
                                                                aCAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACg                            >  1:1224084/1‑91 (MQ=255)
                                                                 cAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGa                           <  1:1429556/91‑1 (MQ=255)
                                                                   aaCCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATg                         >  1:3561550/1‑91 (MQ=255)
                                                                   aaCCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATg                         <  1:3059066/91‑1 (MQ=255)
                                                                   aaCCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATg                         <  2:1913025/91‑1 (MQ=255)
                                                                    aCCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGa                        >  2:3771079/1‑91 (MQ=255)
                                                                    aCCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGa                        >  2:2618347/1‑91 (MQ=255)
                                                                     ccACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGAc                       >  1:3247180/1‑91 (MQ=255)
                                                                          aaaaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTAc                  >  2:1465149/1‑91 (MQ=255)
                                                                           aaaaCTCTAAACATTAGTTTAAACTAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTAcc                 <  2:3664198/91‑1 (MQ=255)
                                                                           aaaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTAcc                 >  1:2701884/1‑91 (MQ=255)
                                                                           aaaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTAcc                 <  2:3146831/91‑1 (MQ=255)
                                                                           aaaGCTCTAAACATTAATTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTAcc                 <  2:3036594/91‑1 (MQ=255)
                                                                            aaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCg                >  2:2083070/1‑91 (MQ=255)
                                                                             aGCTCTAAACATTAGTTTAAACCAATGCTTCGAGCTTTCTAATTATTCTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGt               >  2:741840/1‑91 (MQ=255)
                                                                             aGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGt               <  2:1544613/91‑1 (MQ=255)
                                                                             aGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGt               >  2:2728098/1‑91 (MQ=255)
                                                                              gCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTc              >  2:2908245/1‑91 (MQ=255)
                                                                                 ctAAGCATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCttt           >  1:3201118/1‑91 (MQ=37)
                                                                                 ctAAGCATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCttt           >  2:3944989/1‑91 (MQ=37)
                                                                                    agcATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTCTTTTATGCTTTCAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtact        >  2:347976/3‑89 (MQ=25)
                                                                                    agcATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaat        >  1:877611/3‑91 (MQ=38)
                                                                                    agcATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaat        >  2:2849826/3‑91 (MQ=38)
                                                                                    agcATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaat        >  2:341766/3‑91 (MQ=37)
                                                                                    agcATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaat        >  2:3318389/3‑91 (MQ=38)
                                                                                     aCATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaata       <  1:819878/91‑1 (MQ=255)
                                                                                     aCATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaata       >  1:2679832/1‑91 (MQ=255)
                                                                                      cATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaataa      >  1:1505791/1‑91 (MQ=38)
                                                                                       aTTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAAc     >  2:1253802/1‑91 (MQ=255)
                                                                                       aTTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAAc     <  1:3861366/91‑1 (MQ=255)
                                                                                          aGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACttt  >  2:3495115/1‑91 (MQ=255)
                                                                                          aGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACttt  >  1:827011/1‑91 (MQ=255)
                                                                                          |                                                                                          
AACTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACTTT  >  USA300TCH1516_ALE/2448394‑2448574

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