Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,448,479 A→G 8.2% 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 8.2% 739.2 / 41.6 243intergenic (‑238/+39)yxeP_4/hutIputative hydrolase YxeP/Imidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base A (118/105);  new base G (13/7);  total (131/112)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

Only 100 of 261 total aligned reads displayed.

TGTCAAACTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAAT  >  USA300TCH1516_ALE/2448389‑2448568
                                                                                          |                                                                                         
tgtCAAACTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTaaa                                                                                           <  1:249263/91‑1 (MQ=255)
tgtCAAACTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTaaa                                                                                           >  1:1498452/1‑91 (MQ=255)
tgtCAAACTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTaaa                                                                                           >  1:755221/1‑91 (MQ=255)
tgtCAAACTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTaaa                                                                                           >  1:2828313/1‑91 (MQ=255)
 gtCAAACTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAAc                                                                                          >  2:2177760/1‑91 (MQ=255)
    aaaCTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAActaa                                                                                       <  1:3548828/91‑4 (MQ=255)
    aaaCTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAAGAtt                                                                                       >  1:3494902/1‑91 (MQ=255)
      aCTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAg                                                                                     >  1:3144843/1‑91 (MQ=255)
      aCTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAg                                                                                     >  1:1578463/1‑91 (MQ=255)
       cTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGt                                                                                    <  2:3006530/91‑1 (MQ=255)
        tATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGtt                                                                                   >  2:2051633/1‑91 (MQ=255)
        tATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGtt                                                                                   >  2:1575102/1‑91 (MQ=255)
        tATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGtt                                                                                   <  1:3810930/91‑1 (MQ=255)
          ttCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTa                                                                                 >  1:3557136/1‑91 (MQ=255)
             tAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAAc                                                                              <  2:3250717/91‑1 (MQ=255)
             tAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAAc                                                                              >  1:1442293/1‑91 (MQ=255)
             tAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAAc                                                                              >  1:2696992/1‑91 (MQ=255)
             tAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAAc                                                                              >  1:2549800/1‑91 (MQ=255)
              accTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAGCATTGGTTTAAAct                                                                             <  2:1285540/89‑2 (MQ=255)
              aaCTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAAcc                                                                             <  1:2294578/91‑1 (MQ=255)
              aaCTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAAcc                                                                             <  1:2293933/91‑1 (MQ=255)
              aaCTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAAcc                                                                             <  1:2287147/91‑1 (MQ=255)
                    tataATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCt                                                                       <  2:1103365/91‑1 (MQ=255)
                    tataATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAAAAAAAAAAAGCTCTAAACATTAGTTTAAACCAATGCt                                                                       <  2:3652883/91‑1 (MQ=255)
                      taATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTa                                                                     >  1:1331749/1‑91 (MQ=255)
                      taATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTa                                                                     >  2:3454643/1‑91 (MQ=255)
                       aaTTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTag                                                                    <  2:2586056/91‑1 (MQ=255)
                       aaTTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTag                                                                    >  2:3338412/1‑91 (MQ=255)
                       aaTTGAGACTCTATACAAAAACCTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTag                                                                    <  2:2005441/91‑1 (MQ=25)
                        aTTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTaga                                                                   >  1:1824398/1‑91 (MQ=255)
                         ttGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTagag                                                                  <  1:796557/91‑1 (MQ=255)
                           gagaCTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTAGAGCt                                                                <  2:1014337/91‑1 (MQ=25)
                              aCTCTATACAAAAACGTGTTCTCGAATATTCGTGTTTTTACAAATCACAAAAATCTCGAAACATTAGTTTAAACCAATGCTTAGAGCTTTc                                                             <  1:3102954/91‑1 (MQ=255)
                                tctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTa                                                           >  1:1012434/1‑91 (MQ=255)
                                tctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTa                                                           >  2:2482061/1‑91 (MQ=255)
                                tctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTa                                                           <  2:568084/91‑1 (MQ=255)
                                tctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTa                                                           >  2:634012/1‑91 (MQ=255)
                                tctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTa                                                           <  1:1750339/91‑1 (MQ=255)
                                 ctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTaa                                                          >  2:2637985/1‑91 (MQ=255)
                                 ctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTaa                                                          <  2:3476398/91‑1 (MQ=255)
                                  tataCAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAat                                                         <  2:1983306/91‑1 (MQ=255)
                                   ataCAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAatt                                                        >  1:3824190/1‑91 (MQ=255)
                                      cAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAattatt                                                     >  1:68127/1‑91 (MQ=255)
                                        aaaaCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTAtttt                                                   <  2:2131234/91‑1 (MQ=255)
                                           aCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATg                                                <  2:859059/91‑1 (MQ=255)
                                              tgtTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCtt                                             >  2:1353631/1‑91 (MQ=255)
                                              tgtTCTCGAATATACTTGTTTTGACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCtt                                             <  2:3122472/91‑1 (MQ=255)
                                                ttCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTa                                           >  2:1744282/1‑91 (MQ=255)
                                                ttCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTa                                           <  1:1037294/91‑1 (MQ=255)
                                                  ctcGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCGAATTATTTTATGCTTTaaa                                         <  2:3397758/91‑1 (MQ=255)
                                                   tcGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTaaaa                                        <  1:1399717/91‑1 (MQ=255)
                                                   tcGAATATACGTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTaaaa                                        <  2:1194090/91‑1 (MQ=255)
                                                     gAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGa                                      <  1:368739/91‑1 (MQ=255)
                                                       atatACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATa                                    >  1:1686711/1‑91 (MQ=255)
                                                        tataCTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATAc                                   <  1:2801339/91‑1 (MQ=255)
                                                        tataCTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATAc                                   <  2:2903884/91‑1 (MQ=255)
                                                        tataCTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATAc                                   <  1:2447533/91‑1 (MQ=255)
                                                         ataCTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACt                                  >  1:3048384/1‑91 (MQ=255)
                                                         ataCTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACt                                  >  1:1207483/1‑91 (MQ=255)
                                                          taCTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACtg                                 <  2:2449841/91‑1 (MQ=255)
                                                            cTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACtgtg                               >  2:2154664/1‑91 (MQ=255)
                                                            cTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTATGCTTTAAAAGATACTgtgt                              >  1:100999/1‑91 (MQ=255)
                                                              tGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGtt                             >  2:2149019/1‑91 (MQ=255)
                                                              tGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGtt                             <  1:2647176/91‑1 (MQ=255)
                                                                tttttACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTAt                           >  2:2569660/1‑91 (MQ=255)
                                                                tttttACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTAt                           <  2:1909369/91‑1 (MQ=255)
                                                                tttttACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACCGTGTTAt                           >  2:1755868/1‑91 (MQ=255)
                                                                 ttttACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATc                          <  2:1686711/91‑1 (MQ=255)
                                                                     aCAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACg                      >  1:503393/1‑91 (MQ=255)
                                                                     aCAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACg                      >  2:3631046/1‑91 (MQ=255)
                                                                      cAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGa                     >  1:375224/1‑91 (MQ=255)
                                                                      cAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGa                     >  1:3189231/1‑91 (MQ=255)
                                                                      cAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGa                     <  2:2379946/91‑1 (MQ=255)
                                                                      cAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGa                     >  1:609944/1‑91 (MQ=255)
                                                                        aaCCACAAAAAGCTCGAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTCATGCTTTAAAAGATACTGTGTTATCTACGATg                   >  1:3316969/1‑91 (MQ=255)
                                                                         aCCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCGACGATGa                  >  1:1938569/1‑91 (MQ=255)
                                                                          ccACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGAc                 <  1:970116/91‑1 (MQ=255)
                                                                          ccACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGAc                 <  1:3151138/91‑1 (MQ=255)
                                                                            ttagaaaGCTCTAAGCATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCt               >  2:952639/5‑91 (MQ=21)
                                                                            ttagaaaGCTCTAAGCATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCt               >  1:3007742/5‑91 (MQ=21)
                                                                            acaAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCt               <  2:684995/91‑1 (MQ=255)
                                                                            acaAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCt               <  1:3267849/91‑1 (MQ=255)
                                                                            acaAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCt               <  1:1646399/91‑1 (MQ=255)
                                                                             tagaaaGCTCTAAGCATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCtt              >  1:2765460/4‑91 (MQ=21)
                                                                              aaaaaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTa             >  1:2737056/1‑91 (MQ=255)
                                                                              aaaaaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTa             >  2:528090/1‑91 (MQ=255)
                                                                               aaaaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTAc            <  1:3420390/91‑1 (MQ=255)
                                                                                aaaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTAcc           <  1:2917667/91‑1 (MQ=255)
                                                                                aaaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTAcc           <  2:3342550/91‑1 (MQ=255)
                                                                                aaaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTAcc           >  2:3236642/1‑91 (MQ=255)
                                                                                 aaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCg          <  2:1498452/91‑1 (MQ=255)
                                                                                    ctctAAGCATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCt       >  1:190015/1‑91 (MQ=37)
                                                                                    ctctAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCt       >  1:659045/1‑91 (MQ=255)
                                                                                    ctctAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCt       <  2:2810223/91‑1 (MQ=255)
                                                                                    ctctAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCt       >  2:871345/1‑91 (MQ=255)
                                                                                     tctAAGCATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCtt      >  1:699738/1‑91 (MQ=25)
                                                                                     tctAAGCATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCtt      >  1:213830/1‑91 (MQ=25)
                                                                                       tAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTta    >  2:734164/1‑91 (MQ=255)
                                                                                        aaaCATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaa   >  2:3631951/1‑91 (MQ=255)
                                                                                         aaCATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaat  <  2:3202497/91‑1 (MQ=255)
                                                                                          |                                                                                         
TGTCAAACTATTCTAACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAAT  >  USA300TCH1516_ALE/2448389‑2448568

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: