Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,077,629 C→A 8.5% M171I (ATG→ATT yxlF_3 ← putative ABC transporter ATP‑binding protein YxlF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,077,6290CA 8.5% 677.6 / 42.9 199M171I (ATG→ATTyxlF_3putative ABC transporter ATP‑binding protein YxlF
Reads supporting (aligned to +/- strand):  ref base C (87/95);  new base A (12/5);  total (99/100)
Fisher's exact test for biased strand distribution p-value = 8.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

Only 100 of 211 total aligned reads displayed.

ATCTTGAACAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATCCAAGATT  >  USA300TCH1516_ALE/2077539‑2077719
                                                                                          |                                                                                          
   ttGAACAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCata                                                                                         >  2:2703696/1‑91 (MQ=255)
   ttGAACAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCata                                                                                         >  2:142903/1‑91 (MQ=255)
    tGAACAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCatat                                                                                        >  2:864637/1‑91 (MQ=255)
    tGAACAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCatat                                                                                        <  1:902734/91‑1 (MQ=255)
     gAACAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATc                                                                                       <  2:1180336/91‑1 (MQ=255)
      aaCAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCa                                                                                      <  1:1559313/91‑1 (MQ=255)
       aCAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATGCTCATATCaa                                                                                     <  2:3527672/91‑1 (MQ=255)
       aCAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATGAGAATTCTAATATCaa                                                                                     <  2:1053748/91‑1 (MQ=255)
        cAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAg                                                                                    <  1:209894/91‑1 (MQ=255)
         aaaaTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGt                                                                                   <  1:323859/91‑1 (MQ=255)
         aaaaTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGt                                                                                   >  2:2634214/1‑91 (MQ=255)
          aaaTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGtt                                                                                  >  2:663777/1‑91 (MQ=255)
          aaaTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGtt                                                                                  >  1:1236689/1‑91 (MQ=255)
            aTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCa                                                                                <  1:687961/91‑1 (MQ=255)
            aTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCa                                                                                <  2:3213998/91‑1 (MQ=255)
            aTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCa                                                                                <  2:3330147/91‑1 (MQ=255)
            aTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATGAGAATTCTAATATc                                                                                       >  1:3524843/1‑84 (MQ=255)
              gTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCAtt                                                                              >  1:538868/1‑91 (MQ=255)
                ccGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTa                                                                            <  2:1601764/91‑1 (MQ=255)
                ccGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTa                                                                            >  1:1144571/1‑91 (MQ=255)
                 cGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTAc                                                                           >  2:896794/1‑91 (MQ=255)
                 cGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTAc                                                                           >  1:1379041/1‑91 (MQ=255)
                 cGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTAc                                                                           >  1:1378977/1‑91 (MQ=255)
                 cGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTAc                                                                           >  1:1378946/1‑91 (MQ=255)
                  gTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACt                                                                          >  1:1997816/1‑91 (MQ=255)
                  gTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACt                                                                          >  1:52091/1‑91 (MQ=255)
                   tctcTTAAAAATACAGCTCTATCACGAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTa                                                                         <  2:1496935/91‑1 (MQ=255)
                   tctcTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTa                                                                         >  1:2378040/1‑91 (MQ=255)
                     tctTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTaaa                                                                       <  2:3093785/91‑1 (MQ=255)
                     tctTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTaaa                                                                       <  2:2241165/91‑1 (MQ=255)
                     tctTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTaaa                                                                       <  2:826652/91‑1 (MQ=255)
                      ctTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAg                                                                      <  1:456641/91‑1 (MQ=255)
                      ctTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAg                                                                      <  1:3765624/91‑1 (MQ=255)
                       ttAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGa                                                                     <  1:1529108/91‑1 (MQ=255)
                        tAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGAc                                                                    <  1:2548945/91‑1 (MQ=255)
                          aaaaTACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACtt                                                                  >  2:2503178/1‑91 (MQ=255)
                             aTACAGCTCTATCACAAGTTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAAt                                                               <  2:2316618/91‑1 (MQ=255)
                                 aGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTa                                                           >  1:1195727/1‑91 (MQ=255)
                                  gCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAg                                                          <  2:1756498/91‑1 (MQ=255)
                                   ctctATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGt                                                         <  2:2202718/91‑1 (MQ=255)
                                       aTCACAAATTAATTCAATATCTTCTAACTTATGACTGGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGc                                                     <  2:3508583/91‑1 (MQ=255)
                                         cacaACTTAATTCAATATCTTCTAACTTATGACTTGATATTACAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGcac                                                   >  2:2026095/1‑91 (MQ=255)
                                          acaAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGcaca                                                  >  1:977743/1‑91 (MQ=255)
                                            aaaTTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATc                                                >  1:1104576/1‑91 (MQ=255)
                                              attaattCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCaa                                              <  2:1510182/91‑1 (MQ=255)
                                                taattCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAAtt                                            >  2:3026090/1‑91 (MQ=255)
                                                 aattCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTg                                           >  1:1903982/1‑91 (MQ=255)
                                                 aattCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTg                                           <  1:2650317/91‑1 (MQ=255)
                                                 aattCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTg                                           >  2:1212451/1‑91 (MQ=255)
                                                   ttCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAg                                         >  2:3321600/1‑91 (MQ=255)
                                                    tCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGc                                        >  2:1526122/1‑91 (MQ=255)
                                                    tCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGc                                        <  2:2354099/91‑1 (MQ=255)
                                                    tCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGc                                        <  2:1827382/91‑1 (MQ=255)
                                                     cAATATCTTCTGACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGcc                                       <  2:2958769/91‑1 (MQ=255)
                                                     cAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGcc                                       >  2:1140209/1‑91 (MQ=255)
                                                     cAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGcc                                       <  2:501904/91‑1 (MQ=255)
                                                     cAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGcc                                       >  1:3174164/1‑91 (MQ=255)
                                                     cAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGcc                                       <  2:1293472/91‑1 (MQ=255)
                                                      aaTATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCa                                      <  1:72082/91‑1 (MQ=255)
                                                      aaTATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCa                                      <  2:1901775/91‑1 (MQ=255)
                                                        tatCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATc                                    <  2:2330498/91‑1 (MQ=255)
                                                          tcttctAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTg                                  <  1:2793035/91‑1 (MQ=255)
                                                           cttctAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTgg                                 <  1:2634768/91‑1 (MQ=255)
                                                               tAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATcc                             >  1:2258013/1‑91 (MQ=255)
                                                               tAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATcc                             >  1:1977214/1‑91 (MQ=255)
                                                               tAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATcc                             >  2:438940/1‑91 (MQ=255)
                                                               tAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATcc                             >  1:1766032/1‑91 (MQ=255)
                                                               tAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATcc                             <  2:649002/91‑1 (MQ=255)
                                                                aaCTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCa                            <  1:1005935/91‑1 (MQ=255)
                                                                aaCTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCa                            <  2:3777841/91‑1 (MQ=255)
                                                                aaCTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCa                            >  1:1735433/1‑91 (MQ=255)
                                                                aaCTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCa                            <  1:1616936/91‑1 (MQ=255)
                                                                 aCTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCAt                           <  2:2314169/91‑1 (MQ=255)
                                                                   ttATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATAc                         <  2:171138/91‑1 (MQ=255)
                                                                   ttATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATAc                         <  1:36263/91‑1 (MQ=255)
                                                                   ttATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATAc                         <  1:3662463/91‑1 (MQ=255)
                                                                    tATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATAcc                        <  1:3355763/91‑1 (MQ=255)
                                                                     aTGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCa                       <  2:1008897/91‑1 (MQ=255)
                                                                     aTGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCa                       <  1:3476361/91‑1 (MQ=255)
                                                                      tGACTTGATATTAGAATTCTCATATCAAGTTCATTTATTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCAt                      >  1:3272422/1‑91 (MQ=255)
                                                                       gACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCAtt                     <  1:187939/91‑1 (MQ=255)
                                                                        aCTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCAttt                    >  2:3553898/1‑91 (MQ=255)
                                                                         cTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTg                   >  2:3340659/1‑91 (MQ=255)
                                                                          ttGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGt                  <  1:142903/91‑1 (MQ=255)
                                                                            gATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAg                >  1:2883531/1‑91 (MQ=255)
                                                                             atatGAGAATTCTAATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAgg               >  2:3468036/1‑91 (MQ=255)
                                                                              tatTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGc              <  2:625563/91‑1 (MQ=255)
                                                                                ttAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTc            >  1:2493771/1‑91 (MQ=255)
                                                                                tgaGAATTCTAATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTc            >  1:3379404/3‑91 (MQ=255)
                                                                                 tAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCa           >  1:3844386/1‑91 (MQ=255)
                                                                                 tAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCa           >  1:3845532/1‑91 (MQ=255)
                                                                                   gAATTCTAATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCACTTGAGCCATCTGGATCCATACCATTTGTAGGCTCATc         >  2:1763568/1‑91 (MQ=255)
                                                                                    aaTTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTATTCAGCACATCAATTGAGCCATCTGGCTCCATACCATTTGTAGGCTCATcc        >  2:3597041/1‑91 (MQ=255)
                                                                                    aaTTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATcc        <  1:511617/91‑1 (MQ=255)
                                                                                    aaTTCTAATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATcc        >  2:2880657/1‑91 (MQ=255)
                                                                                      ttCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATCCaa      >  2:273932/1‑91 (MQ=255)
                                                                                        ctcATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATCCAAGa    >  1:3449857/1‑91 (MQ=255)
                                                                                         tcATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATCCAAGAt   <  1:573414/91‑1 (MQ=255)
                                                                                         tcATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATCCAAGAt   <  1:1074248/91‑1 (MQ=255)
                                                                                          cATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATCCAAGAtt  >  2:3687750/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATCTTGAACAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATCCAAGATT  >  USA300TCH1516_ALE/2077539‑2077719

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: