Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 770,877 A→T 9.5% intergenic (+63/‑187) 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,8770AT 9.5% 742.4 / 47.2 237intergenic (+63/‑187)ybaK/glcRCys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK/HTH‑type transcriptional repressor GlcR
Reads supporting (aligned to +/- strand):  ref base A (100/114);  new base T (12/11);  total (112/125)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

CAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAA  >  USA300TCH1516_ALE/770787‑770967
                                                                                          |                                                                                          
cAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATaaaa                                                                                            <  1:1176436/91‑1 (MQ=255)
cAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATaaaa                                                                                            <  1:1176637/91‑1 (MQ=255)
cAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATaaaa                                                                                            <  2:3242140/91‑1 (MQ=255)
 aaaaTTGGAGCAGATATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAt                                                                                           >  2:1316108/1‑91 (MQ=255)
  aaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAtt                                                                                          >  1:2222620/1‑91 (MQ=255)
  aaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAtt                                                                                          <  2:2918293/91‑1 (MQ=255)
    aTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCa                                                                                        <  2:2992591/91‑1 (MQ=255)
     ttGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAt                                                                                       >  2:1756670/1‑91 (MQ=255)
      tGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAtt                                                                                      <  1:2981299/91‑1 (MQ=255)
      tGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAtt                                                                                      <  2:1862805/91‑1 (MQ=255)
       ggAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAttt                                                                                     >  1:690489/1‑91 (MQ=255)
        gAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTc                                                                                    >  2:2463121/1‑91 (MQ=255)
         aGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCt                                                                                   >  1:1732009/1‑91 (MQ=255)
          gCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTg                                                                                  >  1:2645217/1‑91 (MQ=255)
           cAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTgt                                                                                 <  2:2786196/91‑2 (MQ=17)
             gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTgtca                                                                               <  1:3681764/91‑4 (MQ=17)
             gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTgtca                                                                               <  1:598345/91‑4 (MQ=17)
             gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCca                                                                               <  1:1398619/91‑1 (MQ=255)
              ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcac                                                                              >  2:1612294/1‑91 (MQ=255)
              ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcac                                                                              <  2:3872525/91‑1 (MQ=255)
               tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTcaca                                                                             <  1:3571723/91‑1 (MQ=14)
               tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTcaca                                                                             <  1:2057318/91‑1 (MQ=17)
               tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcaca                                                                             <  2:3590599/91‑1 (MQ=255)
               tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcaca                                                                             >  1:621358/1‑91 (MQ=255)
                 tCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACt                                                                           >  2:2711509/1‑91 (MQ=255)
                 tCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACt                                                                           <  2:3746455/91‑1 (MQ=255)
                  ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTc                                                                          <  1:1824883/91‑1 (MQ=255)
                  ccATGAATGATTAAGAACAACAAAAGGTGTGGGGCAAGATTAGGAGTGTGACAGAGATGAGATTTTTATAAAATTCATTTCTGCCACACTc                                                                          <  2:1425969/91‑1 (MQ=255)
                    atgaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCt                                                                        <  2:2117609/91‑1 (MQ=255)
                      gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGATAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCttt                                                                      >  1:3868951/1‑91 (MQ=255)
                      gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCttt                                                                      >  2:2363440/1‑91 (MQ=255)
                      gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCttt                                                                      >  1:936538/1‑91 (MQ=255)
                       aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                                                     >  1:1987714/1‑91 (MQ=255)
                        atgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcac                                                                              >  1:2310612/1‑81 (MQ=255)
                        atgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCttttt                                                                    >  2:662155/1‑91 (MQ=255)
                          gaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttga                                                                  <  1:1065641/91‑1 (MQ=255)
                           aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgat                                                                 <  2:62693/91‑1 (MQ=255)
                            ttAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTATGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgatt                                                                >  1:3275702/1‑91 (MQ=255)
                             taataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattg                                                               >  1:3186472/1‑91 (MQ=255)
                             taataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattg                                                               >  2:19544/1‑91 (MQ=255)
                              aataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattga                                                              >  1:2443889/1‑91 (MQ=255)
                               ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATATTATTCATTTCTGCCACACTCCTTTttgattgta                                                             >  2:1437162/1‑89 (MQ=255)
                               ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAGATTTCATTTCTGCCACACTCCTTTTTGATTGaa                                                             <  2:2876691/91‑1 (MQ=255)
                                           aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTTTTTTTATTAGCATTTTAc                                                 >  2:3069897/1‑91 (MQ=255)
                                           aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                                 >  2:1534218/1‑91 (MQ=255)
                                           aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                                 <  1:1071285/91‑1 (MQ=255)
                                           aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                                 <  1:2081870/91‑1 (MQ=255)
                                             tATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGa                                               >  1:2762081/1‑91 (MQ=255)
                                              aTGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGAt                                              <  2:2645217/91‑1 (MQ=255)
                                               tGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACCCCTTTTTGATTGAATTAGCATTTTACGATc                                             <  2:1559250/91‑1 (MQ=255)
                                                ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                            >  2:1047619/1‑91 (MQ=255)
                                                ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                            >  1:951838/1‑91 (MQ=255)
                                                ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                            >  2:1668204/1‑91 (MQ=255)
                                                  ggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATa                                          <  1:1104513/91‑1 (MQ=255)
                                                  ggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATa                                          <  2:478048/91‑1 (MQ=255)
                                                  ggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATa                                          <  2:2148345/91‑1 (MQ=255)
                                                    cAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaaa                                        <  2:3474697/91‑1 (MQ=255)
                                                    cAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaaa                                        <  2:1269881/91‑1 (MQ=255)
                                                    cAAGATTAGGAGTGTGACAGAGATGAAATTTTGATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaaa                                        <  1:1448426/91‑1 (MQ=255)
                                                       gATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAg                                     >  2:836283/1‑91 (MQ=255)
                                                         ttAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTc                                   <  2:3711233/91‑1 (MQ=255)
                                                         ttAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTc                                   <  2:1441887/91‑1 (MQ=255)
                                                         ttAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGAACATAAACAGTc                                   >  1:3510308/1‑91 (MQ=255)
                                                           aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTTTTTTTATTAGCATTTTACGATCATAAACAGTCAt                                 >  2:909854/1‑91 (MQ=255)
                                                           aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                 >  2:409189/1‑91 (MQ=255)
                                                           aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                 >  1:1520523/1‑91 (MQ=255)
                                                           aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                 >  1:2365357/1‑91 (MQ=255)
                                                            ggAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAtt                                >  1:2819494/1‑91 (MQ=14)
                                                              aGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGTATTAGCATTTTACGATCATAAACAGTCATtat                              >  2:2293872/1‑91 (MQ=14)
                                                              aGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtat                              <  2:1732009/91‑1 (MQ=255)
                                                              aGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtat                              >  2:3875848/1‑91 (MQ=255)
                                                               gtgtgtCAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtata                             <  2:3573350/91‑1 (MQ=255)
                                                                tgtgGCAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATaa                            >  2:794040/1‑91 (MQ=17)
                                                                tgtgACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATaa                            <  2:2456269/91‑1 (MQ=255)
                                                                 gtggcAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt                           >  1:1662013/5‑91 (MQ=14)
                                                                 gtgACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt                           <  2:3402453/91‑1 (MQ=255)
                                                                   ggcAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTg                         >  2:2321562/3‑91 (MQ=17)
                                                                   gACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTg                         <  2:369003/91‑1 (MQ=255)
                                                                     cAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                       <  1:861489/91‑1 (MQ=255)
                                                                     cAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                       <  1:1612294/91‑1 (MQ=255)
                                                                     cAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGGCATTATAATTGAg                       >  1:1757239/1‑91 (MQ=255)
                                                                     cAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                       >  2:604640/1‑91 (MQ=18)
                                                                     cAGAAATGAATTTTATAAAAATTTCATCTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                       >  2:3037199/1‑91 (MQ=21)
                                                                      agagATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGt                      <  1:1298986/91‑1 (MQ=255)
                                                                         gaTGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGtttt                   >  2:2046394/1‑88 (MQ=255)
                                                                         gaTGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAtt                   <  1:2885999/91‑1 (MQ=255)
                                                                           tGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTg                 >  2:36355/1‑91 (MQ=255)
                                                                            gAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGa                >  2:208215/1‑91 (MQ=255)
                                                                             aaaTTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGGGTATTTGaa               <  2:2862985/91‑1 (MQ=255)
                                                                             aaaTTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGaa               >  1:2557675/1‑91 (MQ=255)
                                                                               aTTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACa             <  2:2811971/91‑1 (MQ=255)
                                                                               aTTTTTATAAAATTCATTTCTGCCACACTCCGTTTTGATTGAATGAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACa             <  1:3099285/91‑1 (MQ=255)
                                                                                tttttATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACAt            <  1:2668342/91‑1 (MQ=255)
                                                                                tttttATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACAt            >  1:3007080/1‑91 (MQ=255)
                                                                                 ttttATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATa           >  2:2364596/1‑91 (MQ=255)
                                                                                 ttttATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATa           >  2:1079356/1‑91 (MQ=255)
                                                                                   ttATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATaaa         <  1:2475903/91‑1 (MQ=255)
                                                                                   ttATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATaaa         <  1:346043/91‑1 (MQ=255)
                                                                                       aaaaTTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATg     <  2:1119650/91‑1 (MQ=255)
                                                                                          aTTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTaa  >  1:1274769/1‑91 (MQ=255)
                                                                                          |                                                                                          
CAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAA  >  USA300TCH1516_ALE/770787‑770967

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: