Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 770,888 C→T 6.3% intergenic (+74/‑176) ybaK → / → glcR Cys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK/HTH‑type transcriptional repressor GlcR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE770,8880CT 6.3% 374.5 / 10.0 111intergenic (+74/‑176)ybaK/glcRCys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK/HTH‑type transcriptional repressor GlcR
Reads supporting (aligned to +/- strand):  ref base C (52/52);  new base T (4/3);  total (56/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 112 total aligned reads displayed.

CAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCGTC  >  USA300TCH1516_ALE/770798‑770977
                                                                                            |                                                                                         
cAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGc                                                                                           >  1:2112570/1‑91 (MQ=255)
  gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCca                                                                                         <  1:934096/91‑1 (MQ=255)
   ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCcac                                                                                        <  1:378467/91‑1 (MQ=255)
    tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTcaca                                                                                       <  1:1753057/91‑1 (MQ=14)
     aTCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCacac                                                                                      >  1:1133654/1‑91 (MQ=255)
       ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTc                                                                                    <  1:3239499/91‑1 (MQ=255)
        cATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTCACACTcc                                                                                   <  1:497038/91‑1 (MQ=14)
           gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATGAAA‑‑TTCATTTCTGCCACACTCCttt                                                                                >  1:4052960/1‑91 (MQ=255)
           gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCttt                                                                                <  1:4027861/91‑1 (MQ=255)
           gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCttt                                                                                <  1:2326469/91‑1 (MQ=255)
            aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCtttt                                                                               >  1:2332863/1‑91 (MQ=255)
            aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTCACACTCCtaat                                                                               <  1:2582593/91‑4 (MQ=2)
            aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCtttt                                                                               <  1:3189800/91‑1 (MQ=255)
            aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCtttt                                                                               >  1:1852143/1‑91 (MQ=255)
            aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCtttt                                                                               >  1:1852017/1‑91 (MQ=255)
            aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCtttt                                                                               >  1:243119/1‑91 (MQ=255)
               gaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttga                                                                            >  1:3818351/1‑91 (MQ=255)
               gaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttga                                                                            <  1:1640639/91‑1 (MQ=255)
                aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttgat                                                                           >  1:2520084/1‑91 (MQ=255)
                aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttgat                                                                           <  1:1399048/91‑1 (MQ=255)
                aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttgat                                                                           >  1:3571537/1‑91 (MQ=255)
                   aataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttgattga                                                                        >  1:1018175/1‑91 (MQ=255)
                    ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGaa                                                                       >  1:1970447/1‑91 (MQ=255)
                      aacaacAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAAtt                                                                     <  1:337186/91‑1 (MQ=255)
                      aacaacAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAAtt                                                                     <  1:3112305/91‑1 (MQ=255)
                      aacaacAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAAtt                                                                     <  1:336878/91‑1 (MQ=255)
                       acaacaAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTa                                                                    >  1:497326/1‑91 (MQ=255)
                       acaacaAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTa                                                                    >  1:496747/1‑91 (MQ=255)
                        caacaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATCTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAg                                                                   >  1:3286612/1‑91 (MQ=255)
                           caaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAt                                                                >  1:1737287/1‑91 (MQ=255)
                            aaaaaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAtt                                                               <  1:4020534/91‑1 (MQ=255)
                               aaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTa                                                            >  1:3789443/1‑91 (MQ=255)
                                aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                                           <  1:3160809/91‑1 (MQ=255)
                                aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                                           <  1:2245523/91‑1 (MQ=255)
                                aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                                           <  1:2246148/91‑1 (MQ=255)
                                   aTGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGAt                                                        >  1:2104325/1‑91 (MQ=255)
                                    tGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATc                                                       <  1:2182310/91‑1 (MQ=255)
                                    tGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATc                                                       >  1:1571222/1‑91 (MQ=255)
                                     ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                                      >  1:510250/1‑91 (MQ=255)
                                     ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                                      >  1:2336528/1‑91 (MQ=255)
                                     ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                                      >  1:2593080/1‑91 (MQ=255)
                                     ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                                      >  1:510064/1‑91 (MQ=255)
                                      gggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCAt                                                     <  1:2166780/91‑1 (MQ=255)
                                      gggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCAt                                                     <  1:2166830/91‑1 (MQ=255)
                                       ggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATa                                                    >  1:513040/1‑91 (MQ=255)
                                              ttACGAGTGCGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTc                                             <  1:3727446/91‑1 (MQ=255)
                                                aGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                           <  1:1999169/91‑1 (MQ=255)
                                                aGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                           >  1:2663868/1‑91 (MQ=255)
                                                aGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                           <  1:3772122/91‑1 (MQ=255)
                                                 ggAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAtt                                          >  1:1452989/1‑91 (MQ=255)
                                                 ggAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAtt                                          <  1:911528/91‑1 (MQ=255)
                                                 ggAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAtt                                          >  1:1139591/1‑91 (MQ=255)
                                                 ggAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAtt                                          >  1:575465/1‑91 (MQ=255)
                                                  gAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATta                                         <  1:1178726/91‑1 (MQ=255)
                                                  gAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATta                                         >  1:322715/1‑91 (MQ=255)
                                                   aGTGTGGCAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtat                                        >  1:1802078/1‑91 (MQ=14)
                                                   aGTGTGGCAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtat                                        >  1:83481/1‑91 (MQ=14)
                                                    gtgtGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtata                                       >  1:843355/1‑91 (MQ=255)
                                                      gtgACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt                                     <  1:3953783/91‑1 (MQ=255)
                                                      gtgACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt                                     <  1:4020385/91‑1 (MQ=255)
                                                      gtgACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt                                     <  1:4020387/91‑1 (MQ=255)
                                                       tgACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAtt                                    <  1:3198840/91‑1 (MQ=255)
                                                        gACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTg                                   <  1:2960081/91‑1 (MQ=255)
                                                        gACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTg                                   <  1:1381512/91‑1 (MQ=255)
                                                        gACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTg                                   <  1:1381638/91‑1 (MQ=255)
                                                         aCAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGa                                  <  1:4112789/91‑1 (MQ=255)
                                                          cAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                                 >  1:10473/1‑91 (MQ=255)
                                                          cAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                                 >  1:3009693/1‑91 (MQ=18)
                                                           agagATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGt                                >  1:3580204/1‑91 (MQ=255)
                                                              gaTGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAtt                             >  1:140868/1‑91 (MQ=255)
                                                              gaTGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAtt                             >  1:2184089/1‑91 (MQ=255)
                                                              gaTGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAtt                             >  1:2607177/1‑91 (MQ=255)
                                                                 gAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGa                          >  1:2156992/1‑91 (MQ=255)
                                                                  tgaattttATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTg                           >  1:100679/5‑91 (MQ=21)
                                                                      ttttATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATa                     <  1:3363275/91‑1 (MQ=255)
                                                                      ttttATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATa                     >  1:728943/1‑91 (MQ=255)
                                                                      ttttATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATa                     <  1:2246211/91‑1 (MQ=255)
                                                                         tatataA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATgaaa                  >  1:3271297/1‑91 (MQ=255)
                                                                         tataAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATaaaa                  >  1:3946832/1‑91 (MQ=255)
                                                                         tataAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATaaaa                  >  1:3946893/1‑91 (MQ=255)
                                                                           taAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAAt                >  1:1562805/1‑91 (MQ=255)
                                                                             aaa‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGt              >  1:3057596/1‑91 (MQ=255)
                                                                                    cATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAAttt         <  1:125838/91‑1 (MQ=255)
                                                                                    cATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAAttt         >  1:3247955/1‑91 (MQ=255)
                                                                                      tttCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTa       >  1:2447329/1‑91 (MQ=255)
                                                                                      tttCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTa       <  1:2661618/91‑1 (MQ=255)
                                                                                      tttCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTa       <  1:2660705/91‑1 (MQ=255)
                                                                                       ttCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGGGTATTTGAACATAAAAATGTAATTTTAt      >  1:1786569/1‑91 (MQ=255)
                                                                                        tCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATc     >  1:277024/1‑91 (MQ=255)
                                                                                        tCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATc     <  1:302090/91‑1 (MQ=255)
                                                                                         cTGCTACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCg    <  1:3452359/91‑1 (MQ=255)
                                                                                         cTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCg    <  1:1741211/91‑1 (MQ=255)
                                                                                         cTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCg    <  1:538577/91‑1 (MQ=255)
                                                                                         cTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCg    <  1:1003351/91‑1 (MQ=255)
                                                                                          tGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCGt   >  1:787558/1‑91 (MQ=255)
                                                                                          tGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCGt   >  1:789194/1‑91 (MQ=255)
                                                                                           gCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCGTc  <  1:3004559/91‑1 (MQ=255)
                                                                                           gCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCGTc  <  1:4076245/91‑1 (MQ=255)
                                                                                           gCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCGTc  <  1:1205335/91‑1 (MQ=255)
                                                                                           gCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCGTc  <  1:2164897/91‑1 (MQ=255)
                                                                                            |                                                                                         
CAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCGTC  >  USA300TCH1516_ALE/770798‑770977

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