Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 770,855 A→G 6.1% intergenic (+41/‑209) 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,8550AG 6.1% 737.9 / 24.9 231intergenic (+41/‑209)ybaK/glcRCys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK/HTH‑type transcriptional repressor GlcR
Reads supporting (aligned to +/- strand):  ref base A (105/112);  new base G (3/11);  total (108/123)
Fisher's exact test for biased strand distribution p-value = 5.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

Only 100 of 233 total aligned reads displayed.

ATTTGATTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTG  >  USA300TCH1516_ALE/770765‑770944
                                                                                          |                                                                                         
   tGATTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACag                                                                                        >  2:3724722/1‑91 (MQ=255)
    gATTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAGAAGTATGGGGCAAGATTAGGAGTGTGACaga                                                                                       <  2:2957488/91‑1 (MQ=255)
    gATTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACaga                                                                                       <  2:3585981/91‑1 (MQ=255)
      ttACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAgaga                                                                                     <  1:1029018/91‑1 (MQ=255)
      ttACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAgaga                                                                                     <  2:676091/91‑1 (MQ=255)
      ttACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAgaga                                                                                     >  2:3767639/1‑91 (MQ=255)
       tACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGAt                                                                                    >  2:862782/1‑91 (MQ=255)
         cAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGa                                                                                  <  1:1894753/91‑1 (MQ=255)
          aaTAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAATAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGaa                                                                                 <  2:3776862/91‑1 (MQ=255)
           aTAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGaaa                                                                                >  2:1403475/1‑91 (MQ=255)
           aTAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGaaa                                                                                >  1:1675244/1‑91 (MQ=255)
             aaCTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAAAAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAAtt                                                                              <  2:1652890/91‑1 (MQ=255)
                  aaGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTtat                                                                         >  1:3150480/1‑91 (MQ=255)
                  aaGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTtat                                                                         <  1:3724722/91‑1 (MQ=255)
                  aaGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTtat                                                                         >  2:1737600/1‑91 (MQ=255)
                   aGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTtata                                                                        >  2:2647759/1‑91 (MQ=255)
                   aGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTtata                                                                        <  1:3537067/91‑1 (MQ=255)
                   aGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTtata                                                                        >  2:1532647/1‑91 (MQ=255)
                      cAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATaaaa                                                                     >  2:1702686/1‑91 (MQ=255)
                      cAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATaaaa                                                                     <  1:1176530/91‑1 (MQ=255)
                       aaaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAt                                                                    <  2:2630674/91‑1 (MQ=255)
                       aaaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAt                                                                    >  2:1315886/1‑91 (MQ=255)
                       aaaaTTGGAGCAGTTAGCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAt                                                                    <  2:647672/91‑1 (MQ=255)
                        aaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAtt                                                                   <  2:2918293/91‑1 (MQ=255)
                          aTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCa                                                                 <  2:2992591/91‑1 (MQ=255)
                            tGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAtt                                                               <  1:2981300/91‑1 (MQ=255)
                             ggAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAttt                                                              >  2:804837/1‑91 (MQ=255)
                              gAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTc                                                             >  2:2463121/1‑91 (MQ=255)
                               ggCAGTTAGCCATGAATGAATAAAAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATctct                                                            <  2:2026041/90‑1 (MQ=12)
                               aGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCt                                                            >  1:1732009/1‑91 (MQ=255)
                                gCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTg                                                           >  1:2645217/1‑91 (MQ=255)
                                gCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTg                                                           >  1:1354957/1‑91 (MQ=255)
                                 cAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTgt                                                          <  2:2786196/91‑2 (MQ=17)
                                  aGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGcc                                                         <  2:2740997/91‑1 (MQ=255)
                                   gTTATCCATGAATGATTAATAACACCGAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCca                                                        <  2:398889/91‑1 (MQ=255)
                                   gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTgtca                                                        <  1:2720573/91‑4 (MQ=17)
                                   gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCca                                                        <  1:1398619/91‑1 (MQ=255)
                                    ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcac                                                       >  2:1612294/1‑91 (MQ=255)
                                    tgaTCCATGAATGATAAAAAAAAAAAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcac                                                       <  2:3872442/89‑1 (MQ=255)
                                     tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTcaca                                                      <  2:97937/91‑1 (MQ=14)
                                     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)
                                      aTCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCacac                                                     <  1:2306590/91‑1 (MQ=255)
                                       tCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACt                                                    <  1:1069107/91‑1 (MQ=14)
                                       tCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAGTTCTGCCACACt                                                    <  1:2991387/91‑1 (MQ=255)
                                        ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTATGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTc                                                   <  2:1425979/91‑1 (MQ=255)
                                        ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTc                                                   <  1:1824863/91‑1 (MQ=255)
                                        ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTc                                                   <  2:614850/91‑1 (MQ=255)
                                         cATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTcc                                                  <  2:1136294/91‑1 (MQ=255)
                                            gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGATAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCttt                                               >  1:3868951/1‑91 (MQ=255)
                                            gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCttt                                               >  2:2363440/1‑91 (MQ=255)
                                             aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                              >  1:1987714/1‑91 (MQ=255)
                                             aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                              >  1:3803276/1‑91 (MQ=255)
                                              atgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcac                                                       <  2:2310612/81‑1 (MQ=255)
                                                gaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttga                                           <  1:614109/91‑1 (MQ=255)
                                                 aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgat                                          <  1:2649648/91‑1 (MQ=255)
                                                  ttAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTATGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgatt                                         >  1:3275594/1‑91 (MQ=255)
                                                   taataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattg                                        >  2:19544/1‑91 (MQ=255)
                                                   taataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttg                                            >  1:2355851/1‑87 (MQ=255)
                                                    aataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattga                                       >  1:2443889/1‑91 (MQ=255)
                                                      taaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAAtttt                                                                            >  1:707362/1‑52 (MQ=255)
                                                      taaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAAtttt                                                                            <  2:707362/52‑1 (MQ=255)
                                                       aacaacAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAAtt                                    >  2:2236100/1‑91 (MQ=255)
                                                         caacaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAg                                  <  1:3312877/91‑1 (MQ=255)
                                                          aacaaAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGc                                 <  2:2826343/91‑1 (MQ=255)
                                                           acaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCa                                <  2:409326/91‑1 (MQ=255)
                                                           acaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCa                                >  2:210698/1‑91 (MQ=255)
                                                              aaaaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAttt                             <  1:2850817/91‑1 (MQ=255)
                                                              aaaaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAttt                             >  1:107522/1‑91 (MQ=255)
                                                                 aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                          <  2:1769753/91‑1 (MQ=255)
                                                                 aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                          <  1:2081870/91‑1 (MQ=255)
                                                                 aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                          >  2:1534218/1‑91 (MQ=255)
                                                                  gTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACg                         >  2:3488803/1‑91 (MQ=255)
                                                                  gTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACg                         <  2:1591229/91‑1 (MQ=255)
                                                                   tATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGa                        <  1:970225/91‑1 (MQ=255)
                                                                    aTGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGAt                       <  1:2236100/91‑1 (MQ=255)
                                                                     tGGGGCAATATTAGGCGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATc                      <  2:3371064/91‑1 (MQ=255)
                                                                     tGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATc                      >  1:1254924/1‑91 (MQ=255)
                                                                      ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                     >  1:951401/1‑91 (MQ=255)
                                                                        ggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATa                   <  2:2148345/91‑1 (MQ=255)
                                                                          cAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaaa                 <  1:1118178/91‑1 (MQ=255)
                                                                          cAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaaa                 <  2:51549/91‑1 (MQ=255)
                                                                               ttAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTc            <  2:1441887/91‑1 (MQ=255)
                                                                               ttAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTc            <  2:3711137/91‑1 (MQ=255)
                                                                               ttAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGAACATAAACAGTc            >  1:3510308/1‑91 (MQ=255)
                                                                                 aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTTTTTTTATTAGCATTTTACGATCATAAACAGTCAt          >  2:909854/1‑91 (MQ=255)
                                                                                 aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaa                  >  2:1359074/1‑83 (MQ=255)
                                                                                 aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt          <  2:1260239/91‑1 (MQ=255)
                                                                                 aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt          >  1:1520523/1‑91 (MQ=255)
                                                                                 aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt          >  1:2365373/1‑91 (MQ=255)
                                                                                 aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt          <  2:2410400/91‑1 (MQ=255)
                                                                                  ggAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAtt         <  1:862782/91‑1 (MQ=255)
                                                                                   gAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATta        >  1:1454316/1‑91 (MQ=255)
                                                                                    aGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtat       <  2:2366989/91‑1 (MQ=255)
                                                                                    aGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtat       >  2:3875848/1‑91 (MQ=255)
                                                                                       gtgACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt    <  2:3402453/91‑1 (MQ=255)
                                                                                       gtgACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt    <  2:3135943/91‑1 (MQ=255)
                                                                                         gACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTg  <  2:369003/91‑1 (MQ=255)
                                                                                          |                                                                                         
ATTTGATTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTG  >  USA300TCH1516_ALE/770765‑770944

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: