Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,526,532 T→A 7.0% K471* (AAA→TAA)  dinG_1 ← putative ATP‑dependent helicase DinG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,526,5320TA 7.0% 487.0 / 21.8 142K471* (AAA→TAA) dinG_1putative ATP‑dependent helicase DinG
Reads supporting (aligned to +/- strand):  ref base T (54/78);  new base A (4/6);  total (58/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 146 total aligned reads displayed.

ATTTTTCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCAAACAA  >  USA300TCH1516_ALE/1526442‑1526622
                                                                                          |                                                                                          
aTTTTTCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATAttt                                                                                            >  1:2784206/1‑91 (MQ=255)
aTTTTTCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATAttt                                                                                            >  1:2944229/1‑91 (MQ=255)
  ttttCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTaa                                                                                          <  1:3385728/91‑1 (MQ=255)
   tttCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTAAat                                                                                         <  1:618494/91‑1 (MQ=255)
   tttCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTAAat                                                                                         <  1:2643712/91‑1 (MQ=255)
   tttCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTAAat                                                                                         <  1:2969535/91‑1 (MQ=255)
    ttCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAata                                                                                        <  1:2383896/91‑1 (MQ=255)
    tggaaaaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAata                                                                                        >  1:2256832/5‑91 (MQ=255)
     tCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAatat                                                                                       >  1:2485600/1‑91 (MQ=255)
     tCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTAAatat                                                                                       <  1:2132919/91‑1 (MQ=255)
     tCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCGGATATTAAatat                                                                                       <  1:3203097/91‑1 (MQ=255)
     ggaaaaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAatat                                                                                       >  1:3128366/4‑91 (MQ=255)
      ccAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATc                                                                                      <  1:301903/91‑1 (MQ=255)
      ccAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATc                                                                                      >  1:2592740/1‑91 (MQ=255)
        aaaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAg                                                                                    <  1:2144100/91‑1 (MQ=255)
         aaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGc                                                                                   <  1:1637202/91‑1 (MQ=255)
         aaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGc                                                                                   <  1:142633/91‑1 (MQ=255)
         aaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGc                                                                                   <  1:2353803/91‑1 (MQ=255)
          aaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCa                                                                                  <  1:194318/91‑1 (MQ=255)
             tCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaa                                                                               >  1:2680227/1‑91 (MQ=255)
             tCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaa                                                                               <  1:557821/91‑1 (MQ=255)
                ctggtTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATTTCAGCATaactc                                                                            <  1:3207495/91‑1 (MQ=255)
                ctgctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactc                                                                            >  1:2202174/1‑91 (MQ=255)
                ctgctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactc                                                                            <  1:1140109/91‑1 (MQ=255)
                ctgctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactc                                                                            >  1:2074296/1‑91 (MQ=255)
                ctgctTTTCCAGTTGATCAATTGCCTTTAACAAGTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactc                                                                            <  1:3401351/91‑1 (MQ=255)
                ctgctTTTCCAGTTGATCAATTGCCTTTAACAAGTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactc                                                                            <  1:3401242/91‑1 (MQ=255)
                 tgctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactca                                                                           >  1:685474/1‑91 (MQ=255)
                  gctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactcaa                                                                          >  1:4057711/1‑91 (MQ=255)
                   ctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactcaac                                                                         <  1:1373788/91‑1 (MQ=255)
                   ctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactcaac                                                                         <  1:1460768/91‑1 (MQ=255)
                   ctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactcaac                                                                         <  1:1734246/91‑1 (MQ=255)
                     tttCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactcaactc                                                                       <  1:864405/91‑1 (MQ=255)
                     tttCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactcaactc                                                                       <  1:864404/91‑1 (MQ=255)
                     tttCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactcaactc                                                                       <  1:864302/91‑1 (MQ=255)
                       tCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCAt                                                                     <  1:3119988/91‑1 (MQ=255)
                       tCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCAt                                                                     <  1:3118796/91‑1 (MQ=255)
                        ccAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCAtt                                                                    >  1:3261901/1‑91 (MQ=255)
                        ccAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCAtt                                                                    >  1:2621972/1‑91 (MQ=255)
                           gTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGt                                                                 <  1:2123734/91‑1 (MQ=255)
                            ttGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTc                                                                <  1:221562/91‑1 (MQ=255)
                             tGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCa                                                               <  1:378027/91‑1 (MQ=255)
                              gATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCAc                                                              <  1:2715346/91‑1 (MQ=255)
                               aTCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACt                                                             >  1:1689460/1‑91 (MQ=255)
                                 cAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTg                                                           >  1:2219091/1‑91 (MQ=255)
                                  aaTTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGa                                                          <  1:2717521/91‑1 (MQ=255)
                                   aTTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGAt                                                         <  1:3562715/91‑1 (MQ=255)
                                    ttGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATAAAACCAAGCTGATCTTTAATATCAGCATAACTCAACTCATTTGTCACTTGAtt                                                        >  1:2014012/1‑91 (MQ=255)
                                     tGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGAttt                                                       >  1:1516739/1‑91 (MQ=255)
                                     tGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGAttt                                                       >  1:1517115/1‑91 (MQ=255)
                                       tcTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTaa                                                     >  1:1283865/2‑91 (MQ=255)
                                        cTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAg                                                    >  1:1892853/1‑91 (MQ=255)
                                        cTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAg                                                    <  1:3333310/91‑1 (MQ=255)
                                         tttAACAATTTTTCATTTTCGTTTTCACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAgg                                                   <  1:3870882/91‑1 (MQ=255)
                                          ttAGCAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTCAGGc                                                  <  1:3275965/91‑1 (MQ=255)
                                           tAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCa                                                 >  1:455868/1‑91 (MQ=255)
                                                aTTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTc                                            <  1:455038/91‑1 (MQ=255)
                                                  ttttCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTg                                          >  1:3466594/1‑91 (MQ=255)
                                                  ttttCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTg                                          >  1:3466404/1‑91 (MQ=255)
                                                   tttCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTgg                                         <  1:4102216/91‑1 (MQ=255)
                                                    ttCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGt                                        >  1:2169641/1‑91 (MQ=255)
                                                      cATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTaa                                      <  1:2329957/91‑1 (MQ=255)
                                                      cATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTaa                                      <  1:2329906/91‑1 (MQ=255)
                                                       aTTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAg                                     >  1:3569261/1‑91 (MQ=255)
                                                       aTTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAg                                     <  1:692002/91‑1 (MQ=255)
                                                       aTTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAg                                     >  1:3569728/1‑91 (MQ=255)
                                                        ttttCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGc                                    >  1:2489621/1‑91 (MQ=255)
                                                         tttCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCg                                   <  1:3888461/91‑1 (MQ=255)
                                                         tttCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCg                                   <  1:2924051/91‑1 (MQ=255)
                                                             gTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGAtgt                               >  1:2875264/1‑91 (MQ=255)
                                                                 ttACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCtt                           >  1:2532517/1‑91 (MQ=255)
                                                                 ttACCAATCAAAACAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCtt                           >  1:3196630/1‑91 (MQ=255)
                                                                   aCCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCa                         >  1:313511/1‑91 (MQ=255)
                                                                   aCCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCa                         >  1:4034892/1‑91 (MQ=255)
                                                                      aaTCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATcc                      <  1:2164149/91‑1 (MQ=255)
                                                                        tCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCAc                    <  1:1815597/91‑1 (MQ=255)
                                                                        tCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCAc                    <  1:3711711/91‑1 (MQ=255)
                                                                         cAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACt                   <  1:2287022/91‑1 (MQ=255)
                                                                         cAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACt                   <  1:2286417/91‑1 (MQ=255)
                                                                           aaCCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGGTTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACtat                 <  1:3592975/91‑1 (MQ=255)
                                                                            aCCAAGCTGATATTTAATATCAGCATAACTCAACTCAGTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACtata                >  1:3510297/1‑91 (MQ=255)
                                                                            aCCAAGCTGAAATTTAATATCAGCATAACTCAACTCATTCGTCACTGGATTAAAAGCTTAGTCTGGTAAGCGATATGCTTCATCTACtata                >  1:4063824/1‑91 (MQ=255)
                                                                              cAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACa              >  1:3836188/1‑91 (MQ=255)
                                                                              cAAGCTGATATTTAATATCAGCATAACTCAACTCATTCGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACa              >  1:3505793/1‑91 (MQ=255)
                                                                                aGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAat            <  1:1553566/91‑1 (MQ=255)
                                                                                aGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAat            <  1:1041588/91‑1 (MQ=255)
                                                                                  cTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatca          <  1:3746954/91‑1 (MQ=255)
                                                                                   tGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatcat         <  1:1288535/91‑1 (MQ=255)
                                                                                   tGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatcat         >  1:1198477/1‑91 (MQ=255)
                                                                                   tGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatcat         <  1:2524277/91‑1 (MQ=255)
                                                                                   ggATATTTAATATCAGCATAACTCCACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatcat         <  1:2084830/90‑1 (MQ=255)
                                                                                    gATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatcatc        >  1:4005118/1‑91 (MQ=255)
                                                                                    gATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatcatc        <  1:2813747/91‑1 (MQ=255)
                                                                                    gATATTAAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatcatc        >  1:1637504/1‑91 (MQ=255)
                                                                                    gATATTAAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatcatc        >  1:1637160/1‑91 (MQ=255)
                                                                                     atatTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAAtcatca       <  1:1662692/91‑1 (MQ=255)
                                                                                     atatTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAAtcatca       <  1:1663240/91‑1 (MQ=255)
                                                                                     atatTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAAtcatca       >  1:3833819/1‑91 (MQ=255)
                                                                                       atTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCaaa     <  1:2369287/91‑1 (MQ=255)
                                                                                         ttAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCAAACa   <  1:850294/91‑1 (MQ=255)
                                                                                          |                                                                                          
ATTTTTCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCAAACAA  >  USA300TCH1516_ALE/1526442‑1526622

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