Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,744,473 C→T 100% 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,4730CT99.5% 716.0 / ‑4.8 203A67T (GCT→ACT) aptAdenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (95/106);  total (96/107)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 203 total aligned reads displayed.

TATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTAT  >  USA300TCH1516_ALE/1744385‑1744559
                                                                                        |                                                                                      
tatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTa                                                                                      >  2:1869226/1‑91 (MQ=255)
  ttCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACa                                                                                    >  1:244509/1‑91 (MQ=255)
   tCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAg                                                                                   >  1:930409/1‑91 (MQ=255)
   tCTAGGTCATACGCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAg                                                                                   <  2:3233496/91‑1 (MQ=255)
    cTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAgg                                                                                  <  2:2036634/91‑1 (MQ=255)
      aGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGAc                                                                                <  2:3640575/91‑1 (MQ=255)
      aGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGAc                                                                                <  2:1773565/91‑1 (MQ=255)
      aGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGAc                                                                                <  2:2310036/91‑1 (MQ=255)
       ggTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACa                                                                               >  2:821242/1‑91 (MQ=255)
        gTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAg                                                                              <  1:3260695/91‑1 (MQ=255)
        gTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAg                                                                              <  2:2188212/91‑1 (MQ=255)
        gTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAg                                                                              <  1:2330470/91‑1 (MQ=255)
         tCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGc                                                                             >  1:617440/1‑91 (MQ=255)
          cATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGcc                                                                            <  1:74136/91‑1 (MQ=255)
          cATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGcc                                                                            >  1:1177070/1‑91 (MQ=255)
           aTACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCa                                                                           >  2:2453949/1‑91 (MQ=255)
            tACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCaa                                                                          <  1:3213699/91‑1 (MQ=255)
             aCTCATACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATg                                                                        >  2:3073112/1‑91 (MQ=255)
              cTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATg                                                                        >  2:2864072/1‑91 (MQ=255)
                 aTAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATa                                                                     <  1:1801663/91‑1 (MQ=255)
                 aTAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATa                                                                     >  2:2575816/1‑91 (MQ=255)
                 aTAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATa                                                                     <  2:2403109/91‑1 (MQ=255)
                  tAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATaa                                                                    <  1:2281663/91‑1 (MQ=255)
                  tAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATaa                                                                    >  2:2282754/1‑91 (MQ=255)
                    aCGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAg                                                                  >  1:3751175/1‑91 (MQ=255)
                     cGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGc                                                                 >  2:1939063/1‑91 (MQ=255)
                      gAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGcc                                                                <  2:1073946/91‑1 (MQ=255)
                       aaTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCa                                                               <  1:1995378/91‑1 (MQ=255)
                        aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCAc                                                              >  2:1670143/1‑91 (MQ=255)
                          gACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCAcgc                                                            >  2:1387686/1‑91 (MQ=255)
                             ttCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCt                                                         <  2:2095838/91‑1 (MQ=255)
                              tCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCtt                                                        >  2:493169/1‑91 (MQ=255)
                               cACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTc                                                       >  1:1734309/1‑91 (MQ=255)
                                 cGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAg                                                     <  1:712933/91‑1 (MQ=255)
                                  gAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAgg                                                    <  1:1183837/91‑1 (MQ=255)
                                  gAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAgg                                                    <  1:560505/91‑1 (MQ=255)
                                    ggTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTc                                                  <  1:3081044/91‑1 (MQ=255)
                                     gTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTcc                                                 >  2:739188/1‑91 (MQ=255)
                                     gTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTcc                                                 <  2:3465028/91‑1 (MQ=255)
                                       aaTTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTa                                               <  1:2629012/91‑1 (MQ=255)
                                       aaTTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTa                                               >  2:1929996/1‑91 (MQ=255)
                                       aaTTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTa                                               <  1:1705135/91‑1 (MQ=255)
                                        aTTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTAc                                              >  2:1822724/1‑91 (MQ=255)
                                        aTTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTAc                                              <  1:1667673/91‑1 (MQ=255)
                                        aTTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTAc                                              >  1:2477649/1‑91 (MQ=255)
                                          ttCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACaa                                            <  2:3038891/91‑1 (MQ=255)
                                          ttCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACaa                                            >  1:2695931/1‑91 (MQ=255)
                                            cccTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACg                                          >  1:656840/1‑91 (MQ=255)
                                             ccTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGa                                         <  2:609934/91‑1 (MQ=255)
                                                 cttTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                                     <  1:451973/91‑1 (MQ=255)
                                                   tttCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCaa                                   <  1:2386735/91‑1 (MQ=255)
                                                   tttCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCaa                                   <  2:331520/91‑1 (MQ=255)
                                                     tCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACa                                 >  1:2305838/1‑91 (MQ=255)
                                                     tCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACa                                 >  1:3068724/1‑91 (MQ=255)
                                                      cTACCAGGTGCAAAGCCAATCCCCATGGAGTAAGTTACAGAGCAGCCAATGATAAGGCCACGCGCTTCAGGACCTACAACGATATGAAccc                                >  2:178076/1‑89 (MQ=255)
                                                       tAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtc                               >  2:3621441/1‑91 (MQ=255)
                                                        caCAGGTGCAAAGCCAATCCCCATTGAATAAGATACAGGAAAGCCAATGATAAAGCCACGCTCTTCAGGTCCTACAACGATATCAACAtct                              <  2:1628679/90‑1 (MQ=255)
                                                        aaCAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtct                              >  2:2721732/1‑91 (MQ=255)
                                                        aaCAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtct                              <  1:496143/91‑1 (MQ=255)
                                                        aaCAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtct                              <  2:1744587/91‑1 (MQ=255)
                                                         aCAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtctc                             >  1:3116539/1‑91 (MQ=255)
                                                          cAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATctct                            <  2:1688742/91‑1 (MQ=255)
                                                          cAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATctct                            <  2:2749300/91‑1 (MQ=255)
                                                          cAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATctct                            <  1:821800/91‑1 (MQ=255)
                                                           aGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTg                           <  2:2714566/91‑1 (MQ=255)
                                                                cAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCttt                      <  2:761541/91‑1 (MQ=255)
                                                                 aaaGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCtttt                     >  1:1700854/1‑91 (MQ=255)
                                                                 aaaGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCtttt                     >  2:2252341/1‑91 (MQ=255)
                                                                  aaGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTg                    >  1:1584068/1‑91 (MQ=255)
                                                                  aaGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTg                    >  1:2808160/1‑91 (MQ=255)
                                                                  aaGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTg                    >  1:52330/1‑91 (MQ=255)
                                                                    gCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCg                  >  2:1408930/1‑91 (MQ=255)
                                                                      cAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTa                >  1:3644341/1‑91 (MQ=255)
                                                                      cAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTa                >  2:3475771/1‑91 (MQ=255)
                                                                       aaTCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTAt               <  1:2110753/91‑1 (MQ=255)
                                                                        aTCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTAtt              >  1:1163717/1‑91 (MQ=255)
                                                                        aTCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTAtt              >  2:560652/1‑91 (MQ=255)
                                                                         tCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTc             <  1:1300300/91‑1 (MQ=255)
                                                                         tCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTc             <  1:1768108/91‑1 (MQ=255)
                                                                         tCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTc             >  2:3536180/1‑91 (MQ=255)
                                                                         tCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTc             <  1:1366135/91‑1 (MQ=255)
                                                                          ccccATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCt            <  1:3501826/91‑1 (MQ=255)
                                                                          ccccATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCt            >  2:457675/1‑91 (MQ=255)
                                                                          ccccATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCt            >  2:457849/1‑91 (MQ=255)
                                                                           cccATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTa           <  1:2744864/91‑1 (MQ=255)
                                                                            ccATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTAc          >  1:1271935/1‑91 (MQ=255)
                                                                             cATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACa         <  2:1163717/91‑1 (MQ=255)
                                                                             cATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACa         <  2:431135/91‑1 (MQ=255)
                                                                               ttGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACact       >  2:3546592/1‑89 (MQ=255)
                                                                               ttGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAt       >  2:2542009/1‑91 (MQ=255)
                                                                                tGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTATAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAtt      <  2:3189097/91‑1 (MQ=255)
                                                                                tGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAtt      >  1:895658/1‑91 (MQ=255)
                                                                                tGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAtt      >  1:1550739/1‑91 (MQ=255)
                                                                                tGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTGCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAtt      <  2:1648226/91‑1 (MQ=255)
                                                                                 gAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAttt     <  1:1431318/91‑1 (MQ=255)
                                                                                  aaTAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAtttt    <  1:1236786/91‑1 (MQ=255)
                                                                                  aaTAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAtttt    >  1:3525763/1‑91 (MQ=255)
                                                                                   aTAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATctct                            <  1:1009411/66‑1 (MQ=255)
                                                                                   aTAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATctct                            >  2:1009411/1‑66 (MQ=255)
                                                                                    tAAGTTACAGGACAGCCAATGATGAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTAt  <  2:3792233/91‑1 (MQ=255)
                                                                                        |                                                                                      
TATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTAT  >  USA300TCH1516_ALE/1744385‑1744559

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: