Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,309,807 Δ1 bp 10.0% coding (62/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,8070G.10.0% 368.8 / 34.6 110coding (62/1032 nt)yfiZ_2putative siderophore transport system permease protein YfiZ
Reads supporting (aligned to +/- strand):  ref base G (49/50);  new base . (7/4);  total (56/54)
Fisher's exact test for biased strand distribution p-value = 5.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 127 total aligned reads displayed.

ATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTT  >  USA300TCH1516_ALE/2309717‑2309897
                                                                                            |                                                                                          
atatCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAgtg                                                                                            <  1:1296775/91‑1 (MQ=255)
atatCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAgtg                                                                                            <  1:190869/91‑1 (MQ=255)
  atCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGtgtt                                                                                          >  1:1634348/1‑91 (MQ=255)
   tCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGtgttg                                                                                         >  1:3660081/1‑91 (MQ=255)
   tCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGtgttg                                                                                         >  1:2846979/1‑91 (MQ=255)
   tCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGtgttg                                                                                         >  1:1807953/1‑91 (MQ=255)
   tCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGtgttg                                                                                         <  1:1631083/91‑1 (MQ=255)
    cGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGtgttgt                                                                                        <  1:549229/91‑1 (MQ=255)
     gCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATAGTGAGTTtgtg                                                                                       <  1:1693739/91‑1 (MQ=25)
     gCTAAAATTAATCTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTtgtg                                                                                       >  1:3199174/1‑91 (MQ=255)
     gCTAAAATTAATCTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTtgtg                                                                                       >  1:3198640/1‑91 (MQ=255)
      cTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGc                                                                                      >  1:342936/1‑91 (MQ=255)
       tAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGACGCACAACACTCACGTATATAGTGAGTTTGTGctt                                                                                    <  1:2415038/91‑3 (MQ=17)
       tAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGACGCACAACACTCACGTATATAGTGAGTTTGTGctt                                                                                    <  1:1932623/91‑3 (MQ=17)
        aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATAGTGAGTTTGTGctt                                                                                    <  1:1450632/91‑3 (MQ=21)
         aaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGACGCACAACACTCACGTATATAGTGAGTTTGTGctttc                                                                                  <  1:1052232/91‑5 (MQ=9)
          aattaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtct                                                                                  <  1:622983/91‑1 (MQ=255)
          aattaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtct                                                                                  <  1:1963773/91‑1 (MQ=255)
           attaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctt                                                                                 <  1:3360773/91‑1 (MQ=255)
            ttaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtcttt                                                                                <  1:113510/91‑1 (MQ=255)
              aattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtcttttt                                                                              >  1:1333460/1‑91 (MQ=255)
              aattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtcttttt                                                                              >  1:2859769/1‑91 (MQ=255)
                ttttCGAAGAACCAATCGCCATATTTAAATATATATAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttct                                                                            >  1:2652277/1‑91 (MQ=255)
                ttttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttct                                                                            >  1:1767664/1‑91 (MQ=255)
                ttttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttct                                                                            >  1:2653359/1‑91 (MQ=255)
                 tttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttctt                                                                           <  1:1604221/91‑1 (MQ=255)
                  ttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttcttt                                                                          >  1:2171759/1‑91 (MQ=255)
                   tCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttctttt                                                                         <  1:1359783/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:491508/1‑91 (MQ=255)
                           ccAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTa                                                                 <  1:3274478/91‑1 (MQ=255)
                           ccAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTa                                                                 >  1:1347201/1‑91 (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)
                                        ttAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCt                                                    <  1:3442777/91‑1 (MQ=255)
                                          aaaTATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTAc                                                  >  1:1110139/1‑91 (MQ=255)
                                          aaaTATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTAc                                                  >  1:4004141/1‑91 (MQ=255)
                                          aaaTATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGCGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTAc                                                  <  1:1445249/91‑1 (MQ=255)
                                           aaTATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACt                                                 >  1:3348928/1‑91 (MQ=255)
                                             tatataCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTa                                               >  1:784605/1‑91 (MQ=255)
                                                ataCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAttt                                            <  1:2768699/91‑1 (MQ=255)
                                                ataCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAttt                                            <  1:3741249/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)
                                                          aaagaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGc                                  >  1:2037426/1‑91 (MQ=255)
                                                           aagaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                                 >  1:3692775/1‑91 (MQ=255)
                                                            agaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCa                                <  1:745243/91‑1 (MQ=255)
                                                            agaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCa                                <  1:745256/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)
                                                               aagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAttt                             >  1:2921087/1‑91 (MQ=255)
                                                               aagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAttt                             >  1:3642670/1‑91 (MQ=255)
                                                                 gCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTAc                           >  1:4050708/1‑91 (MQ=255)
                                                                 gCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTAc                           >  1:4050544/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)
                                                                    caAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTc                        >  1:1948769/1‑91 (MQ=255)
                                                                      aaCTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCa                      >  1:3982294/1‑91 (MQ=255)
                                                                      aaCTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAATTACTTCCa                      >  1:3982914/1‑91 (MQ=255)
                                                                        aCTCACGTATATAGTGAGTTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCAcc                    >  1:2709528/1‑91 (MQ=38)
                                                                         cTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCAcc                    <  1:2575229/91‑1 (MQ=255)
                                                                         cTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCAcc                    >  1:2709898/1‑91 (MQ=255)
                                                                          tCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCt                   <  1:3455929/91‑1 (MQ=255)
                                                                          tCACGATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCt                   <  1:3457231/91‑1 (MQ=255)
                                                                          tCACGATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCt                   <  1:3456106/91‑1 (MQ=255)
                                                                           cacgtatataGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCtt                  >  1:1643862/5‑91 (MQ=255)
                                                                           cACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCtt                  <  1:1405369/91‑1 (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:1522985/91‑1 (MQ=255)
                                                                                atataCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACtt              >  1:657/1‑91 (MQ=255)
                                                                                  ataCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACtttt            >  1:3900610/1‑91 (MQ=255)
                                                                                  ataCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACtttt            >  1:2558803/1‑91 (MQ=255)
                                                                                   taCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTg           <  1:1837645/91‑1 (MQ=255)
                                                                                   taCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTg           <  1:865926/91‑1 (MQ=255)
                                                                                   taCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTg           <  1:3054195/91‑1 (MQ=255)
                                                                                   taCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTg           <  1:3054920/91‑1 (MQ=255)
                                                                                    aCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGt          >  1:2820790/1‑91 (MQ=255)
                                                                                    aCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGt          <  1:3250173/91‑1 (MQ=255)
                                                                                    aCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGt          <  1:355113/91‑1 (MQ=255)
                                                                                    aCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGt          <  1:871719/91‑1 (MQ=255)
                                                                                       tGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCt       >  1:1339827/1‑91 (MQ=255)
                                                                                       tGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCt       >  1:1339798/1‑91 (MQ=255)
                                                                                           tgttgtGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCttttt   <  1:3758500/91‑1 (MQ=255)
                                                                                           tgttgtGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCttttt   <  1:3758606/91‑1 (MQ=255)
                                                                                            |                                                                                          
ATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTT  >  USA300TCH1516_ALE/2309717‑2309897

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: