Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,701,267 C→T 5.0% V198M (GTG→ATG)  bacF ← Transaminase BacF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,701,2670CT 5.0% 664.8 / 13.8 197V198M (GTG→ATG) bacFTransaminase BacF
Reads supporting (aligned to +/- strand):  ref base C (103/84);  new base T (7/3);  total (110/87)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 202 total aligned reads displayed.

TTCAATCGCAACATCTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTGGATACGT  >  USA300TCH1516_ALE/2701178‑2701357
                                                                                         |                                                                                          
ttCAATCGCAACATCTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACg                                                                                           >  2:2289524/1‑91 (MQ=255)
ttCAATCGCAACATCTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACg                                                                                           >  2:3812539/1‑91 (MQ=255)
 tCAATCGCAACATCTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGa                                                                                          <  2:496776/91‑1 (MQ=255)
   aaTCGCAACATCTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCGTGCATGAtt                                                                                        <  1:3395708/91‑1 (MQ=255)
     tCGCAACATCTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGAtttt                                                                                      <  2:742416/91‑1 (MQ=255)
     tCGCAACATCTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGAtttt                                                                                      <  2:2700894/91‑1 (MQ=255)
      cGCAACATCTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTa                                                                                     >  2:1014054/1‑91 (MQ=255)
       gCAACATCTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAg                                                                                    <  1:1894561/91‑1 (MQ=255)
       gCAACATCTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAg                                                                                    <  1:3048125/91‑1 (MQ=255)
       gCAACATCTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAg                                                                                    >  1:1495654/1‑91 (MQ=255)
        cAACATCTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAgt                                                                                   >  1:1807735/1‑91 (MQ=255)
        cAACATCTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAgt                                                                                   >  2:3125996/1‑91 (MQ=255)
         aaCATCTTTCCCATTTTCCGAAGCAAGTATGCTAGGATTTTTAGCGGCGAAACCGAGAGCCAAATAAGCAAAATCAAGCACGATTTTAgtg                                                                                  <  2:397212/91‑1 (MQ=255)
         aaCATCTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAgtg                                                                                  >  2:563159/1‑91 (MQ=255)
         aaCATCTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTATAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAgtg                                                                                  >  2:3244574/1‑91 (MQ=255)
              cTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGt                                                                             <  1:2100651/91‑1 (MQ=255)
              cTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGt                                                                             >  2:669344/1‑91 (MQ=255)
              cTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGt                                                                             >  1:2638809/1‑91 (MQ=255)
                ttCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTAc                                                                           >  1:1799752/1‑91 (MQ=255)
                ttCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTAc                                                                           <  1:850332/91‑1 (MQ=255)
                 tCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTAcc                                                                          >  2:2846656/1‑91 (MQ=255)
                 tCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTAcc                                                                          >  2:1847201/1‑91 (MQ=255)
                  cccATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCt                                                                         <  1:1038504/91‑1 (MQ=255)
                  cccATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCt                                                                         <  2:615949/91‑1 (MQ=255)
                  cccATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCt                                                                         <  1:2279357/91‑1 (MQ=255)
                   ccATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCtt                                                                        <  1:2839627/91‑1 (MQ=255)
                   ccATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCtt                                                                        >  2:292621/1‑91 (MQ=255)
                   ccATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCtt                                                                        <  2:2377069/91‑1 (MQ=255)
                        ttCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAAt                                                                   >  1:2047323/1‑91 (MQ=255)
                         tCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCGTGCACGATTTTAGTGTCTGTACCTTTAAAtt                                                                  <  2:3021132/91‑1 (MQ=255)
                         tCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAAtt                                                                  >  1:2516269/1‑91 (MQ=255)
                           cGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTa                                                                <  2:3199105/91‑1 (MQ=255)
                            gAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAg                                                               <  2:2815486/91‑1 (MQ=255)
                             aaGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGc                                                              <  2:272943/91‑1 (MQ=255)
                              aGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCt                                                             >  1:2121971/1‑91 (MQ=255)
                               gCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTa                                                            >  1:740393/1‑91 (MQ=255)
                                cAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTAt                                                           >  1:39258/1‑91 (MQ=255)
                                cAAGTATACTAGGATTTGTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTAt                                                           <  1:3384994/91‑1 (MQ=255)
                                   gTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGc                                                        <  1:601066/91‑1 (MQ=255)
                                      taCTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTtc                                                     <  1:1237939/91‑1 (MQ=255)
                                        cTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTtcat                                                   <  1:139058/91‑1 (MQ=255)
                                        cTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTtcat                                                   <  1:3874688/91‑1 (MQ=255)
                                           ggATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCaa                                                >  1:943455/1‑91 (MQ=255)
                                               ttttAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAAc                                            >  1:2020433/1‑91 (MQ=255)
                                               ttttAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAAc                                            >  1:2417403/1‑91 (MQ=255)
                                                tttAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAAct                                           <  1:3638443/91‑1 (MQ=255)
                                                 ttAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAActt                                          >  2:3873383/1‑91 (MQ=255)
                                                 ttAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAActt                                          >  1:3626302/1‑91 (MQ=255)
                                                 ttAGCGGCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAActt                                          <  2:1077728/91‑1 (MQ=255)
                                                  tAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCCCGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAActtc                                         >  1:2604130/1‑91 (MQ=255)
                                                  tAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAActtc                                         >  1:2261523/1‑91 (MQ=255)
                                                  tAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAActtc                                         >  2:3744438/1‑91 (MQ=255)
                                                  tAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAActtc                                         >  2:2973/1‑91 (MQ=255)
                                                   aGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAActtct                                        >  2:1042505/1‑91 (MQ=255)
                                                    gCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAActtctt                                       >  2:3700033/1‑91 (MQ=255)
                                                     cGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCttt                                      <  1:1298319/91‑1 (MQ=255)
                                                      gTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTc                                     <  2:3505019/91‑1 (MQ=255)
                                                      gTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTc                                     <  2:1192479/91‑1 (MQ=255)
                                                      gTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTc                                     >  1:2511015/1‑91 (MQ=255)
                                                       tCGAAACAGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCg                                    <  1:1968984/91‑1 (MQ=255)
                                                          aaaCCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAg                                 <  2:3449309/91‑1 (MQ=255)
                                                             ccGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTg                              <  2:2216801/91‑1 (MQ=255)
                                                             ccGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTg                              >  1:1169363/1‑91 (MQ=255)
                                                               gAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTc                            <  2:2049552/91‑1 (MQ=255)
                                                                aaaGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAGTTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCg                           <  2:392419/91‑1 (MQ=255)
                                                                 aaGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGa                          <  2:3213016/91‑1 (MQ=255)
                                                                 aaGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGa                          <  1:711676/91‑1 (MQ=255)
                                                                  aGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGAt                         >  1:2894572/1‑91 (MQ=255)
                                                                  aGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGAt                         <  1:1162038/91‑1 (MQ=255)
                                                                   gCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATc                        >  2:3819035/1‑91 (MQ=255)
                                                                   gCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATc                        >  2:2593298/1‑91 (MQ=255)
                                                                   gCACCATAAGCAAAATCATGCACGATTTTAGAGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATc                        >  2:2546049/1‑91 (MQ=255)
                                                                    cACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATcc                       <  1:2946531/91‑1 (MQ=255)
                                                                     aCCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCa                      >  1:1870106/1‑91 (MQ=255)
                                                                      ccATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAg                     >  1:1596541/1‑91 (MQ=255)
                                                                        aTAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGtt                   >  1:3010983/1‑91 (MQ=255)
                                                                        aTAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGtt                   <  2:2069291/91‑1 (MQ=255)
                                                                        aTAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGtt                   <  2:2069000/91‑1 (MQ=255)
                                                                         tAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAGATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTg                  <  2:410719/91‑1 (MQ=255)
                                                                             cAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGact              >  2:1051390/1‑89 (MQ=255)
                                                                              aaaaTCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGatta             <  2:67339/91‑1 (MQ=255)
                                                                                aaTCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGattatt           <  1:2822479/91‑1 (MQ=255)
                                                                                aaTCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGattatt           >  2:3299060/1‑91 (MQ=255)
                                                                                 atcgtGCATGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTAttt          >  1:351004/5‑91 (MQ=255)
                                                                                   cgtGCATGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTgg        >  2:2096023/3‑91 (MQ=255)
                                                                                    gtGCATGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTGGa       >  2:1591054/2‑91 (MQ=255)
                                                                                    aTGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTGGa       <  2:150136/91‑1 (MQ=255)
                                                                                     tGCTCGATTTTAGTGTCTGTTCCTTTAAATTTAGCTATTGCTTTTTTAATACCTTCTTTCGTAGTTGTTGATCTAGTTGGATTATTTGGAt      >  2:1370731/1‑91 (MQ=255)
                                                                                     tGCATGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTGGAt      >  1:3345638/1‑91 (MQ=255)
                                                                                     tGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTGGAt      >  2:3693381/1‑91 (MQ=255)
                                                                                      gCATGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTGGATa     >  1:3543157/1‑91 (MQ=255)
                                                                                      gCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTGGATa     >  2:2317258/1‑91 (MQ=255)
                                                                                      gCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTGGATa     >  1:3747766/1‑91 (MQ=255)
                                                                                      gCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTGGATa     >  2:1262813/1‑91 (MQ=255)
                                                                                       cACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTGGATAc    >  1:761790/1‑91 (MQ=255)
                                                                                        aCGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCCCAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTGGATACg   >  2:162372/1‑91 (MQ=255)
                                                                                        aCGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTGGATACg   <  1:3940490/91‑1 (MQ=255)
                                                                                         cGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTGGATACGt  <  1:5001/91‑1 (MQ=255)
                                                                                         cGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTGGATACGt  <  2:3791245/91‑1 (MQ=255)
                                                                                         cGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTGGATACGt  <  1:1001360/91‑1 (MQ=255)
                                                                                         |                                                                                          
TTCAATCGCAACATCTTTCCCATTTTCCGAAGCAAGTATACTAGGATTTTTAGCGTCGAAACCGAAAGCACCATAAGCAAAATCATGCACGATTTTAGTGTCTGTACCTTTAAATTTAGCTATCGCTTCATCAAAAACTTCTTTCGTAGCTGTCGATCCAGTTGGATTATTTGGATACGT  >  USA300TCH1516_ALE/2701178‑2701357

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