Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 79,619 T→A 6.9% intergenic (+228/‑32) USA300HOU_RS00355 → / → cmoB hypothetical protein/tRNA U34 carboxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE79,6190TA 6.9% 822.0 / 40.1 250intergenic (+228/‑32)USA300HOU_RS00355/cmoBhypothetical protein/tRNA U34 carboxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (108/122);  new base A (9/8);  total (117/133)
Fisher's exact test for biased strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01

Only 100 of 270 total aligned reads displayed.

GACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGC  >  USA300TCH1516_ALE/79529‑79709
                                                                                          |                                                                                          
gACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAatt                                                                                            <  1:675202/91‑1 (MQ=255)
gACCACAGCTTAATTCTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAACGTACATAGTGTAGTACTATTCAAAACGTAatt                                                                                            >  1:3370723/1‑91 (MQ=255)
   cacaGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAattatt                                                                                         <  1:1573066/91‑1 (MQ=255)
       gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTGCTATTCAAAACGTAATTATTACGa                                                                                     >  2:3680444/1‑91 (MQ=255)
       gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGa                                                                                     <  2:2006326/91‑1 (MQ=255)
        cTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGAt                                                                                    >  1:1368119/1‑91 (MQ=255)
         ttAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATTATTACGAtt                                                                                   >  2:3041167/1‑91 (MQ=255)
         ttAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGAtt                                                                                   <  1:2963623/91‑1 (MQ=255)
         ttAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGAtt                                                                                   <  2:3149693/91‑1 (MQ=255)
          tAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGtttt                                                                                  <  1:3115130/91‑1 (MQ=39)
          tAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGAttt                                                                                  <  2:1291472/91‑1 (MQ=255)
          tAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGAttt                                                                                  <  1:3737362/91‑1 (MQ=255)
           aaTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTg                                                                                 <  2:1702574/91‑1 (MQ=255)
           aaTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTg                                                                                 >  1:2740464/1‑91 (MQ=255)
           aaTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTg                                                                                 <  1:636973/91‑1 (MQ=255)
             tttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGaa                                                                               <  1:597347/91‑1 (MQ=38)
             tttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGaa                                                                               >  2:14704/1‑91 (MQ=255)
             tttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGaa                                                                               <  2:3761454/91‑1 (MQ=255)
              ttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAAt                                                                              >  2:3556118/1‑91 (MQ=255)
               tttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAAtt                                                                             >  1:518767/1‑91 (MQ=255)
                  ggTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAg                                                                          >  1:1037220/1‑91 (MQ=255)
                  ggTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAg                                                                          <  1:3359202/91‑1 (MQ=255)
                      ttcattcaTGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgag                                                                      >  2:1605781/1‑91 (MQ=255)
                      ttcattcaTGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgag                                                                      >  1:1991835/1‑91 (MQ=255)
                       tcattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAGTTATTACGATTTGAATTAAGCgagg                                                                     >  1:3686663/1‑91 (MQ=255)
                        cattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgagga                                                                    <  2:3627362/91‑1 (MQ=255)
                        cattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgagga                                                                    >  2:2173787/1‑91 (MQ=255)
                         attcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgaggag                                                                   <  2:379373/91‑1 (MQ=255)
                         attcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgaggag                                                                   <  1:1824052/91‑1 (MQ=255)
                            catGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaat                                                                >  2:633320/1‑91 (MQ=255)
                              tGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatga                                                              <  2:1155444/91‑1 (MQ=255)
                               gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                             >  1:3661186/1‑91 (MQ=255)
                               gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                             >  2:1901294/1‑91 (MQ=255)
                               gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                             <  2:2277959/91‑1 (MQ=255)
                                gctgtAAAAAGTCCTTAGTGCATAACGGGACATGGTGCAGTGCGATTCAAAGCGGAGTCATTACGATTTGAATTAAGCGAGGATgatgaaa                                                            <  2:780797/87‑1 (MQ=255)
                                 cTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaat                                                           <  1:1085476/91‑1 (MQ=255)
                                  tttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatg                                                          <  2:2704131/91‑1 (MQ=255)
                                  tttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatg                                                          >  2:3489806/1‑91 (MQ=255)
                                  tttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatg                                                          <  1:3103995/91‑1 (MQ=255)
                                   ttAAAAAATGCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatga                                                         <  1:3545176/91‑1 (MQ=255)
                                   ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatga                                                         >  1:681732/1‑91 (MQ=255)
                                     aaaaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACt                                                       <  2:3396540/91‑1 (MQ=255)
                                     aaaaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACt                                                       <  2:1008615/91‑1 (MQ=255)
                                         aaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGa                                                   >  1:2856798/1‑91 (MQ=255)
                                          aTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGAc                                                  <  1:2920646/91‑1 (MQ=255)
                                            ccTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACtt                                                >  1:1290145/1‑91 (MQ=255)
                                             cTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTa                                               <  1:3343525/91‑1 (MQ=255)
                                             cTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTa                                               >  1:2823453/1‑91 (MQ=255)
                                             cTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTa                                               <  1:1851799/91‑1 (MQ=255)
                                              ttATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTAt                                              >  1:3838771/1‑91 (MQ=255)
                                               tATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATg                                             <  2:2610979/91‑1 (MQ=255)
                                               tATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATg                                             >  1:2465645/1‑91 (MQ=255)
                                                aTTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGa                                            <  2:1736814/91‑1 (MQ=255)
                                                aTTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGa                                            <  2:2446623/91‑1 (MQ=255)
                                                 ttGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGAc                                           >  2:495149/1‑91 (MQ=255)
                                                 ttGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGAc                                           <  1:3015430/91‑1 (MQ=255)
                                                  tGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACg                                          <  2:1183251/91‑1 (MQ=255)
                                                        aaTGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtggt                                    >  2:563210/1‑91 (MQ=255)
                                                        aaTGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtggt                                    <  2:2167901/91‑1 (MQ=255)
                                                          tGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtggtgg                                  <  1:3644565/91‑1 (MQ=255)
                                                            tACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCa                                <  2:3828112/91‑1 (MQ=255)
                                                              cATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCaaa                              >  2:3880672/1‑91 (MQ=255)
                                                              cATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCaaa                              <  1:2839203/91‑1 (MQ=255)
                                                               aTAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCaaaa                             >  2:3591781/1‑91 (MQ=255)
                                                                  gtgtAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAgg                          <  2:1937245/91‑1 (MQ=255)
                                                                  gtgtAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAgg                          <  1:1204485/91‑1 (MQ=255)
                                                                   tgtAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGt                         >  2:2751645/1‑91 (MQ=255)
                                                                   tgtAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGt                         <  1:1499713/91‑1 (MQ=255)
                                                                    gtagtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTc                        <  2:700332/91‑1 (MQ=255)
                                                                    gtagtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTc                        <  2:97960/91‑1 (MQ=255)
                                                                     tagtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGCGGAGCATGAAATGACTAAGACTTATGACGTTTGGAGGCAAAAAGGTCa                       >  1:3952423/1‑91 (MQ=255)
                                                                     tagtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCa                       <  1:3323473/91‑1 (MQ=255)
                                                                      agtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCaa                      <  1:2114725/91‑1 (MQ=255)
                                                                      agtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCaa                      <  1:1767000/91‑1 (MQ=255)
                                                                      agtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCaa                      >  1:1629968/1‑91 (MQ=255)
                                                                       gtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAg                     <  2:2366123/91‑1 (MQ=255)
                                                                       gtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAg                     <  1:912577/91‑1 (MQ=255)
                                                                        taCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGa                    >  2:3358933/1‑91 (MQ=255)
                                                                        taCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGa                    <  1:785425/91‑1 (MQ=255)
                                                                        taCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGa                    >  1:1895110/1‑91 (MQ=255)
                                                                         aCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGaa                   <  1:1641413/91‑1 (MQ=255)
                                                                          cTATTCAAAACGTAATTATTACGATGTGAATTAAGCGAGGAGAATGAAATGACGAAGACGTATGACGTTTGGTGGCAAAAAGGTCAAGaaa                  >  2:2840589/1‑90 (MQ=255)
                                                                            aTTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCa                >  2:1486782/1‑91 (MQ=255)
                                                                            aTTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCa                >  2:1521051/1‑91 (MQ=255)
                                                                             ttCAAATCGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAg               >  1:3250321/1‑91 (MQ=38)
                                                                              tCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAga              >  1:3065152/1‑91 (MQ=255)
                                                                              tCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAga              >  1:2576690/1‑91 (MQ=255)
                                                                               cAAATCGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgat             >  2:593958/1‑91 (MQ=38)
                                                                               cAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgat             >  2:292555/1‑91 (MQ=255)
                                                                                 aatcGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatga           >  1:2274353/4‑91 (MQ=39)
                                                                                 aaaCGTAATTATTACGATTTGAATTAAGCGAGGAAAATGAAATGACTAGACTTATGACGAATGGTGGCAAAAAGGTCAAGAATCAgatgaa          >  2:3262848/1‑90 (MQ=255)
                                                                                  aaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgat          >  1:1417709/1‑91 (MQ=255)
                                                                                  aaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgat          <  1:2361593/91‑1 (MQ=255)
                                                                                  aaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgat          <  2:1663704/91‑1 (MQ=255)
                                                                                   aCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgatg         >  2:3698347/1‑91 (MQ=255)
                                                                                      tAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGata      >  1:1884929/1‑91 (MQ=255)
                                                                                       aaTTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGatat     <  1:2173787/91‑1 (MQ=255)
                                                                                       aaTTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGatat     >  2:929835/1‑91 (MQ=255)
                                                                                         ttattaCGATTTGAATGAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATgg   <  2:3285588/91‑1 (MQ=255)
                                                                                          tattaCGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGc  <  1:885629/91‑1 (MQ=255)
                                                                                          |                                                                                          
GACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGC  >  USA300TCH1516_ALE/79529‑79709

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: