Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 6,838 A→G 5.9% E593G (GAA→GGA)  gyrB → DNA gyrase subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE6,8380AG 5.9% 200.1 / 3.7 68E593G (GAA→GGA) gyrBDNA gyrase subunit B
Reads supporting (aligned to +/- strand):  ref base A (36/28);  new base G (2/2);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

ACTGACACAAGGTAAACAAAAGTATTATGTATACAATGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCGTAGACAATTTATAGAAGATAATGCAGTTTA  >  USA300TCH1516_ALE/6698‑6976
                                                                                                                                            |                                                                                                                                          
aCTGACACAAGGTAAACAAAAGTATTATGTATACAATGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGa                                                                                                                                            <  1:6330/141‑1 (MQ=255)
   tacacAAGGTAAACAAAAGTATTATGTATACAATGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAaac                                                                                                                                         <  2:37871/140‑1 (MQ=255)
   gACACAAGGTAAACAAAAGTATTATGTATACAATGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGaaaa                                                                                                                                         >  1:623747/1‑140 (MQ=255)
      acaAGGTAAACAAAAGGATTATGTATACAATGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAaacaac                                                                                                                                      >  1:414516/1‑141 (MQ=255)
          ggTAAACAAAAGTATTATGTATAAAATGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATg                                                                                                                                  <  1:140192/141‑1 (MQ=255)
            tAAACAAAAGTATTATGTATACAATGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACAACAAAATGGTCTATTGCACGATACTAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGAAAACAACAATGaa                                                                                                                                >  2:306830/1‑141 (MQ=255)
                    aGTATTATGTATACAATGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGc                                                                                                                        >  1:628830/1‑141 (MQ=255)
                      tattatGTATACAATGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGGTACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCAc                                                                                                                      <  2:592833/141‑1 (MQ=255)
                        ttatGTATACAATGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGCTCAATGATGGGAAACAGCAATGAACCCTGAGCAcgg                                                                                                                    >  1:234437/1‑139 (MQ=255)
                         tatGTATACAATGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACcgc                                                                                                                   <  2:259088/141‑1 (MQ=255)
                          atGTATACAATGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACcgcg                                                                                                                  <  1:475965/141‑1 (MQ=255)
                            gTATACAATGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAATCCTGAGCACCgcgcc                                                                                                                >  1:580677/1‑140 (MQ=255)
                              ataCAATGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGctct                                                                                                              <  1:123969/141‑1 (MQ=255)
                              ataCAATGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGctct                                                                                                              <  1:123970/141‑1 (MQ=255)
                                  aaTGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTa                                                                                                          <  2:508953/141‑1 (MQ=255)
                                  aaTGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTa                                                                                                          <  1:650523/141‑1 (MQ=255)
                                    tGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACa                                                                                                        >  1:548241/1‑141 (MQ=255)
                                     gATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACaa                                                                                                       >  1:685363/1‑141 (MQ=255)
                                        aGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGCTACAAAGGTCTTGGAGAAATGAATGCAGGTCAACGATGGGAAGCAACAATGAACCCTGAGCACCGCGCTCTTGCACAAGTa                                                                                                    >  1:300891/1‑141 (MQ=255)
                                        aGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACGACAATGAACCCTGAGCACCGCGCCCTTTTACAAGTa                                                                                                    >  1:190749/1‑141 (MQ=255)
                                         gggAACTTGATAAACTTAAATCTGAATTGAATCCAGCACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACGATGAACCGTGAGCGCCGCGCTCTTTTACAAGTaa                                                                                                   >  1:455031/1‑141 (MQ=255)
                                                  aTAAACTTAAATCTGAGTTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCGGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTGAAACTTGAAg                                                                                          >  1:346366/1‑141 (MQ=255)
                                                  aTAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCACTTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAg                                                                                          >  2:129719/1‑141 (MQ=255)
                                                  aTAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAg                                                                                          >  2:136296/1‑141 (MQ=255)
                                                   tAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGa                                                                                         >  1:581676/1‑141 (MQ=255)
                                                     aaCTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATg                                                                                       >  1:627467/1‑141 (MQ=255)
                                                       cTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCg                                                                                     <  1:122153/141‑1 (MQ=255)
                                                        ttAAATCTGAATTGAATCCAACACCAATATGGTCGTTTGCACGGTACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCTCCGCGCTCTTTTACAAGTATAACTTGAAGATGCGa                                                                                    <  2:494253/141‑1 (MQ=255)
                                                           aaTCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTg                                                                                 <  1:342647/141‑1 (MQ=255)
                                                              cTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAg                                                                              <  1:557172/141‑1 (MQ=255)
                                                                 aattgaatCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGGAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCgg                                                                           >  1:62277/1‑141 (MQ=255)
                                                                 aattgaatCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGGAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCgg                                                                           >  1:98000/1‑141 (MQ=255)
                                                                 aattgaatCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCgg                                                                           >  1:169642/1‑141 (MQ=255)
                                                                 aattgaatCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCgg                                                                           >  2:568426/1‑141 (MQ=255)
                                                                 aattgaatCCAACACCAAAATGGTCTATTGCACGATACAAAGGACTTGGAGAAATGAATGCAGATCAGTTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTCGAAGATGCGATTGAAGCgg                                                                           >  2:405577/1‑141 (MQ=255)
                                                                    tgaatCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGAcc                                                                        <  2:271641/141‑1 (MQ=255)
                                                                     gaatCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAATAAAACTTGAAGATGCGATTGAAGCGGACCa                                                                       >  2:280363/1‑141 (MQ=255)
                                                                           aaCACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACAtt                                                                 <  1:58419/141‑1 (MQ=255)
                                                                           aaCACCAAAATGGTCAATCGCACGATACAAAAGAATAAGAAACAAGACTGCAAAAAAATTAAAGAAAACTAGAAAGAGGCCTCAGTATCGCTCTCGTTTTCAAGTGAGGCGTGGAGATTTGATTGAAGCGGACCAAACAtt                                                                 <  2:683837/141‑1 (MQ=255)
                                                                                cAAAATGGTCTATTACACGATACAATGGTCTTGGAGAAATGAATACAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGaaa                                                            >  1:381967/1‑141 (MQ=255)
                                                                                   aaTGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTACAACTTGACGGTGCGGTTGAAGCTGACCTAACGAGTGAAATGt                                                         >  2:652753/1‑141 (MQ=255)
                                                                                      ggTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTaa                                                      >  2:618619/1‑141 (MQ=255)
                                                                                       gTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGACACGACCATGAACCCTGAGCACCGCGCTCTTTTAACAGTCTTACACGTAGAGGCGATTGAAGCGGACCAAACACTTGAACTGTTAAt                                                     >  2:244564/1‑141 (MQ=255)
                                                                                       gTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAAt                                                     >  1:577775/1‑141 (MQ=255)
                                                                                       gTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAAt                                                     >  1:629726/1‑141 (MQ=255)
                                                                                         cTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTTAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATgg                                                   <  2:258100/141‑1 (MQ=255)
                                                                                          taatGCACGATACAAAGGTCATGGAGAAATGAATGCAGATCAATTATGGGAAACAACAACGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTCAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATggg                                                  <  2:4204/138‑1 (MQ=255)
                                                                                          tATTGCACAATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATggg                                                  <  2:523491/141‑1 (MQ=255)
                                                                                            ttGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTg                                                >  2:273868/1‑141 (MQ=255)
                                                                                                      aaaaGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGAGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGa                                      <  1:306830/140‑1 (MQ=255)
                                                                                                            tccttgagaACTGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCg                                <  2:515684/136‑1 (MQ=255)
                                                                                                            tCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCg                                <  1:531701/141‑1 (MQ=255)
                                                                                                               tGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAGAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCGTAg                             <  2:685363/141‑1 (MQ=255)
                                                                                                                ggAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCGTAGa                            <  1:393017/141‑1 (MQ=255)
                                                                                                                  agaAAAGAACGCACAACAATTAAGGGAAACAACAATGAACCCTGAGAACGACTTTCATATACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCGTAGACa                          <  1:19891/141‑1 (MQ=255)
                                                                                                                   gaAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTCGAAAACCGTAGACaa                         >  2:541754/1‑141 (MQ=255)
                                                                                                                   gaAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCGTAGACaa                         >  2:399229/1‑141 (MQ=255)
                                                                                                                   gaAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCGTAGACaa                         >  2:541795/1‑141 (MQ=255)
                                                                                                                   gaAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCGTAGACaa                         >  1:252154/1‑141 (MQ=255)
                                                                                                                    aaaTGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCGTAGACAAt                        <  2:468901/141‑1 (MQ=255)
                                                                                                                     aatgaatgCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGCTTGAAGCGGACCAACCATTTGAAAGATTAATGGGTGACGTTGTAGAAAACCGTAGACAAtt                       >  2:17020/1‑141 (MQ=255)
                                                                                                                      atgaatgCAGATCAATTATGGGGAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCGTAGACAAttt                      <  2:62277/141‑1 (MQ=255)
                                                                                                                      atgaatgCAGATCAATTATGGGGAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCGTAGACAAttt                      <  2:98000/141‑1 (MQ=255)
                                                                                                                         aatgCAGATCAATTATGGGAAACAGCAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACAGTTGAAAGGTTAATGGGTGCCGTGGGAGAAAACCGTAGCCAATTtgta                   >  2:386890/1‑141 (MQ=255)
                                                                                                                          atgCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCGTAGACAATTTATag                  <  1:632219/141‑1 (MQ=255)
                                                                                                                                 acAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCGTAGACAATTTATAGAAGATaa           <  2:582174/140‑1 (MQ=255)
                                                                                                                                  cAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCGTAGACAATTTATCGAAGATaag          >  1:310455/1‑140 (MQ=255)
                                                                                                                                   aaTTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCGTAGACAATTTATAGAAGATAATg         >  1:653605/1‑141 (MQ=255)
                                                                                                                                     ttATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTCATGGGTGACGTTGTAGAAAACCGTAGACAATTTATAGAAGATAATGCa       >  1:319756/1‑141 (MQ=255)
                                                                                                                                          ggAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCGTAGACAATTTATAGAAGATAATGAAGTTTa  >  2:295496/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                          
ACTGACACAAGGTAAACAAAAGTATTATGTATACAATGATAGGGAACTTGATAAACTTAAATCTGAATTGAATCCAACACCAAAATGGTCTATTGCACGATACAAAGGTCTTGGAGAAATGAATGCAGATCAATTATGGGAAACAACAATGAACCCTGAGCACCGCGCTCTTTTACAAGTAAAACTTGAAGATGCGATTGAAGCGGACCAAACATTTGAAATGTTAATGGGTGACGTTGTAGAAAACCGTAGACAATTTATAGAAGATAATGCAGTTTA  >  USA300TCH1516_ALE/6698‑6976

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