Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 770,883 T→C 9.8% intergenic (+69/‑181) 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,8830TC 9.8% 681.8 / 43.2 232intergenic (+69/‑181)ybaK/glcRCys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK/HTH‑type transcriptional repressor GlcR
Reads supporting (aligned to +/- strand):  ref base T (97/111);  new base C (13/10);  total (111/121)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 244 total aligned reads displayed.

TGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTAT  >  USA300TCH1516_ALE/770793‑770973
                                                                                          |                                                                                          
tGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAtt                                                                                            <  1:2981300/91‑1 (MQ=255)
tGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAtt                                                                                            <  2:1862805/91‑1 (MQ=255)
 ggAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAttt                                                                                           >  1:813227/1‑91 (MQ=255)
  gAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTc                                                                                          >  2:2463069/1‑91 (MQ=255)
    gCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTg                                                                                        >  1:2645217/1‑91 (MQ=255)
    gCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTg                                                                                        >  1:1354957/1‑91 (MQ=255)
      aGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGcc                                                                                      <  2:2740997/91‑1 (MQ=255)
       gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTgtca                                                                                     <  1:3681764/91‑4 (MQ=17)
       gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCca                                                                                     <  2:398888/91‑1 (MQ=255)
       gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCca                                                                                     <  1:1398619/91‑1 (MQ=255)
        ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcac                                                                                    >  2:1600255/1‑91 (MQ=255)
        ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcac                                                                                    <  2:3872525/91‑1 (MQ=255)
         tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTcaca                                                                                   <  1:3571723/91‑1 (MQ=14)
         tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcaca                                                                                   >  1:621358/1‑91 (MQ=255)
           tCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTTGGTGTGTGAAAGATATGGCGTTCTCAGAGAACTCCTCTCTGCCACACt                                                                                 <  2:3676653/91‑1 (MQ=255)
           tCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACt                                                                                 <  1:3113137/91‑1 (MQ=255)
            ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTc                                                                                <  1:1824863/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)
                    gaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttga                                                                        <  1:614109/91‑1 (MQ=255)
                    gaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttga                                                                        <  1:1065641/91‑1 (MQ=255)
                     aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgat                                                                       <  2:62693/91‑1 (MQ=255)
                       taataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattg                                                                     >  2:19544/1‑91 (MQ=255)
                       taataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattg                                                                     >  1:3186472/1‑91 (MQ=255)
                       taataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattg                                                                     >  2:2668342/1‑91 (MQ=255)
                        aataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattga                                                                    >  1:2443889/1‑91 (MQ=255)
                         ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATATTATTCATTTCTGCCACACTCCTTTttgattgta                                                                   >  2:1437162/1‑89 (MQ=255)
                         ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGaa                                                                   <  2:2877971/91‑1 (MQ=255)
                           aacaacAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAAtt                                                                 <  2:3744699/91‑1 (MQ=255)
                             caacaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAg                                                               <  1:3312877/91‑1 (MQ=255)
                               acaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCa                                                             <  1:1756670/91‑1 (MQ=255)
                                 aaaaaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAtt                                                           <  2:1533432/91‑1 (MQ=255)
                                  aaaaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAttt                                                          <  1:2850817/91‑1 (MQ=255)
                                     aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                                       <  2:1769753/91‑1 (MQ=255)
                                     aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                                       <  1:1071285/91‑1 (MQ=255)
                                       tATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGa                                                     >  1:2762081/1‑91 (MQ=255)
                                        aTGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGAt                                                    <  1:2236100/91‑1 (MQ=255)
                                         tGGGGCAATATTAGGCGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATc                                                   <  2:3371064/91‑1 (MQ=255)
                                         tGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACCCCTTTTTGATTGAATTAGCATTTTACGATc                                                   <  2:1559250/91‑1 (MQ=255)
                                          ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                                  >  1:2077274/1‑91 (MQ=255)
                                          ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                                  >  2:1668204/1‑91 (MQ=255)
                                          ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                                  >  1:951838/1‑91 (MQ=255)
                                          ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                                  >  2:1047619/1‑91 (MQ=255)
                                           gggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCAt                                                 <  2:1393197/91‑1 (MQ=255)
                                            ggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATa                                                <  2:2148345/91‑1 (MQ=255)
                                            ggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATa                                                <  2:478048/91‑1 (MQ=255)
                                            ggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATa                                                <  1:1104513/91‑1 (MQ=255)
                                              cAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaaa                                              <  2:1269881/91‑1 (MQ=255)
                                              cAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaaa                                              >  1:1693782/1‑91 (MQ=255)
                                              cAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaaa                                              <  2:3474697/91‑1 (MQ=255)
                                                aGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACa                                            >  2:2471623/1‑91 (MQ=255)
                                                   ttAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTc                                         <  2:1441887/91‑1 (MQ=255)
                                                   ttAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTc                                         >  2:3589922/1‑91 (MQ=255)
                                                   ttAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTc                                         <  2:3711137/91‑1 (MQ=255)
                                                   ttAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGAACATAAACAGTc                                         >  1:3510308/1‑91 (MQ=255)
                                                     aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                       >  1:1520523/1‑91 (MQ=255)
                                                     aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                       >  2:677235/1‑91 (MQ=255)
                                                     aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                       >  1:2365357/1‑91 (MQ=255)
                                                     aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                       <  2:2410400/91‑1 (MQ=255)
                                                     aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                       >  2:409189/1‑91 (MQ=255)
                                                       gAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATta                                     <  2:856237/91‑1 (MQ=255)
                                                       gAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATta                                     >  1:1454316/1‑91 (MQ=255)
                                                        aGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGTATTAGCATTTTACGATCATAAACAGTCATtat                                    >  2:2293872/1‑91 (MQ=14)
                                                           gtgACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt                                 <  2:3135943/91‑1 (MQ=255)
                                                             ggcAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTg                               >  2:2321562/3‑91 (MQ=17)
                                                               cAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                             <  1:861489/91‑1 (MQ=255)
                                                               cAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                             <  1:1612294/91‑1 (MQ=255)
                                                               cAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                             >  2:604640/1‑91 (MQ=18)
                                                               cAGAAATGAATTTTATAAAAATTTCATCTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                             >  2:3037331/1‑91 (MQ=21)
                                                               cAGAAATGAATTTTATAAAAATTTCATCTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                             >  2:3037164/1‑91 (MQ=21)
                                                                agaaaTGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGt                            >  2:3535043/5‑91 (MQ=21)
                                                                   gaTGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGtttt                         >  2:2046394/1‑88 (MQ=255)
                                                                   gaTGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAtt                         >  1:1182461/1‑91 (MQ=255)
                                                                   gaTGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAtt                         <  1:2885999/91‑1 (MQ=255)
                                                                      gAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGa                      >  2:208215/1‑91 (MQ=255)
                                                                       aaaTTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGaa                     <  2:2862882/91‑1 (MQ=255)
                                                                       aaaTTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGaa                     >  2:80216/1‑91 (MQ=255)
                                                                       aaaTTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGaa                     >  1:2557675/1‑91 (MQ=255)
                                                                       aaaTTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGaa                     <  2:2863315/91‑1 (MQ=255)
                                                                         aTTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACa                   <  2:2811971/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:3627948/91‑1 (MQ=255)
                                                                           ttttATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATa                 >  2:1079356/1‑91 (MQ=255)
                                                                             ttATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATaaa               <  1:346043/91‑1 (MQ=255)
                                                                             ttATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATaaa               <  1:2475903/91‑1 (MQ=255)
                                                                                 aaaaTTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATg           <  2:1119650/91‑1 (MQ=255)
                                                                                    aTTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTaa        >  1:1274769/1‑91 (MQ=255)
                                                                                     ttCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAAt       <  2:950016/91‑1 (MQ=255)
                                                                                      tCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAAtt      >  2:1703579/1‑91 (MQ=255)
                                                                                      tCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAAtt      >  1:2921388/1‑91 (MQ=255)
                                                                                       cATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAAttt     >  1:507382/1‑91 (MQ=255)
                                                                                        aTTTCTGCCACACTCCTTTTTTTTTGTTTTTTCATTTTACGATCATAAACATTCATTATAATTGAGTATTTGAACATAAAAATGTAAtttt    >  2:278875/1‑91 (MQ=255)
                                                                                          ttCTGCCACACTCCTTTTTGCTTGAATTAGCATTTTACGATCATAAACAGTCATAATAATTGAGTATTTGAACATAAAAATGTAATTTTac  >  2:3835593/1‑90 (MQ=255)
                                                                                          ttCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTAt  >  1:164572/1‑91 (MQ=255)
                                                                                          |                                                                                          
TGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTAT  >  USA300TCH1516_ALE/770793‑770973

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: