Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,744,473 C→T A67T (GCT→ACT)  apt ← Adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,744,4730CT97.4% 532.8 / 1.6 156A67T (GCT→ACT) aptAdenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (83/68);  total (86/70)
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.78e-01

Only 100 of 156 total aligned reads displayed.

ATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATC  >  USA300TCH1516_ALE/1744386‑1744560
                                                                                       |                                                                                       
  tCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAg                                                                                    >  2:725367/1‑91 (MQ=255)
   cTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAgg                                                                                   <  1:767771/91‑1 (MQ=255)
     aGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGAc                                                                                 >  1:237094/1‑91 (MQ=255)
     aGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGAc                                                                                 >  1:2202616/1‑91 (MQ=255)
     aGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGAc                                                                                 <  2:1894515/91‑1 (MQ=255)
     aGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGAc                                                                                 <  2:1894035/91‑1 (MQ=255)
      ggTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACa                                                                                >  2:2001292/1‑91 (MQ=255)
      ggTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACa                                                                                >  2:2864108/1‑91 (MQ=255)
      ggTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACa                                                                                >  1:2798406/1‑91 (MQ=255)
      ggTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACa                                                                                >  2:2240775/1‑91 (MQ=255)
        tCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGc                                                                              <  2:446540/91‑1 (MQ=255)
        tCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGc                                                                              <  2:313512/91‑1 (MQ=255)
         cATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGcc                                                                             <  2:2705674/91‑1 (MQ=255)
          aTACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCa                                                                            >  2:1684798/1‑91 (MQ=255)
          aTACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCa                                                                            >  2:1606055/1‑91 (MQ=255)
           tACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCaa                                                                           <  2:844559/91‑1 (MQ=255)
             cTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATg                                                                         >  1:1258772/1‑91 (MQ=255)
              tCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGa                                                                        >  1:29880/1‑91 (MQ=255)
              tCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAAGAAGTTACAGGACAGCCAATGa                                                                        <  2:1548215/91‑1 (MQ=255)
                aTAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATa                                                                      >  1:1601406/1‑91 (MQ=255)
                 tAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATaa                                                                     >  2:1697978/1‑91 (MQ=255)
                  aaCGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATaaa                                                                    >  2:1630678/1‑91 (MQ=255)
                   aCGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAg                                                                   >  2:1154023/1‑91 (MQ=255)
                   aCGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAg                                                                   <  1:525434/91‑1 (MQ=255)
                     gAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGcc                                                                 <  1:102490/91‑1 (MQ=255)
                     gAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGcc                                                                 >  2:1336773/1‑91 (MQ=255)
                     gAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGcc                                                                 >  2:510269/1‑91 (MQ=255)
                     gAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGcc                                                                 <  1:2401784/91‑1 (MQ=255)
                     gAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGcc                                                                 <  1:2401942/91‑1 (MQ=255)
                       aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAc                                                               <  2:2260374/91‑1 (MQ=255)
                         gACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCAcgc                                                             <  2:3039402/91‑1 (MQ=255)
                         gACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCAcgc                                                             >  2:2831142/1‑91 (MQ=255)
                          aCTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCAcgcg                                                            <  2:1892809/91‑1 (MQ=255)
                           cTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACgcgc                                                           <  2:517462/91‑1 (MQ=255)
                           cTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACgcgc                                                           >  1:1007325/1‑91 (MQ=255)
                            ttCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCt                                                          >  1:2814751/1‑91 (MQ=255)
                            ttCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCt                                                          <  1:1630678/91‑1 (MQ=255)
                             tCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCtt                                                         >  2:1533340/1‑91 (MQ=255)
                             tCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCtt                                                         >  1:1489844/1‑91 (MQ=255)
                               aCGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCa                                                       >  2:847495/1‑91 (MQ=255)
                                 gAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAgg                                                     >  1:1578400/1‑91 (MQ=255)
                                 gAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAgg                                                     >  1:1046040/1‑91 (MQ=255)
                                 gAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAgg                                                     >  1:2675708/1‑91 (MQ=255)
                                   ggTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTc                                                   <  2:137619/91‑1 (MQ=255)
                                     tAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCt                                                 >  2:1217554/1‑91 (MQ=255)
                                      aaTTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTa                                                >  1:2414327/1‑91 (MQ=255)
                                      aaTTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTa                                                >  2:776389/1‑91 (MQ=255)
                                       aTTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTAc                                               >  1:898892/1‑91 (MQ=255)
                                        tttCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACa                                              >  1:2347531/1‑91 (MQ=255)
                                        tttCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACa                                              >  1:292729/1‑91 (MQ=255)
                                         ttCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACaa                                             <  2:3036587/91‑1 (MQ=255)
                                         ttCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACaa                                             <  1:2771469/91‑1 (MQ=255)
                                           cccTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACg                                           >  1:126723/1‑91 (MQ=255)
                                            ccTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGa                                          <  2:502779/91‑1 (MQ=255)
                                            ccTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGa                                          >  2:2765857/1‑91 (MQ=255)
                                               tcttTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGatat                                       >  2:1353083/1‑91 (MQ=255)
                                                cttTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                                      >  2:2157056/1‑91 (MQ=255)
                                                cttTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                                      <  1:1191687/91‑1 (MQ=255)
                                                 ttttCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCa                                     >  2:968819/1‑91 (MQ=255)
                                                   ttCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAAc                                   <  1:1397923/91‑1 (MQ=255)
                                                     cTAACAGGTGCAAAGCCAATCCCCATTGACTAAGGCCCAGGCCAGCCCATGATCACGCCCCGCGCTCCACAGCCTACACCGACGTCGACAt                                 >  2:1424117/1‑91 (MQ=255)
                                                     cTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAt                                 >  2:1424233/1‑91 (MQ=255)
                                                       aaCAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtct                               >  1:3076484/1‑91 (MQ=255)
                                                       aaCAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtct                               <  2:1756035/91‑1 (MQ=255)
                                                       aaCAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtct                               >  2:1472558/1‑91 (MQ=255)
                                                        aCAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtctc                              >  2:212193/1‑91 (MQ=255)
                                                        aCAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtctc                              >  1:38193/1‑91 (MQ=255)
                                                        aCAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtctc                              >  1:38121/1‑91 (MQ=255)
                                                         cAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATctct                             <  2:2348072/91‑1 (MQ=255)
                                                          aGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTg                            <  2:740855/91‑1 (MQ=255)
                                                          aGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTg                            <  1:554702/91‑1 (MQ=255)
                                                              gCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCtt                        >  1:2160741/1‑91 (MQ=255)
                                                              gCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCtt                        >  1:2163189/1‑91 (MQ=255)
                                                                aaaGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCtttt                      >  1:2493976/1‑91 (MQ=255)
                                                                aaaGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCtttt                      >  1:2494786/1‑91 (MQ=255)
                                                                 aaGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTg                     >  1:1353801/1‑91 (MQ=255)
                                                                 aaGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTg                     <  1:1687926/91‑1 (MQ=255)
                                                                 aaGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTg                     >  2:406034/1‑91 (MQ=255)
                                                                     cAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTa                 >  2:263577/1‑91 (MQ=255)
                                                                     cAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTa                 <  2:2531582/91‑1 (MQ=255)
                                                                     cAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTa                 >  1:835213/1‑91 (MQ=255)
                                                                     cAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTa                 >  1:1967002/1‑91 (MQ=255)
                                                                      aaTCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTAt                <  1:3058133/91‑1 (MQ=255)
                                                                        tCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTc              >  2:1330070/1‑91 (MQ=255)
                                                                        tCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTATGCGTATTc              >  1:2974111/1‑91 (MQ=255)
                                                                        tCCCCATAGAATAAGTTACAGGACAGCCAATGATAACGCCAAGCGCTTCAGGTCCTACAACGATATAGACATCTCTCTCTCTTGCGTATTc              >  2:2879972/1‑91 (MQ=255)
                                                                         ccccATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCt             <  2:326684/91‑1 (MQ=255)
                                                                         ccccATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCt             <  1:1922655/91‑1 (MQ=255)
                                                                         ccccATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCt             <  1:1922644/91‑1 (MQ=255)
                                                                          cccATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTa            >  2:103747/1‑91 (MQ=255)
                                                                          cccATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTa            <  2:119732/91‑1 (MQ=255)
                                                                          cccATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTa            >  2:104083/1‑91 (MQ=255)
                                                                          cccATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTa            <  1:1703303/91‑1 (MQ=255)
                                                                           ccATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTAc           <  2:21663/91‑1 (MQ=255)
                                                                             aTTGAATAAGTTACAGGACAGCCAATGATAATGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACaa         >  2:2370961/1‑91 (MQ=255)
                                                                             aTTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACaa         >  2:593437/1‑91 (MQ=255)
                                                                               tGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAtt       <  1:1572831/91‑1 (MQ=255)
                                                                               tGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAtt       >  2:1728162/1‑91 (MQ=255)
                                                                               tGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACACCTCTGTCTTTCGCGTATTCTACAAtt       >  2:2683137/1‑91 (MQ=255)
                                                                                gAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAttt      <  1:2422025/91‑1 (MQ=255)
                                                                                       |                                                                                       
ATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATC  >  USA300TCH1516_ALE/1744386‑1744560

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: