Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,666,504 A→G 5.1% S56S (TCT→TCC znuB ← High‑affinity zinc uptake system membrane protein ZnuB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,666,5040AG 5.1% 338.9 / 5.3 112S56S (TCT→TCCznuBHigh‑affinity zinc uptake system membrane protein ZnuB
Reads supporting (aligned to +/- strand):  ref base A (57/49);  new base G (4/2);  total (61/51)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

TAGCAATTTCTTGGTAAGCAGTGTATGACGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAAGAGTACCTCATAAAATCA  >  USA300TCH1516_ALE/1666365‑1666639
                                                                                                                                           |                                                                                                                                       
tAGCAATTTCTTGGTAAGCAGTGTATGACGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAg                                                                                                                                        <  1:82279/141‑1 (MQ=255)
tAGCAATTTCTTGGTAAGCAGTGTATGACGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAg                                                                                                                                        >  1:313296/1‑141 (MQ=255)
tAGCAATTTCGTGGTAAGCAGTGTATGACGTTCTTAATATTTCAATTAGAAGCGCAACTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATCGTAAGTAAAAACATACCGAAAg                                                                                                                                        <  2:965073/141‑1 (MQ=255)
   cAATTTCTTGGTAAGCAGTGTATGACGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGata                                                                                                                                     >  3:71233/1‑141 (MQ=255)
     aTTTCTTGGTAAGCAGTGTATGACGTTCTTACATTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGAtata                                                                                                                                   >  1:303504/1‑141 (MQ=255)
      tttCTTGGTAAGCAGTGTATGACGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATac                                                                                                                                  >  1:460067/1‑141 (MQ=255)
      tttATTGGTAAGAAGTGTATGACGTTCTAAATTTTTCAATTAGAAACGCAACTACTATTGCAAATAATATTTCAAACAACATTGGATTAATAAATACTAGTGTTTTCATAATAGTAAGTAAAAACATACCGAAAGATATac                                                                                                                                  >  4:254054/1‑141 (MQ=255)
        tCTTGGTAAGCAGTGTATGACGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCTAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAATCATACCGAAGGATATacca                                                                                                                                <  4:221679/141‑1 (MQ=255)
         cTTGGTAAGCAGTGTATGACGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATaccac                                                                                                                               >  2:414901/1‑141 (MQ=255)
         cTTGGTAAGCAGTGTATGACGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATaccac                                                                                                                               <  1:766259/141‑1 (MQ=255)
            ggTAAGCAGTGTATGACGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTa                                                                                                                            >  2:277497/1‑141 (MQ=255)
             gTAAGCAGTGTATGACGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTaa                                                                                                                           >  1:1084874/1‑141 (MQ=255)
             gTAAGCAGTGTATGACGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTaa                                                                                                                           >  1:1063325/1‑141 (MQ=255)
                  cAGTGTATGACGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATGATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTa                                                                                                                      <  2:871097/141‑1 (MQ=255)
                  cAGTGTATGACGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTa                                                                                                                      <  1:348155/141‑1 (MQ=255)
                   aGTGTATGACGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTAc                                                                                                                     >  4:434669/1‑141 (MQ=255)
                         tGACGTTCTTAATTTTTCAATTAGAAGCGCCCCTACTATTGCGAATAAGATTCCAAACCACATGGGATTAAGAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACt                                                                                                               <  1:757100/141‑1 (MQ=255)
                         tGACGTTCTTAATTTCTCAATTAGAAGCGCACCTACTAGTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACt                                                                                                               <  2:724771/141‑1 (MQ=255)
                           aCGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTa                                                                                                             >  1:594706/1‑141 (MQ=255)
                           aCGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTa                                                                                                             >  3:244195/1‑141 (MQ=255)
                             gTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTaga                                                                                                           >  4:426575/1‑141 (MQ=255)
                                  tAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATc                                                                                                      >  3:80561/1‑141 (MQ=255)
                                   aaTTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCa                                                                                                     >  2:727359/1‑141 (MQ=255)
                                     tttttCAATTAGAAGCGCACCTACTATTGCTAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGc                                                                                                   <  2:227915/141‑1 (MQ=255)
                                     tttttCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGc                                                                                                   <  3:127859/141‑1 (MQ=255)
                                      ttttCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCGCATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTATAGctgatgct                                                                                                  <  3:144331/141‑8 (MQ=255)
                                      ttttCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCt                                                                                                  <  1:758143/141‑1 (MQ=255)
                                        ttCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAAATATACCACCTAAAGTTACATGACTTAGAGCATCAGCtat                                                                                                >  4:44720/1‑141 (MQ=255)
                                           aaTTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACAAAATTGGATTAAAACATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGAATTAGAGCCTCAGCTATAAg                                                                                             >  2:373080/1‑141 (MQ=255)
                                              tAGAAGCGCACCTCCTATTGCGAATAAGCTTCCAAACCACAATGGACTAATAAAGACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGa                                                                                          <  3:431468/141‑1 (MQ=255)
                                              tAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGa                                                                                          <  3:404825/141‑1 (MQ=255)
                                               tttttcgcACCTACTACTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTATTTAAAAACATACCGAAGGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGAt                                                                                         >  2:116423/6‑141 (MQ=255)
                                               aGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAGGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGAt                                                                                         >  2:116492/1‑141 (MQ=255)
                                               aGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAGGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGAt                                                                                         >  1:172343/1‑141 (MQ=255)
                                                 aaGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAg                                                                                       >  2:456993/1‑141 (MQ=255)
                                                   gcgcACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAgtc                                                                                     >  3:164892/1‑141 (MQ=255)
                                                    cgcACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATTgtcg                                                                                    >  3:294524/1‑141 (MQ=255)
                                                    cgcACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAgtcg                                                                                    <  1:1047047/141‑1 (MQ=255)
                                                    cgcACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAgtcg                                                                                    <  3:354501/141‑1 (MQ=255)
                                                      cACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAgtcgtc                                                                                  <  2:6759/141‑1 (MQ=255)
                                                          tactaTTGCGAATAAGATTCCAAACCACATTGGACTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTaac                                                                              <  1:1005681/141‑1 (MQ=255)
                                                               ttGCGAATAAGATTCCAAACCACATTATATTAATAATCTCTATTTTTTTCATATTAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGa                                                                         <  1:753707/141‑1 (MQ=255)
                                                                tGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGCAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTGTAAGTGATAGTCGTCTAACAACGaa                                                                        >  2:333111/1‑140 (MQ=255)
                                                                  cGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATaa                                                                      >  4:299943/1‑141 (MQ=255)
                                                                    aaTAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATaaaa                                                                    >  2:1009193/1‑141 (MQ=255)
                                                                     aTAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAg                                                                   >  1:461662/1‑141 (MQ=255)
                                                                         gATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCAcc                                                               <  2:68809/141‑1 (MQ=255)
                                                                          agtCCAAACCACATTGGATTAATAAATACTAGTGTTGTCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCg                                                              <  2:220616/139‑1 (MQ=255)
                                                                           ttCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGa                                                             >  2:611968/1‑141 (MQ=255)
                                                                           ttCCAAACCACACTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGa                                                             >  2:330915/1‑141 (MQ=255)
                                                                            tCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGAt                                                            <  4:406704/141‑1 (MQ=255)
                                                                             ccAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGAtt                                                           >  3:143810/1‑141 (MQ=255)
                                                                             ccAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGAtt                                                           >  1:1009662/1‑141 (MQ=255)
                                                                               aaaCCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCTCCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTag                                                         <  1:269414/141‑1 (MQ=255)
                                                                               aaaCCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTag                                                         <  4:48169/141‑1 (MQ=255)
                                                                                aaCCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTaga                                                        >  1:312784/1‑141 (MQ=255)
                                                                                   cacaTTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAGGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGgc                                                     <  1:266618/141‑1 (MQ=255)
                                                                                   cacaTTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAGGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGgc                                                     >  2:214242/1‑141 (MQ=255)
                                                                                   cacaTTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGgc                                                     <  4:164892/141‑1 (MQ=255)
                                                                                    acaTTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGgcg                                                    <  2:716917/141‑1 (MQ=255)
                                                                                    acaTTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGgcg                                                    <  1:498787/141‑1 (MQ=255)
                                                                                      aTTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGAGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCa                                                  <  2:1063028/141‑1 (MQ=255)
                                                                                       ttcgATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCaa                                                 <  4:392443/138‑1 (MQ=255)
                                                                                       ttGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCaa                                                 >  2:477616/1‑141 (MQ=255)
                                                                                        tGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAAt                                                >  2:496666/1‑141 (MQ=255)
                                                                                        tGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAAt                                                >  2:412256/1‑141 (MQ=255)
                                                                                        tGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAAt                                                >  1:984573/1‑141 (MQ=255)
                                                                                          gATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACTTTGCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATaa                                              >  3:331423/1‑141 (MQ=255)
                                                                                           aTTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATaaa                                             >  4:218535/1‑141 (MQ=255)
                                                                                           aTTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATaaa                                             <  2:22010/141‑1 (MQ=255)
                                                                                            ttAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATaaaa                                            <  2:391281/141‑1 (MQ=255)
                                                                                              aataaataCTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAAcc                                          <  1:980837/141‑1 (MQ=255)
                                                                                               ataaataCTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCt                                         <  2:190781/141‑1 (MQ=255)
                                                                                                taaataCTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAGGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTa                                        >  1:568948/1‑141 (MQ=255)
                                                                                                taaataCTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTa                                        <  1:766168/141‑1 (MQ=255)
                                                                                                taaataCTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTa                                        <  2:468764/141‑1 (MQ=255)
                                                                                                  aataCTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATc                                      <  1:169737/141‑1 (MQ=255)
                                                                                                   ataCTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCa                                     >  4:135097/1‑141 (MQ=255)
                                                                                                   ataCTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCACAATAAAACCTATCa                                     >  2:1024690/1‑141 (MQ=255)
                                                                                                   ataCTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCACAATAAAACCTATCa                                     >  2:1024674/1‑141 (MQ=255)
                                                                                                           gtTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCa                             <  1:1056075/141‑1 (MQ=255)
                                                                                                               gCATAATTGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTaa                         <  1:257828/141‑1 (MQ=255)
                                                                                                                cATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAAt                        <  1:1037366/141‑1 (MQ=255)
                                                                                                                  taataGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTa                      >  4:3802/1‑141 (MQ=255)
                                                                                                                     taGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAag                   >  4:282100/1‑141 (MQ=255)
                                                                                                                     taGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAag                   >  1:909747/1‑141 (MQ=255)
                                                                                                                      agtaagtaAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAaga                  <  3:329730/141‑1 (MQ=255)
                                                                                                                         aagtaaAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAAGAGTa               >  1:358054/1‑141 (MQ=255)
                                                                                                                          agtaaAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAAGAGTAc              >  2:6950/1‑141 (MQ=255)
                                                                                                                           gtaaAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAAGAGTAcc             >  1:529477/1‑141 (MQ=255)
                                                                                                                            taaAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAAGAGTACCt            >  1:748517/1‑141 (MQ=255)
                                                                                                                            taaAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGCCTAACAACGAAAAAAGCACCGATTAGAGGCGCAATACAACATATCACGATACCCCTAATTAAAGAGTACCt            >  4:131840/1‑141 (MQ=255)
                                                                                                                             aaaaaCATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAAGAGTACCTc           <  1:338134/141‑1 (MQ=255)
                                                                                                                              aaaaCATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAAGAGTACCTCa          >  1:794862/1‑141 (MQ=255)
                                                                                                                               aaaCATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAAGAGTACCTCAt         <  2:478729/141‑1 (MQ=255)
                                                                                                                                   aTACCGAAAGATATACCACCTAAAGTTATATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAAGAGTACCTCATaaaa     <  3:403327/141‑1 (MQ=255)
                                                                                                                                     aCCGAAAGATATCCCACTTAACGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCAAAAAACGATAAAAGCACAGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAAGAGTACCTCATAAAATc   <  4:428941/141‑1 (MQ=255)
                                                                                                                                     aCCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAAGAGTACCTCATAAAATc   <  3:334872/141‑1 (MQ=255)
                                                                                                                                      ccGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAAGAGTACCTCATAAAATCa  >  1:410548/1‑141 (MQ=255)
                                                                                                                                      ccGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAAGAGTACCTCATAAAATCa  >  3:314574/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                       
TAGCAATTTCTTGGTAAGCAGTGTATGACGTTCTTAATTTTTCAATTAGAAGCGCACCTACTATTGCGAATAAGATTCCAAACCACATTGGATTAATAAATACTAGTGTTGGCATAATAGTAAGTAAAAACATACCGAAAGATATACCACCTAAAGTTACATGACTTAGAGCATCAGCTATAAGTGATAGTCGTCTAACAACGATAAAAGCACCGATTAGAGGCGCAATAAAACCTATCAAGATACCACTAATTAAAGAGTACCTCATAAAATCA  >  USA300TCH1516_ALE/1666365‑1666639

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