Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,667,524 G→T 6.1% intergenic (‑35/+91) znuC_2 ← / ← nfo High‑affinity zinc uptake system ATP‑binding protein ZnuC/putative endonuclease 4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,667,5240GT 6.1% 920.7 / 30.4 263intergenic (‑35/+91)znuC_2/nfoHigh‑affinity zinc uptake system ATP‑binding protein ZnuC/putative endonuclease 4
Reads supporting (aligned to +/- strand):  ref base G (127/119);  new base T (2/15);  total (129/134)
Fisher's exact test for biased strand distribution p-value = 1.76e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 270 total aligned reads displayed.

TAGTAATTGACATTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAAGAGAGCATTATCATT  >  USA300TCH1516_ALE/1667434‑1667614
                                                                                          |                                                                                          
tagtaATTGACATTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAg                                                                                            >  1:2624575/1‑91 (MQ=255)
tagtaATTGACATTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAg                                                                                            >  1:236629/1‑91 (MQ=255)
 agtaATTGACATTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGt                                                                                           >  2:593463/1‑91 (MQ=255)
 agtaATTGACATTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGt                                                                                           <  1:694640/91‑1 (MQ=255)
   taATTGACATTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACATTac                                                                                         <  2:2672636/91‑2 (MQ=255)
   taATTGACATTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACATTac                                                                                         <  2:2894143/91‑2 (MQ=255)
      ttGACATTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGt                                                                                      <  1:1090889/91‑1 (MQ=255)
      ttGACATTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGt                                                                                      >  1:2670865/1‑91 (MQ=255)
        gACATTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGttt                                                                                    >  1:3604837/1‑91 (MQ=255)
         aCATTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACATTACTGTTTa                                                                                   <  1:130100/91‑1 (MQ=255)
         aCATTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTa                                                                                   >  1:871189/1‑91 (MQ=255)
          cATTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACATTACTGTTTAc                                                                                  <  1:109152/91‑1 (MQ=255)
          cATTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACATTACTGTTTAc                                                                                  <  1:1821528/91‑1 (MQ=255)
           aTTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACg                                                                                 <  1:3779560/91‑1 (MQ=255)
           aTTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACg                                                                                 <  2:1011337/91‑1 (MQ=255)
           aTTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACg                                                                                 <  2:1522948/91‑1 (MQ=255)
            tttttCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGa                                                                                <  2:3238640/91‑1 (MQ=255)
            tttttCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGa                                                                                >  2:559123/1‑91 (MQ=255)
              tttCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACt                                                                              >  2:2604241/1‑91 (MQ=255)
                 caattcaaAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCaa                                                                           <  1:1533211/91‑1 (MQ=255)
                     tcaaAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAAtt                                                                       <  1:1520241/91‑1 (MQ=255)
                      caaAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTa                                                                      <  2:2558580/91‑1 (MQ=255)
                       aaaGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAg                                                                     >  1:2141741/1‑91 (MQ=255)
                       aaaGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAg                                                                     >  1:1145923/1‑91 (MQ=255)
                         aGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTCGAc                                                                   >  1:3560402/1‑91 (MQ=255)
                         aGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGAc                                                                   <  1:3541683/91‑1 (MQ=255)
                           aCTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACaa                                                                 <  2:2817201/91‑1 (MQ=255)
                             tGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACaaaa                                                               <  1:2746387/91‑1 (MQ=255)
                             tGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACaaaa                                                               <  2:3073879/91‑1 (MQ=255)
                             tGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACaaaa                                                               <  1:1279685/91‑1 (MQ=255)
                              ggTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAAt                                                              >  2:322813/1‑91 (MQ=255)
                                  tCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAag                                                          >  2:3902441/1‑91 (MQ=255)
                                    aTCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAagag                                                        >  2:3160030/1‑91 (MQ=255)
                                     tCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAgaga                                                       >  2:1144335/1‑91 (MQ=255)
                                     tCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAgaga                                                       <  1:1385582/91‑1 (MQ=255)
                                      cGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGat                                                      <  2:1676829/91‑1 (MQ=255)
                                      cGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGat                                                      >  1:3826585/1‑91 (MQ=255)
                                        tATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGatat                                                    >  1:579406/1‑91 (MQ=255)
                                        tATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGatat                                                    >  2:3132372/1‑91 (MQ=255)
                                         aTTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACACAATAAAGAGATATg                                                   >  1:1406141/1‑91 (MQ=255)
                                         aTTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATg                                                   >  2:2938227/1‑91 (MQ=255)
                                         aTTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATg                                                   <  2:266279/91‑1 (MQ=255)
                                           tctcACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGcc                                                 >  1:1270194/1‑91 (MQ=255)
                                           tctcACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGcc                                                 >  1:294675/1‑91 (MQ=255)
                                           tctcACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGcc                                                 >  1:2513479/1‑91 (MQ=255)
                                           tctcACCTCGCATTCAACTATACAACACCTAGTAACATCTGTAAACAGTGATGTTTACGACTCAATATTAGACCAAATAAAGAGATATGcc                                                 <  2:1754332/91‑1 (MQ=255)
                                             tcACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGcccc                                               <  2:3828748/91‑1 (MQ=255)
                                             tcACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGcccc                                               >  2:2431516/1‑91 (MQ=255)
                                             tcACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGcccc                                               <  1:2044093/91‑1 (MQ=255)
                                             tcACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGACATGcccc                                               >  1:1957689/1‑91 (MQ=255)
                                              cACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGccccc                                              <  1:3591035/91‑1 (MQ=255)
                                              cACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGccccc                                              <  2:1831032/91‑1 (MQ=255)
                                              cACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGccccc                                              >  1:2915369/1‑91 (MQ=255)
                                               aCCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCt                                             >  1:879076/1‑91 (MQ=255)
                                               aCCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCt                                             <  2:1451278/91‑1 (MQ=255)
                                                 cTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTc                                           <  2:716915/91‑1 (MQ=255)
                                                 cTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTc                                           >  1:3288093/1‑91 (MQ=255)
                                                     cATTCAACTTTACAAGTCCTAGTAATGTATGTAATCCGTAATGTTTACGCCTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGt                                       <  2:3427409/91‑1 (MQ=255)
                                                     cATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGt                                       <  2:3359038/91‑1 (MQ=255)
                                                     cATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGt                                       <  2:2175924/91‑1 (MQ=255)
                                                      aTTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGtt                                      >  2:136435/1‑91 (MQ=255)
                                                      aTTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGACATATGCCCCCTTCAAGtt                                      <  2:2123800/91‑1 (MQ=255)
                                                       ttCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAGATTAGACAAAATAAAGAGATATGCCCCCTTCAAGttt                                     <  1:3351492/91‑1 (MQ=255)
                                                        tCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGtttt                                    >  2:3719114/1‑91 (MQ=255)
                                                          aaCTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTtttat                                  >  1:3724915/1‑91 (MQ=255)
                                                           aCTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTtttatt                                 >  1:937125/1‑91 (MQ=255)
                                                             tataCAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTtttattta                               >  1:1918056/1‑91 (MQ=255)
                                                             tataCAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTtttattta                               <  2:2998162/91‑1 (MQ=255)
                                                              ataCAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTttatttat                              <  2:1604849/91‑1 (MQ=255)
                                                               taCAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATc                             >  2:2309737/1‑91 (MQ=255)
                                                               taCAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATc                             >  2:3895805/1‑91 (MQ=255)
                                                               taCAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATc                             >  2:542473/1‑91 (MQ=255)
                                                                aCAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCg                            <  1:1597723/91‑1 (MQ=255)
                                                                aCAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCg                            >  1:633180/1‑91 (MQ=255)
                                                                 cAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGc                           >  1:323277/1‑91 (MQ=255)
                                                                 cAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGc                           >  2:620434/1‑91 (MQ=255)
                                                                  aaCTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCa                          <  1:2160899/91‑1 (MQ=255)
                                                                   aCTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCAt                         >  1:506557/1‑91 (MQ=255)
                                                                   aCTCCTAGTAAAATATGTAAAAAGTAATGTTTACGACTCAAAATTACACAAAATAAAGAAATATAACTCCATCAAGTTTTATTTATCACAt                         >  2:3530546/1‑91 (MQ=255)
                                                                    cTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCAtt                        <  2:2083442/91‑1 (MQ=255)
                                                                      ccTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTc                      >  1:3041518/1‑91 (MQ=255)
                                                                        tAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCtt                    <  1:2984951/91‑1 (MQ=255)
                                                                           tAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGaa                 <  2:2516372/91‑1 (MQ=255)
                                                                           tAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGaa                 >  2:1779224/1‑91 (MQ=255)
                                                                             aCATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAaga               <  1:2266261/91‑1 (MQ=255)
                                                                             aCATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAaga               <  2:1369875/91‑1 (MQ=255)
                                                                              cATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAagag              >  2:2179566/1‑91 (MQ=255)
                                                                               atatGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAagaga             <  1:2810391/91‑1 (MQ=255)
                                                                               atatGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAagaga             <  2:422439/91‑1 (MQ=255)
                                                                                tatGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAagagag            >  1:1689518/1‑91 (MQ=255)
                                                                                tatGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAagagag            >  2:1192980/1‑91 (MQ=255)
                                                                                 atGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAAGAGAGc           >  2:1310533/1‑91 (MQ=255)
                                                                                  tGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAAGAGAGCa          >  1:3410857/1‑91 (MQ=255)
                                                                                   gTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAAGAGAGCAt         >  2:1807868/1‑91 (MQ=255)
                                                                                    tAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAAGAGAGCAtt        <  1:3531901/91‑1 (MQ=255)
                                                                                      aaCAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAAGAGAGCATTAt      <  1:628167/91‑1 (MQ=255)
                                                                                      aaCAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAAGAGAGCATTAt      <  2:906148/91‑1 (MQ=255)
                                                                                        cAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAAGAGAGCATTATCa    >  2:2861392/1‑91 (MQ=255)
                                                                                        cAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAAGAGAGCATTATCa    >  2:2086834/1‑91 (MQ=255)
                                                                                          gTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAAGAGCGCATTATCAtt  >  2:831250/1‑91 (MQ=255)
                                                                                          |                                                                                          
TAGTAATTGACATTTTTCAATTCAAAGACTGGTGTCATCGTATTCTCACCTCGCATTCAACTATACAACTCCTAGTAACATATGTAAACAGTAATGTTTACGACTCAAAATTAGACAAAATAAAGAGATATGCCCCCTTCAAGTTTTATTTATCGCATTTCTTGAAGAGAGCATTATCATT  >  USA300TCH1516_ALE/1667434‑1667614

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