Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,664,318 T→C 5.3% I30V (ATA→GTA)  pbpH ← Penicillin‑binding protein H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,664,3180TC 5.3% 286.3 / 4.9 95I30V (ATA→GTA) pbpHPenicillin‑binding protein H
Reads supporting (aligned to +/- strand):  ref base T (45/45);  new base C (2/3);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

TTCTGTCTAAAATACGACCTCTTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCACTTTAAATCATATATGAATTTA  >  USA300TCH1516_ALE/1664180‑1664456
                                                                                                                                            |                                                                                                                                          
ttCTGTCTAAAATACGACCTCTTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATACaa                                                                                                                                          >  2:818558/1‑141 (MQ=255)
 tCTGTCTAAAATACGACCTCTTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGCAAAAAAATAAAtata                                                                                                                                           >  3:307567/1‑141 (MQ=255)
 tCTGTCTAAAATACGACCTCTTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAAt                                                                                                                                         >  1:825112/1‑141 (MQ=255)
   tGTCTAAAATACGACCTCTTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATca                                                                                                                                       >  3:278484/1‑141 (MQ=255)
    gTCTAAAATACGACCTCTTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATACAATcac                                                                                                                                      >  2:422023/1‑141 (MQ=255)
       tAACATACGACCTCTTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACAcc                                                                                                                                   <  3:93081/141‑1 (MQ=255)
       tAAAATACGACCTCTTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACAcc                                                                                                                                   <  2:652602/141‑1 (MQ=255)
       tAAAATACGACCTCTTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACAcc                                                                                                                                   <  1:264760/141‑1 (MQ=255)
        aaaaTACGACCTCTTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCa                                                                                                                                  >  2:171684/1‑141 (MQ=255)
         aaaTACGCCCTCTTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCaa                                                                                                                                 >  3:359746/1‑141 (MQ=255)
              cGACCTCTTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAAtata                                                                                                                            >  4:64255/1‑141 (MQ=255)
              cGACCTCTTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACATAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAAtata                                                                                                                            >  1:584506/1‑141 (MQ=255)
               gACCTCTTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAACATAAATATAATAACACCAAATATaa                                                                                                                           >  2:1081940/1‑141 (MQ=255)
                aCCTCTTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATaaa                                                                                                                          >  4:232586/1‑141 (MQ=255)
                 ccTCTTGGCACAGGGTGGTTCACGGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAg                                                                                                                         <  2:267258/141‑1 (MQ=255)
                  ctctTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGt                                                                                                                        <  2:690182/141‑1 (MQ=255)
                          actgACTCTTTCACTGTAATGGTTTCATCTTTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AATAATAAATATAATCACTCCAAATATAGAGTTAATtctc                                                                                                                <  2:587102/138‑1 (MQ=255)
                          acaGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATtctc                                                                                                                <  2:87652/141‑1 (MQ=255)
                            aGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCtt                                                                                                              <  1:313933/141‑1 (MQ=255)
                              aCTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGt                                                                                                            >  1:356217/1‑141 (MQ=255)
                              aCTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGt                                                                                                            >  2:1049506/1‑141 (MQ=255)
                                tCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTa                                                                                                          >  4:418625/1‑141 (MQ=255)
                                 cATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTaa                                                                                                         <  3:277305/141‑1 (MQ=255)
                                 cATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTaa                                                                                                         <  1:127376/141‑1 (MQ=255)
                                   ttCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGCGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCCCACCAAATATAAAGTTAATTCTCTTGTTAAtt                                                                                                       >  2:460610/1‑141 (MQ=255)
                                   ttCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAAtt                                                                                                       <  1:692234/141‑1 (MQ=255)
                                   ttCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAAtt                                                                                                       >  1:713613/1‑141 (MQ=255)
                                    tCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTg                                                                                                      >  3:193577/1‑141 (MQ=255)
                                    tCACTCTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTg                                                                                                      <  2:1054823/141‑1 (MQ=255)
                                           aaTGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTg                                                                                               <  1:850167/141‑1 (MQ=255)
                                           aaTGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATACAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTg                                                                                               <  1:472758/141‑1 (MQ=255)
                                              gTTTTCATCATTTTTTATAATTTGTTTATAATGTGATCCTTGTGCGTTTTTTAAATAACCTAAACGTAGTACTAGTACTGC‑‑TAAAATAAATATAATAACACCAAATATAAAGTTAATTATCTAGTTAAATGAATTTTAAAc                                                                                            >  4:413024/1‑141 (MQ=255)
                                               ttttCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAGTTCTCTTGTTAATTGTATTTTGAACg                                                                                           >  2:857965/1‑141 (MQ=255)
                                                 ttCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGAt                                                                                         >  3:301530/1‑141 (MQ=255)
                                                  tcatcaTTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACTCCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGAtt                                                                                        <  2:801434/141‑1 (MQ=255)
                                                   catcatTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGAttt                                                                                       >  2:123086/1‑141 (MQ=255)
                                                   catcatTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACATAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAGTTCTCTTGTTAATTGTATTTTGAACGAttt                                                                                       >  1:662924/1‑141 (MQ=255)
                                                    atcgtTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTtc                                                                                      <  2:142749/141‑1 (MQ=255)
                                                    atcatTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTtc                                                                                      <  1:228368/141‑1 (MQ=255)
                                                      catTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTtcat                                                                                    <  1:622973/141‑1 (MQ=255)
                                                        ttttttATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTtcatca                                                                                  >  2:985219/1‑141 (MQ=255)
                                                          ttttATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCAtt                                                                                >  1:208060/1‑141 (MQ=255)
                                                          ttttATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAAGATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCAtt                                                                                >  2:574138/1‑141 (MQ=255)
                                                           tttATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCAttt                                                                               <  2:743122/141‑1 (MQ=255)
                                                            ttATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTg                                                                              >  4:451027/1‑141 (MQ=255)
                                                            ttATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTg                                                                              >  1:225/1‑141 (MQ=255)
                                                            ttATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTg                                                                              <  2:418398/141‑1 (MQ=255)
                                                              ataATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGAt                                                                            >  1:113538/1‑141 (MQ=255)
                                                               taATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGAtt                                                                           <  2:234813/141‑1 (MQ=255)
                                                                    tGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATACAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCt                                                                      <  1:818558/141‑1 (MQ=255)
                                                                      tttATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCGTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCttt                                                                    <  2:652169/141‑1 (MQ=255)
                                                                      tttATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTAATAGAAATGA‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCttt                                                                    <  2:525712/141‑1 (MQ=255)
                                                                       ttATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCtttt                                                                   >  1:891614/1‑141 (MQ=255)
                                                                        tataATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTa                                                                  >  4:335494/1‑141 (MQ=255)
                                                                         ataATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAg                                                                 <  2:982243/141‑1 (MQ=255)
                                                                         ataATGTGAGCCTTGTGCGATTTGTAAATAACCAAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAAATCTCATCTTAATTGTATTTTGAACGCTTTCATCATTTGATTTTTCTTTTAg                                                                 >  2:265110/1‑141 (MQ=255)
                                                                              ttgAGCCTTGTGCGATTTGTAAATAACTTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTTTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCttt                                                            <  2:1083530/140‑1 (MQ=255)
                                                                              gtgAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCttt                                                            <  1:750208/141‑1 (MQ=255)
                                                                                 aGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTa                                                         >  2:306667/1‑141 (MQ=255)
                                                                                  gCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTaa                                                        <  2:669902/141‑1 (MQ=255)
                                                                                  gCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTaa                                                        <  2:795867/141‑1 (MQ=255)
                                                                                   ccTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAAc                                                       >  4:292419/1‑141 (MQ=255)
                                                                                   ccTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGATAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAAc                                                       <  1:286255/141‑1 (MQ=255)
                                                                                       gtgCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTAtatttcgtgaaa                                                   >  2:39898/1‑134 (MQ=255)
                                                                                        tgCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAAc                                                  <  3:16978/141‑1 (MQ=255)
                                                                                        tgCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAAc                                                  <  1:171684/141‑1 (MQ=255)
                                                                                         gCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAAct                                                 >  1:105992/1‑141 (MQ=255)
                                                                                         gCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAAct                                                 >  1:186242/1‑141 (MQ=255)
                                                                                             tttGTAAATGACCTAGACGTAGGACAAGTACTGC‑‑AAAAATAAATATAACCACACCAAATAAAAAGTTAAATCACTAGTTAATTGTATTTTACACGATATCATCATATGATCTTTCTTTTAGTCTTTCTAACAAACctacca                                             >  2:380046/1‑140 (MQ=255)
                                                                                                 tAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTacctacct                                         <  1:279659/141‑1 (MQ=255)
                                                                                                 tAAATAACCTAAACGTAGTACTAATACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTacctacct                                         <  2:411811/141‑1 (MQ=255)
                                                                                                  aaaTAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCAAACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACctc                                        >  3:378445/1‑141 (MQ=255)
                                                                                                      aaCCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTAATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTAtt                                    >  2:1035515/1‑141 (MQ=255)
                                                                                                       aCCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTc                                   <  2:340788/141‑1 (MQ=255)
                                                                                                       aCCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTc                                   <  2:507228/141‑1 (MQ=255)
                                                                                                          tacaCGTAGTACTAGTACTGC‑‑AAACATAACTATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCaaa                                <  2:270405/138‑1 (MQ=255)
                                                                                                            aaCGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTTACAAAACTACCTACCTCTATTCAAAGt                              <  4:359746/141‑1 (MQ=255)
                                                                                                            aaCGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGt                              >  1:1074807/1‑141 (MQ=255)
                                                                                                            aaCGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGt                              <  4:101875/141‑1 (MQ=255)
                                                                                                                 agtcctagtactGC‑‑AAAAATAAATATAATCACACCAAATCTAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTc                         >  2:98559/1‑141 (MQ=255)
                                                                                                                 agtactagtactGCAAAAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCtt                           <  4:307567/141‑1 (MQ=255)
                                                                                                                  gtactagtactGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCa                        >  1:232068/1‑141 (MQ=255)
                                                                                                                  gtactagtactGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCa                        >  3:160272/1‑141 (MQ=255)
                                                                                                                   tactagtactGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCAc                       >  3:92279/1‑141 (MQ=255)
                                                                                                                   tactagtactGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCAc                       <  2:216195/141‑1 (MQ=255)
                                                                                                                    actagtactGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCACt                      <  1:531986/141‑1 (MQ=255)
                                                                                                                      tagtactGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCACttt                    >  1:279008/1‑141 (MQ=255)
                                                                                                                          actGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCACTTTAAAt                >  3:26887/1‑141 (MQ=255)
                                                                                                                              c‑‑AAAAATAAATACAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCACTTTAAATCata            <  3:82909/141‑1 (MQ=255)
                                                                                                                                aaacaataaataTAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCACTTTAAATCata            <  2:943602/137‑1 (MQ=255)
                                                                                                                                 aaaaaTAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACGTCTATTCAAAGTCTTTCACTTTAAATCatat           >  1:616445/1‑141 (MQ=255)
                                                                                                                                 aaaaaTAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCACTTTAAATCatat           <  3:281827/141‑1 (MQ=255)
                                                                                                                                  aaaaTAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCACTTTAAATCatatt          >  4:254262/1‑140 (MQ=255)
                                                                                                                                  aaaaTAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCACTTTAAATCatata          <  4:125216/141‑1 (MQ=255)
                                                                                                                                    aataaataTAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCACTTTAAATCATATATg        <  3:337138/141‑1 (MQ=255)
                                                                                                                                        aataTAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCACTTTAAATCATATATGAAtt    <  4:131756/141‑1 (MQ=255)
                                                                                                                                         atatAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTATTTTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCACTTTAAATCATATATGAAttt   >  2:893168/1‑141 (MQ=255)
                                                                                                                                         atatAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCACTTTAAATCATATATGAAttt   <  4:407557/141‑1 (MQ=255)
                                                                                                                                          tataATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCACTTTAAATCATATATGAATTTa  >  2:252565/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                          
TTCTGTCTAAAATACGACCTCTTGGCACAGACTCATTCACTGTAATGTTTTCATCATTTTTTATAATTTGTTTATAATGTGAGCCTTGTGCGATTTGTAAATAACCTAAACGTAGTACTAGTACTGC‑‑AAAAATAAATATAATCACACCAAATATAAAGTTAATTCTCTTGTTAATTGTATTTTGAACGATTTCATCATTTGATTTTTCTTTTAGTCTTTTTAACAAAACTACCTACCTCTATTCAAAGTCTTTCACTTTAAATCATATATGAATTTA  >  USA300TCH1516_ALE/1664180‑1664456

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: