Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,309,800 Δ1 bp 12.0% coding (69/1032 nt) yfiZ_2 ← putative siderophore transport system permease protein YfiZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,309,8000C.12.0% 445.1 / 53.3 133coding (69/1032 nt)yfiZ_2putative siderophore transport system permease protein YfiZ
Reads supporting (aligned to +/- strand):  ref base C (56/61);  new base . (9/7);  total (65/68)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 135 total aligned reads displayed.

GTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGT  >  USA300TCH1516_ALE/2309710‑2309889
                                                                                            |                                                                                         
gTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATAc                                                                                           >  1:588723/1‑91 (MQ=255)
gTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATAc                                                                                           >  1:3491809/1‑91 (MQ=255)
  gAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGt                                                                                         <  1:487449/91‑1 (MQ=255)
  gAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGt                                                                                         <  1:3896641/91‑1 (MQ=255)
  gAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGt                                                                                         <  1:2827013/91‑1 (MQ=255)
  gAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGt                                                                                         <  1:1445797/91‑1 (MQ=255)
    aTGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCGCTATATACGTGa                                                                                       <  1:3198828/91‑1 (MQ=255)
    aTGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGa                                                                                       >  1:2027502/1‑91 (MQ=255)
     tGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAg                                                                                      >  1:2149014/1‑91 (MQ=255)
     tGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAg                                                                                      >  1:2152469/1‑91 (MQ=255)
      gATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATAGTGAg                                                                                      <  1:1515931/91‑1 (MQ=38)
      gATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAgt                                                                                     <  1:1325602/91‑1 (MQ=255)
       atatCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATAGTGAgt                                                                                     <  1:1828032/91‑1 (MQ=38)
       atatCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAgtg                                                                                    <  1:190869/91‑1 (MQ=255)
       atatCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAgtg                                                                                    <  1:1296775/91‑1 (MQ=255)
         atCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGtgtt                                                                                  >  1:1634348/1‑91 (MQ=255)
          tCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGtgttg                                                                                 >  1:3352990/1‑91 (MQ=255)
          tCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGtgttg                                                                                 >  1:1807782/1‑91 (MQ=255)
          tCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGtgttg                                                                                 <  1:1631083/91‑1 (MQ=255)
          tCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGtgttg                                                                                 >  1:2846979/1‑91 (MQ=255)
           cGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGtgttgt                                                                                <  1:549229/91‑1 (MQ=255)
            gCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATAGTGAGTTtgtg                                                                               <  1:1693739/91‑1 (MQ=25)
            gCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTtgtg                                                                               <  1:1676460/91‑1 (MQ=255)
            gCTAAAATTAATCTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTtgtg                                                                               >  1:3199174/1‑91 (MQ=255)
             cTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGc                                                                              >  1:342936/1‑91 (MQ=255)
              tAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGACGCACAACACTCACGTATATAGTGAGTTTGTGctt                                                                            <  1:1932623/91‑3 (MQ=17)
               aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATAGTGAGTTTGTGctt                                                                            <  1:1450632/91‑3 (MQ=21)
                aaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGACGCACAACACTCACGTATATAGTGAGTTTGTGctttc                                                                          <  1:1052232/91‑5 (MQ=9)
                 aattaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtct                                                                          <  1:1963773/91‑1 (MQ=255)
                  attaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctt                                                                         <  1:3360773/91‑1 (MQ=255)
                    taattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttt                                                                       <  1:2897822/91‑1 (MQ=255)
                     aattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtcttttt                                                                      >  1:1333460/1‑91 (MQ=255)
                      attTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttc                                                                     <  1:661767/91‑1 (MQ=255)
                      attTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttc                                                                     >  1:1437042/1‑91 (MQ=255)
                       ttttCGAAGAACCAATCGCCATATTTAAATATATATAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttct                                                                    >  1:2652277/1‑91 (MQ=255)
                       ttttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttct                                                                    >  1:2653359/1‑91 (MQ=255)
                        tttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttctt                                                                   <  1:1604221/91‑1 (MQ=255)
                            gaagaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTCAGTGTTGTGCGTCTTTTTCTTTTAt                                                               >  1:3870321/1‑91 (MQ=255)
                             aagaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATg                                                              <  1:3104959/91‑1 (MQ=255)
                               gaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGAt                                                            <  1:1840373/91‑1 (MQ=255)
                                aaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATc                                                           >  1:899480/1‑91 (MQ=255)
                                 aCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCt                                                          >  1:267627/1‑91 (MQ=255)
                                  ccAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTa                                                         >  1:1346758/1‑91 (MQ=255)
                                  ccAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTa                                                         >  1:1347201/1‑91 (MQ=255)
                                  ccAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTa                                                         <  1:3274478/91‑1 (MQ=255)
                                   cAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTAt                                                        >  1:3947165/1‑91 (MQ=255)
                                     aTCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATgg                                                      >  1:3470902/1‑91 (MQ=255)
                                       cGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCa                                                    >  1:2158551/1‑91 (MQ=255)
                                         ccATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGa                                                  <  1:1792116/91‑1 (MQ=255)
                                          cATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGAt                                                 <  1:1889061/91‑1 (MQ=255)
                                            tatTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGCGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTg                                               >  1:253254/1‑91 (MQ=255)
                                              tttAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGc                                             <  1:1392145/91‑1 (MQ=255)
                                                 aaaTATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTAc                                          >  1:1110139/1‑91 (MQ=255)
                                                    tatataCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTa                                       >  1:784605/1‑91 (MQ=255)
                                                      tataCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAtt                                     <  1:1952725/91‑1 (MQ=255)
                                                       ataCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAttt                                    <  1:2768699/91‑1 (MQ=255)
                                                       ataCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAttt                                    <  1:3741249/91‑1 (MQ=255)
                                                       ataCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAttt                                    <  1:3739560/91‑1 (MQ=255)
                                                          cAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCtt                                 >  1:24540/1‑91 (MQ=255)
                                                              tAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTtat                             <  1:519779/91‑1 (MQ=255)
                                                               aaaaaGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTtata                            <  1:1025995/91‑1 (MQ=255)
                                                               aaaaaGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTtata                            <  1:3170257/91‑1 (MQ=255)
                                                                aaaaGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAg                           >  1:1261987/1‑91 (MQ=255)
                                                                  aagaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                         >  1:3692775/1‑91 (MQ=255)
                                                                  aagaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                         >  1:819221/1‑91 (MQ=255)
                                                                   agaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCa                        <  1:745243/91‑1 (MQ=255)
                                                                   agaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCa                        <  1:2753651/91‑1 (MQ=255)
                                                                   agaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCa                        <  1:745751/91‑1 (MQ=255)
                                                                     gacgCACAACACTCACGTATATAGTGAGTTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAtt                      >  1:3892070/4‑91 (MQ=21)
                                                                     aaatCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGACAtt                      <  1:3272931/91‑1 (MQ=255)
                                                                      aagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCGAttt                     >  1:3456101/1‑91 (MQ=255)
                                                                      aagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAttt                     >  1:2921087/1‑91 (MQ=255)
                                                                      aagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAttt                     >  1:3642670/1‑91 (MQ=255)
                                                                        gCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTAc                   >  1:3259370/1‑91 (MQ=255)
                                                                        gCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTAc                   >  1:4050544/1‑91 (MQ=255)
                                                                        gCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTAc                   >  1:4050708/1‑91 (MQ=255)
                                                                          acaAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACtt                 <  1:2152867/91‑1 (MQ=255)
                                                                           caAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTc                <  1:2132608/91‑1 (MQ=255)
                                                                           caAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTc                <  1:2132592/91‑1 (MQ=255)
                                                                             aaCTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCa              >  1:3982294/1‑91 (MQ=255)
                                                                             aaCTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAATTACTTCCa              >  1:3982914/1‑91 (MQ=255)
                                                                               aCTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCAc             >  1:2322500/1‑91 (MQ=255)
                                                                                cTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCAcc            >  1:2709898/1‑91 (MQ=255)
                                                                                 tCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCt           <  1:3455929/91‑1 (MQ=255)
                                                                                 tCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCt           <  1:3626975/91‑1 (MQ=255)
                                                                                 tCACGATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCt           <  1:3456106/91‑1 (MQ=255)
                                                                                  cACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCtt          <  1:1405369/91‑1 (MQ=255)
                                                                                   acgtatataGTGAGTTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAc        >  1:3162937/4‑91 (MQ=255)
                                                                                   acgtatataGTGAGTTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAc        >  1:3631710/4‑91 (MQ=255)
                                                                                   acgtatataGTGAGTTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAc        >  1:1801376/4‑91 (MQ=255)
                                                                                    cgtatataGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAc        >  1:218015/3‑91 (MQ=255)
                                                                                     gtatataGTGAGTTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTGATAGCCATTTACTTCCACCTTACtt      >  1:3869599/2‑91 (MQ=255)
                                                                                       atataCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACtt      >  1:657/1‑91 (MQ=255)
                                                                                         ataCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACtttt    >  1:2558803/1‑91 (MQ=255)
                                                                                         ataCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACtttt    >  1:3900610/1‑91 (MQ=255)
                                                                                          taCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTg   <  1:3054195/91‑1 (MQ=255)
                                                                                          taCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTg   <  1:1837645/91‑1 (MQ=255)
                                                                                           aCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGt  >  1:2820790/1‑91 (MQ=255)
                                                                                           aCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGt  <  1:355113/91‑1 (MQ=255)
                                                                                           aCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGt  <  1:871719/91‑1 (MQ=255)
                                                                                            |                                                                                         
GTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGT  >  USA300TCH1516_ALE/2309710‑2309889

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: