Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 770,870 T→A 9.3% intergenic (+56/‑194) 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,8700TA 9.3% 761.0 / 50.3 237intergenic (+56/‑194)ybaK/glcRCys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK/HTH‑type transcriptional repressor GlcR
Reads supporting (aligned to +/- strand):  ref base T (99/116);  new base A (12/10);  total (111/126)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 252 total aligned reads displayed.

TAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAA  >  USA300TCH1516_ALE/770781‑770960
                                                                                         |                                                                                            
tAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAAttttt                                                                                             >  1:2682668/1‑91 (MQ=255)
  aaGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTtat                                                                                           <  2:96433/91‑1 (MQ=255)
  aaGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTtat                                                                                           >  1:3150480/1‑91 (MQ=255)
   aGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTtata                                                                                          >  2:2647759/1‑91 (MQ=255)
   aGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTtata                                                                                          <  1:3537067/91‑1 (MQ=255)
      cAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATaaaa                                                                                       <  2:3242140/91‑1 (MQ=255)
      cAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATaaaa                                                                                       >  2:1702686/1‑91 (MQ=255)
      cAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATaaaa                                                                                       <  1:1176436/91‑1 (MQ=255)
       aaaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑t                                                                                    <  2:2630674/91‑1 (MQ=255)
       aaaaTTGGAGCAGATATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑t                                                                                    >  2:1316108/1‑91 (MQ=255)
        aaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑tt                                                                                   >  1:2081103/1‑91 (MQ=255)
        aaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑tt                                                                                   >  1:2222620/1‑91 (MQ=255)
          aTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCa                                                                                 <  1:200088/91‑1 (MQ=255)
            tGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCAtt                                                                               <  2:1862805/91‑1 (MQ=255)
            tGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCAtt                                                                               <  1:2981300/91‑1 (MQ=255)
             ggAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCAttt                                                                              >  1:813227/1‑91 (MQ=255)
             ggAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCAttt                                                                              >  2:804837/1‑91 (MQ=255)
              gAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAGAA‑‑TTCATTTc                                                                             <  2:1054507/91‑1 (MQ=255)
              gAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTc                                                                             >  2:2463069/1‑91 (MQ=255)
                 cAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTgt                                                                          <  2:2786196/91‑2 (MQ=17)
                  aGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGcc                                                                         <  2:2740997/91‑1 (MQ=255)
                    ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCcac                                                                       >  2:1600255/1‑91 (MQ=255)
                    ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCcac                                                                       <  2:3872525/91‑1 (MQ=255)
                     tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTcaca                                                                      <  1:3571723/91‑1 (MQ=14)
                     tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCcaca                                                                      >  1:621358/1‑91 (MQ=255)
                     tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCcaca                                                                      <  2:3590599/91‑1 (MQ=255)
                      aTCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCacac                                                                     <  1:3398957/91‑1 (MQ=255)
                       tCCATGAATGATTTATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACt                                                                    >  2:2132856/1‑91 (MQ=255)
                       tCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACt                                                                    <  2:3746455/91‑1 (MQ=255)
                        ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTc                                                                   <  1:1824863/91‑1 (MQ=255)
                            gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGATAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCttt                                                               >  1:3868951/1‑91 (MQ=255)
                            gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCttt                                                               >  2:2363890/1‑91 (MQ=255)
                             aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCtttt                                                              >  1:3803276/1‑91 (MQ=255)
                              atgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCcac                                                                       >  1:2310612/1‑81 (MQ=255)
                              atgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCcac                                                                       <  2:2310612/81‑1 (MQ=255)
                              atgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCttttt                                                             >  2:662155/1‑91 (MQ=255)
                                gaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttga                                                           <  1:1065641/91‑1 (MQ=255)
                                gaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttga                                                           <  1:614109/91‑1 (MQ=255)
                                 aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttgat                                                          >  2:1641423/1‑91 (MQ=255)
                                 aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttgat                                                          <  1:2646717/91‑1 (MQ=255)
                                 aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttgat                                                          <  2:62693/91‑1 (MQ=255)
                                  ttAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttgatt                                                         >  2:2015534/1‑91 (MQ=255)
                                  ttAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTATGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttgatt                                                         >  1:3275594/1‑91 (MQ=255)
                                   taataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttg                                                            <  2:2355851/87‑1 (MQ=255)
                                    aataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttgattga                                                       >  1:2443889/1‑91 (MQ=255)
                                     ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATATTA‑‑TTCATTTCTGCCACACTCCTTTttgattgta                                                      >  2:1437162/1‑89 (MQ=255)
                                     ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGaa                                                      <  2:2877971/91‑1 (MQ=255)
                                     ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGaa                                                      <  1:1576720/91‑1 (MQ=255)
                                      taaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAAtttt                                                                                              <  2:707362/52‑1 (MQ=255)
                                       aacaacAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAAtt                                                    >  2:2236100/1‑91 (MQ=255)
                                       aacaacAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAAtt                                                    >  1:1295121/1‑91 (MQ=255)
                                         caacaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAg                                                  >  1:249255/1‑91 (MQ=255)
                                          aacaaAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGc                                                 <  2:2826343/91‑1 (MQ=255)
                                          aacaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGc                                                 <  1:1954382/91‑1 (MQ=255)
                                           acaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCa                                                >  2:210698/1‑91 (MQ=255)
                                           acaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCa                                                <  1:1756670/91‑1 (MQ=255)
                                             aaaaaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAtt                                              <  2:1533432/91‑1 (MQ=255)
                                              aaaaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAttt                                             <  1:2850817/91‑1 (MQ=255)
                                                 aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTTTTTTTATTAGCATTTTAc                                          >  2:3069897/1‑91 (MQ=255)
                                                 aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                          <  2:1769753/91‑1 (MQ=255)
                                                 aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                          <  1:1071285/91‑1 (MQ=255)
                                                     tGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACCCCTTTTTGATTGAATTAGCATTTTACGATc                                      <  2:1558371/91‑1 (MQ=255)
                                                      ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                     >  1:951838/1‑91 (MQ=255)
                                                       gggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCAt                                    <  2:1393197/91‑1 (MQ=255)
                                                        ggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATa                                   <  2:478048/91‑1 (MQ=255)
                                                          cAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaaa                                 <  2:3474697/91‑1 (MQ=255)
                                                          cAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaaa                                 <  1:1118178/91‑1 (MQ=255)
                                                          cAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaaa                                 >  1:1693782/1‑91 (MQ=255)
                                                          cAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaaa                                 <  2:1269881/91‑1 (MQ=255)
                                                               ttAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTc                            <  2:1441887/91‑1 (MQ=255)
                                                               ttAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTc                            <  2:3711233/91‑1 (MQ=255)
                                                                 aGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaa                                  <  1:1359074/83‑1 (MQ=255)
                                                                 aGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                          <  2:2410400/91‑1 (MQ=255)
                                                                 aGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                          >  2:409189/1‑91 (MQ=255)
                                                                 aGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                          >  1:2365357/1‑91 (MQ=255)
                                                                  ggAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAtt                         >  1:3503289/1‑91 (MQ=255)
                                                                   gAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATta                        >  1:1457310/1‑91 (MQ=255)
                                                                   gAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATta                        <  2:856237/91‑1 (MQ=255)
                                                                    aGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtat                       >  2:3875848/1‑91 (MQ=255)
                                                                     gtgtgtCAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtata                      <  2:3573350/91‑1 (MQ=255)
                                                                      tgtgACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACGCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATaa                     >  2:2991847/1‑91 (MQ=255)
                                                                       gtggcAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt                    >  1:2888936/5‑91 (MQ=14)
                                                                       gtgACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt                    <  2:3402453/91‑1 (MQ=255)
                                                                        tggcAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAtt                   >  1:946100/4‑91 (MQ=17)
                                                                         gACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTg                  <  2:369003/91‑1 (MQ=255)
                                                                           cAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                <  1:861489/91‑1 (MQ=255)
                                                                           cAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                <  2:1734138/91‑1 (MQ=255)
                                                                           cAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                >  2:604640/1‑91 (MQ=18)
                                                                           cAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                >  2:3037164/1‑91 (MQ=21)
                                                                               gaTGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGtttt            >  2:2046394/1‑88 (MQ=255)
                                                                               gaTGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAtt            >  1:1182461/1‑91 (MQ=255)
                                                                                 tGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTg          >  2:36355/1‑91 (MQ=255)
                                                                                   aaaTTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGaa        <  2:2863315/91‑1 (MQ=255)
                                                                                   aaaTTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGaa        >  1:2562273/1‑91 (MQ=255)
                                                                                     aattttATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGaa        >  1:2786620/3‑91 (MQ=25)
                                                                                     aTTTTTATAAAA‑‑TTCATTTCTGCCACACTCCGTTTTGATTGAATGAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACa      <  1:3099285/91‑1 (MQ=255)
                                                                                      tttttATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACAt     >  1:3007080/1‑91 (MQ=255)
                                                                                      tttttATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACAt     >  1:144697/1‑91 (MQ=255)
                                                                                       ttttATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATa    >  2:1079356/1‑91 (MQ=255)
                                                                                         ttATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATaaa  <  1:1641423/91‑1 (MQ=255)
                                                                                         |                                                                                            
TAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAA  >  USA300TCH1516_ALE/770781‑770960

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: