Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,448,479 A→G 6.4% intergenic (‑238/+39) 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,4790AG 6.4% 660.0 / 23.1 216intergenic (‑238/+39)yxeP_4/hutIputative hydrolase YxeP/Imidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base A (121/80);  new base G (3/11);  total (125/91)
Fisher's exact test for biased strand distribution p-value = 9.21e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

Only 100 of 226 total aligned reads displayed.

TGTCAAACTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATA  >  USA300TCH1516_ALE/2448389‑2448569
                                                                                          |                                                                                          
tgtCAAACTATTCTAACTTTTATGATTGAGACTCTATCCAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTaaa                                                                                            <  2:2517770/91‑1 (MQ=255)
tgtCAAACTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTaaa                                                                                            >  1:2906128/1‑91 (MQ=255)
tgtCAAACTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTaaa                                                                                            >  1:2179040/1‑91 (MQ=255)
  tCAAACTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACa                                                                                          <  2:104850/91‑1 (MQ=255)
   cAAACTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACAt                                                                                         <  2:3637738/91‑1 (MQ=255)
    aaaCTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCAAAACAtt                                                                                        >  1:1961694/1‑91 (MQ=255)
      aCTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACACAAAGCTCTAAACCTTAg                                                                                      >  2:287941/1‑91 (MQ=255)
      aCTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAg                                                                                      >  1:575207/1‑91 (MQ=255)
       cTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGt                                                                                     >  2:2665829/1‑91 (MQ=255)
          ttCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACACAAAGCTCTAAACATTCGTTTa                                                                                  >  2:2935504/1‑91 (MQ=255)
           tCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTaa                                                                                 >  1:856954/1‑91 (MQ=255)
             tAACTTTTATAATTGAGACTCTATACAATAACGTGTTCTCGAGTATACTTGTTTTTTCAAACCACAACAAGCTCTACACGTTAGTTTtaat                                                                               >  2:907297/1‑90 (MQ=255)
             tAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACATACCACAAAAAGCTCTAAACATTAGTTTAAAc                                                                               >  2:131553/1‑91 (MQ=255)
             tAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAAc                                                                               <  2:1978632/91‑1 (MQ=255)
              aaCTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACACAAAGCTCTAAACATTAGTTTAAAcc                                                                              >  2:2627494/1‑91 (MQ=255)
              aaCTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAAcc                                                                              >  2:3869887/1‑91 (MQ=255)
              aaCTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAAcc                                                                              >  2:3488915/1‑91 (MQ=255)
              aaCTTTTATAATTGAGACTCTATAAAAAAACGTAGTCTCGAATATACTTGTTTTTACAAACCACAAAAAGATATAAACATTAGTTTAAAca                                                                              >  2:1475800/1‑90 (MQ=255)
               aCTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCa                                                                             <  1:3231882/91‑1 (MQ=255)
                cTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCaa                                                                            <  2:857499/91‑1 (MQ=255)
                    tataATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCt                                                                        <  1:638432/91‑1 (MQ=255)
                    tataATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCt                                                                        >  1:3742700/1‑91 (MQ=255)
                    tataATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCt                                                                        >  1:1018053/1‑91 (MQ=255)
                    tataATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCt                                                                        <  2:554825/91‑1 (MQ=255)
                    tataATTGAGACTCTAAACAAAAACGTGTTCTCGAATATACTTGTTTTAAAAAAAAAAAAAAAGCTCTAAACATTAGTTTAAACCAATGCt                                                                        <  2:1386713/91‑1 (MQ=255)
                     ataATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCtt                                                                       <  2:6725/91‑1 (MQ=255)
                      taATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCCCAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTa                                                                      >  2:3112482/1‑91 (MQ=255)
                      taATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTa                                                                      <  2:1179715/91‑1 (MQ=37)
                      taATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTa                                                                      >  1:2800822/1‑91 (MQ=255)
                       aaTTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTag                                                                     <  2:3276743/91‑1 (MQ=255)
                         ttGAGACTCTATACAAAAGCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTagag                                                                   <  2:2742661/91‑1 (MQ=255)
                         ttGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTagag                                                                   <  2:2536027/91‑1 (MQ=25)
                         ttGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTagag                                                                   <  1:2698465/91‑1 (MQ=25)
                         ttGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTagag                                                                   >  2:1609453/1‑91 (MQ=255)
                         ttGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTagag                                                                   >  2:2657426/1‑91 (MQ=255)
                          tGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGc                                                                  >  1:1331425/1‑91 (MQ=255)
                          tGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGc                                                                  >  1:2703586/1‑91 (MQ=255)
                             gaCTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACTACAAAAAGCTCTAAACATTAGTTTACACCAATGCTTAGAGCttt                                                               >  2:900610/1‑91 (MQ=255)
                             gaCTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCttt                                                               >  2:1560298/1‑91 (MQ=255)
                               ctctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCt                                                             >  1:1021814/1‑91 (MQ=255)
                                tctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACTACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTa                                                            >  2:1773385/1‑91 (MQ=255)
                                tctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTa                                                            <  1:1386217/91‑1 (MQ=255)
                                  tataCAAAAACGTGTTCTCGAATATACTTGTTTTTACAACCCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAat                                                          >  1:1360743/1‑91 (MQ=255)
                                     aCAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAattat                                                       >  2:2586784/1‑91 (MQ=255)
                                      cAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAattatt                                                      >  1:1907576/1‑91 (MQ=255)
                                      cAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAattatt                                                      >  2:77514/1‑91 (MQ=255)
                                        aaaaCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTAtttt                                                    <  1:1650975/91‑1 (MQ=255)
                                        aaaaCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTAtttt                                                    >  1:2583615/1‑91 (MQ=255)
                                         aaaCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTa                                                   <  2:3855738/91‑1 (MQ=255)
                                          aaCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTAt                                                  <  2:3358596/91‑1 (MQ=255)
                                          aaCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTAt                                                  <  1:1733738/91‑1 (MQ=255)
                                           aCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATg                                                 <  1:2886083/91‑1 (MQ=255)
                                           aCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATg                                                 >  1:1353423/1‑91 (MQ=255)
                                            cGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGc                                                >  2:485154/1‑91 (MQ=255)
                                             gtgtTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCt                                               <  2:3147152/91‑1 (MQ=255)
                                              tgtTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCtt                                              <  2:267536/91‑1 (MQ=255)
                                              tgtTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCtt                                              <  1:3019867/91‑1 (MQ=255)
                                                ttCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTa                                            >  2:2014201/1‑91 (MQ=255)
                                                ttCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTa                                            >  1:2773260/1‑91 (MQ=255)
                                                ttCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTa                                            >  2:2343036/1‑91 (MQ=255)
                                                   tcGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTaaaa                                         <  2:792176/91‑1 (MQ=255)
                                                        tataCTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATAc                                    <  1:1457723/91‑1 (MQ=255)
                                                          taCTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACtg                                  <  2:2991907/91‑1 (MQ=255)
                                                          taCTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACtg                                  <  1:345685/91‑1 (MQ=255)
                                                           aCTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACtgt                                 >  1:332420/1‑91 (MQ=255)
                                                            cTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACtgtg                                >  1:2575302/1‑91 (MQ=255)
                                                             ttGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTgtgt                               >  1:2078910/1‑91 (MQ=255)
                                                             ttGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTgtgt                               >  1:879668/1‑91 (MQ=255)
                                                              tGTTTTTACAAACCACAAAAAGCTCTAAATATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGtt                              <  1:3167442/91‑1 (MQ=255)
                                                              tGTTTTTACAAACCACAAAAAGCTCTAAAGATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGtt                              <  1:3167479/91‑1 (MQ=255)
                                                              tGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGtt                              <  2:2087078/91‑1 (MQ=255)
                                                              tGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGtt                              >  1:2554663/1‑91 (MQ=255)
                                                              tGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGtt                              >  1:1418504/1‑91 (MQ=255)
                                                              tGTTTTTACAAACCACAAAAAGCGCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGtt                              <  2:3215981/91‑1 (MQ=255)
                                                               gTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTa                             <  2:2823447/91‑1 (MQ=255)
                                                                  tttACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCt                          <  1:1685764/91‑1 (MQ=255)
                                                                  tttACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCt                          <  1:2544805/91‑1 (MQ=255)
                                                                  tttACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCt                          <  1:973790/91‑1 (MQ=255)
                                                                  tttACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCt                          >  1:598807/1‑91 (MQ=255)
                                                                     aCAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACg                       >  1:1224084/1‑91 (MQ=255)
                                                                     aCAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACg                       <  1:1825333/91‑1 (MQ=255)
                                                                      cAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGa                      >  1:3768088/1‑91 (MQ=255)
                                                                      cAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGa                      >  1:2269533/1‑91 (MQ=255)
                                                                        aaCCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATg                    <  2:1913025/91‑1 (MQ=255)
                                                                        aaCCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATg                    >  1:3561550/1‑91 (MQ=255)
                                                                         aCCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGa                   >  2:3771079/1‑91 (MQ=255)
                                                                           cacaAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGAcc                 <  1:282554/91‑1 (MQ=255)
                                                                           cacaAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGAcc                 <  2:3264350/91‑1 (MQ=255)
                                                                             caAAAAGCTCTAAACATTAGTTTAAACCAATGGTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCtt               <  2:2928620/91‑1 (MQ=255)
                                                                              aaaaaGCTTTAACAACGAGCTTAAATCAATGCTAAGAGCTTTCAAACTATGCTATACATTAACAAATACTGTGTTATCTACGATGACCTTa              >  2:329220/1‑91 (MQ=255)
                                                                                aaaaCTCTAAACATTAGTTTAAACTAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTAcc            <  2:3664198/91‑1 (MQ=255)
                                                                                aaaGCTCTAAACATTATTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTAcc            <  2:3036330/91‑1 (MQ=255)
                                                                                aaaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTAcc            <  2:3146779/91‑1 (MQ=255)
                                                                                 aaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCg           >  2:2083070/1‑91 (MQ=255)
                                                                                  aGCTCTAAACATTAGTTTAAACCAATGCTTCGAGCTTTCTAATTATTCTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGt          >  2:741840/1‑91 (MQ=255)
                                                                                  aGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGt          >  1:779802/1‑91 (MQ=255)
                                                                                     tctAAGCATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCtt       >  1:745456/1‑91 (MQ=25)
                                                                                        aaaCATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaa    <  2:1453980/91‑1 (MQ=255)
                                                                                        aaaCATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaa    >  1:3860625/1‑91 (MQ=255)
                                                                                          aCATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaata  >  1:2679832/1‑91 (MQ=255)
                                                                                          |                                                                                          
TGTCAAACTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATA  >  USA300TCH1516_ALE/2448389‑2448569

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: