Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,797,361 T→C 6.1% S503S (TCT→TCC)  USA300HOU_RS14370 → N‑acetylmuramoyl‑L‑alanine amidase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,797,3610TC 6.1% 202.3 / 4.1 73S503S (TCT→TCC) USA300HOU_RS14370N‑acetylmuramoyl‑L‑alanine amidase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (34/33);  new base C (3/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-01

AGAAATTAAACAGTATTATTAAACACTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACATAATTGGAATAATCGAGCTCAATACCGTGATTATC  >  USA300TCH1516_ALE/2797222‑2797494
                                                                                                                                           |                                                                                                                                     
aGAAATTACACAGTATTATTACACACTATCAATTAACTCAGTTGGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTg                                                                                                                                      <  2:741394/141‑1 (MQ=255)
       aaaCGGAATGATTAAACACTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTTTATCAAGGATTATGATGATTCATCAGATGAATTCAACCCTTTCCGTGAGGTATCTGATAGTAt                                                                                                                               <  2:676371/141‑1 (MQ=255)
       aaaCAGTATTATTAAACACTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTAt                                                                                                                               >  2:348890/1‑141 (MQ=255)
        aaCAGTATTATTAAACACTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGTATTCAAACCTTTCCGTGAGGTATCTGATAGTATg                                                                                                                              <  1:701302/141‑1 (MQ=255)
        aaCAGTATTATTAAACACTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATg                                                                                                                              >  1:519336/1‑141 (MQ=255)
         aCAGTATTATTAAACACTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGc                                                                                                                             <  2:748577/141‑1 (MQ=255)
         aCAGTATTATTAAACACTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGc                                                                                                                             >  1:562844/1‑141 (MQ=255)
         aCAGTATTATTAAACACTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGc                                                                                                                             >  1:266809/1‑141 (MQ=255)
             tattatTAAACACTATCAATTAACTCAGTTTTGCGATGAACGCATGCCAGATTTACATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCata                                                                                                                         <  2:232984/141‑1 (MQ=255)
             tattatTAAACACTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCata                                                                                                                         <  1:301407/141‑1 (MQ=255)
              attattAAACACTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCatat                                                                                                                        <  2:157163/141‑1 (MQ=255)
              attattAAACACTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATAACTCATCACATGAATTCATACCTTGCCTTGAGGTATCTGATAGTATGCCatat                                                                                                                        <  1:112401/141‑1 (MQ=255)
                 attaAACACTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCca                                                                                                                     >  2:778239/1‑141 (MQ=255)
                   taAACACTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCcaca                                                                                                                   <  1:623039/141‑1 (MQ=255)
                    aaaCACTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACAt                                                                                                                  <  1:641210/141‑1 (MQ=255)
                    aaaCACTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACAt                                                                                                                  >  1:463858/1‑141 (MQ=255)
                      acacTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCCTTCCGTGAGGTATCTGATAGTATGCCATATCCACATgg                                                                                                                <  1:700832/141‑1 (MQ=255)
                           aTCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAAt                                                                                                           >  1:542530/1‑141 (MQ=255)
                           aTCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAAt                                                                                                           >  2:238007/1‑141 (MQ=255)
                             cAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGt                                                                                                         <  1:673904/141‑1 (MQ=255)
                             cAATTAACTCAGATTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCCCATGGCCAATGt                                                                                                         >  1:409176/1‑141 (MQ=255)
                              aaTTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTa                                                                                                        >  1:520834/1‑141 (MQ=255)
                              aaTTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTa                                                                                                        <  2:544476/141‑1 (MQ=255)
                               actAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTAc                                                                                                       <  1:57304/139‑1 (MQ=255)
                               aTTTACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTAc                                                                                                       <  2:387342/141‑1 (MQ=255)
                                 tAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGGGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACtt                                                                                                     >  1:489213/1‑141 (MQ=255)
                                 tAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACtt                                                                                                     >  2:606942/1‑141 (MQ=255)
                                   aCTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCAGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTgg                                                                                                   <  2:515073/141‑1 (MQ=255)
                                        gTTTGACGATGAACCCATGCCAGATGTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAACCCGTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGt                                                                                              <  2:339930/141‑1 (MQ=255)
                                         tttGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAAAATTTCGGCGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGta                                                                                             >  1:149205/1‑141 (MQ=255)
                                            gACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCCGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGtatat                                                                                          >  1:718342/1‑141 (MQ=255)
                                            gACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCCGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGtatat                                                                                          >  1:487406/1‑141 (MQ=255)
                                            gACGATGAACACATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTACGTGAAGCATATGATGTTATGCCATATCCACAGAGCCAATGAACATGGAACGtagat                                                                                          >  1:299737/1‑141 (MQ=255)
                                                 tGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCCGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCg                                                                                     >  2:100141/1‑141 (MQ=255)
                                                  gAACGCATGACAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTGCGTCTATAACCGt                                                                                    >  2:507770/1‑141 (MQ=255)
                                                       cATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGaa                                                                               <  2:789874/141‑1 (MQ=255)
                                                         agCCAGATTCAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAAc                                                                             <  1:422089/140‑1 (MQ=255)
                                                              aaTTTAGATAAATATACACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAAttt                                                                        <  2:317476/140‑1 (MQ=255)
                                                                 ttAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATACCCCTATAAACCAATTTGGt                                                                     >  2:3355/1‑141 (MQ=255)
                                                                      tAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACatc                                                                >  1:170634/1‑141 (MQ=255)
                                                                               aCGTTCTATCAAGGATTATGATGATTCATCAGATTAATTCAAACCTTTCCGTGAGCACAGCCCGATTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAgg                                                       <  1:182564/141‑1 (MQ=255)
                                                                               aCGTTCTATCAAGGAGTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGCGGAATTTGTTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATACCCGTATGAAACAATTTGGTACATCTATCTCAgg                                                       <  1:182542/141‑1 (MQ=255)
                                                                                cGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAggg                                                      >  2:209285/1‑140 (MQ=255)
                                                                                cGTTCTATCAAGGATCATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGt                                                      <  2:647113/141‑1 (MQ=255)
                                                                                 gTTCTATCAAGGATTATGATGATTCATCAGATGAATTGAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATGCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTg                                                     >  2:588350/1‑141 (MQ=255)
                                                                                  ttCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATGGCCACATGGAGATTTGTTTAGGGACGGATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGa                                                    >  1:343288/1‑141 (MQ=255)
                                                                                      aTCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCAAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCGGGTGATTTa                                                >  2:109703/1‑141 (MQ=255)
                                                                                         aaGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAggg                                             >  2:677008/1‑140 (MQ=255)
                                                                                             aTTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATg                                         <  2:341690/141‑1 (MQ=255)
                                                                                               tATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGGATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGca                                       >  2:581144/1‑141 (MQ=255)
                                                                                                atgatgATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCCGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGcac                                      <  1:100141/141‑1 (MQ=255)
                                                                                                 tgatgaTTCATCAGATGAATTCAAACCTTTCCGTTAGGTATCTTATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGcaca                                     <  1:209217/141‑1 (MQ=255)
                                                                                                 tgatgaTTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGcaca                                     <  1:632754/141‑1 (MQ=255)
                                                                                                  tatgttTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTCGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACAt                                    <  2:255898/140‑1 (MQ=255)
                                                                                                  gatgatTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACAt                                    >  1:696427/1‑141 (MQ=255)
                                                                                                   atgatTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTGCTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATgacggta                                   >  2:274965/1‑135 (MQ=255)
                                                                                                         catcaGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCAGATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAACTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACATAATTGGa                             <  2:592604/141‑1 (MQ=255)
                                                                                                           tcaGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATCTAACCGTATGAACCAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACATAATTGGaat                           >  2:287544/1‑141 (MQ=255)
                                                                                                           tcaGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACATAATTGGaat                           >  1:221820/1‑141 (MQ=255)
                                                                                                           tcaGATGAATTCAAACCTTTCCGTGAGGTATATGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACTGGATGAAAAAATTTGGTACATCTATCTCAGGCGCTTTAGGTGATGCAAATAATAGGaat                           >  2:575960/1‑141 (MQ=255)
                                                                                                              gATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACATAATTGGaataat                        >  1:360928/1‑141 (MQ=255)
                                                                                                              gATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACATAATTGGaataat                        >  1:723573/1‑141 (MQ=255)
                                                                                                                    ttCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACCTGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACATAATTGGAATAATCGAGCt                  <  1:63048/141‑1 (MQ=255)
                                                                                                                    ttCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACATAATTGGAATAATCGAGCt                  <  1:680869/141‑1 (MQ=255)
                                                                                                                        aaCCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACATAATTGGAATAATCGAGCTCAAt              >  1:703802/1‑141 (MQ=255)
                                                                                                                        aaCCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGACACTGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGCTTTAGGTGATGCACATAATTGGAATAATCGAGCTCAAt              >  1:457491/1‑141 (MQ=255)
                                                                                                                          ccTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAAACGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGATGATGCACATAATGTGAATAATCGGACTCAATAc            >  2:587891/1‑141 (MQ=255)
                                                                                                                            tttCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACATAATTGGAATAATCGAGCTCAATACCg          <  2:519336/141‑1 (MQ=255)
                                                                                                                             ttCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACATAATTGGAATAATCGAGCTCAATACCGt         <  1:348890/141‑1 (MQ=255)
                                                                                                                             ttCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACATAATTGGAATAATCGAGCTCAATACCGt         <  2:779100/141‑1 (MQ=255)
                                                                                                                               ccGTGGGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACATAATTGGAATAATCGAGCTCAATACCGTGa       <  2:414496/141‑1 (MQ=255)
                                                                                                                                 gTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCCCAGGTGATTTAGGTGATGCACATAATTGGAATAATCGAGCTCAATACCGTGAtt     >  2:628511/1‑141 (MQ=255)
                                                                                                                                    aGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACATAATTGGAATAATCGAGCTCAATACCGTGATTATc  <  1:486064/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                     
AGAAATTAAACAGTATTATTAAACACTATCAATTAACTCAGTTTGACGATGAACGCATGCCAGATTTAGATAAATATGAACGTTCTATCAAGGATTATGATGATTCATCAGATGAATTCAAACCTTTCCGTGAGGTATCTGATAGTATGCCATATCCACATGGCCAATGTACTTGGTACGTATATAACCGTATGAAACAATTTGGTACATCTATCTCAGGTGATTTAGGTGATGCACATAATTGGAATAATCGAGCTCAATACCGTGATTATC  >  USA300TCH1516_ALE/2797222‑2797494

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