Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 171,069 A→T 6.5% intergenic (+418/‑32) USA300HOU_RS07235 → / → adhE hypothetical protein/Aldehyde‑alcohol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE171,0690AT 6.5% 617.9 / 20.5 184intergenic (+418/‑32)USA300HOU_RS07235/adhEhypothetical protein/Aldehyde‑alcohol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (64/108);  new base T (8/4);  total (72/112)
Fisher's exact test for biased strand distribution p-value = 6.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 207 total aligned reads displayed.

TGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTAT  >  USA300TCH1516_ALE/170979‑171159
                                                                                          |                                                                                          
 gATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAaat                                                                                           <  2:587486/91‑1 (MQ=255)
   tatTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAaataa                                                                                         <  2:162503/91‑1 (MQ=255)
   tatTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAaataa                                                                                         <  2:3088555/91‑1 (MQ=255)
    atTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAaataat                                                                                        <  1:2568733/91‑1 (MQ=255)
      ttGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTc                                                                                      <  1:266537/91‑1 (MQ=255)
       tGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTca                                                                                     >  1:1086435/1‑91 (MQ=255)
         ttCAAATAATAGTCGAAATTAAATTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGAGATGTACTATAATGTTTGTGAAATAATTcaca                                                                                   <  2:2961344/91‑1 (MQ=255)
         ttCAAATAATAGTCGAAATTAAATTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGAGATGTACTATAATGTTTGTGAAATAATTcaca                                                                                   <  2:2959740/91‑1 (MQ=255)
          tCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACaa                                                                                  <  1:676737/91‑1 (MQ=255)
          tCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACaa                                                                                  <  1:2995359/91‑1 (MQ=255)
            aaaTAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAg                                                                                <  2:3564996/91‑1 (MQ=255)
             aataatATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTcac                                                                                    <  2:2425895/86‑1 (MQ=255)
             aataatATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGt                                                                               <  1:2522124/91‑1 (MQ=255)
              ataataTTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGta                                                                              >  1:183373/1‑91 (MQ=255)
               taataTTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAATTATTTCACaaacatt                                                                             <  1:3665618/91‑5 (MQ=39)
                aataTTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGtata                                                                            <  2:1341986/91‑1 (MQ=255)
                aataTTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGtata                                                                            <  2:788540/91‑1 (MQ=255)
                 atatTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATaa                                                                           >  2:127765/1‑91 (MQ=255)
                  tatTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATaaa                                                                          <  2:962108/91‑1 (MQ=255)
                  tatTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATaaa                                                                          <  1:1118074/91‑1 (MQ=255)
                  tatTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATAGCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATaaa                                                                          <  2:2059982/91‑1 (MQ=255)
                   atTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAg                                                                         <  2:91824/91‑1 (MQ=255)
                   atTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAg                                                                         <  2:2333075/91‑1 (MQ=255)
                    ttCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAgg                                                                        >  1:642487/1‑91 (MQ=255)
                     tCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGa                                                                       >  1:2227212/1‑91 (MQ=255)
                     tCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGa                                                                       <  2:2614366/91‑1 (MQ=255)
                     tCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGa                                                                       >  1:266983/1‑91 (MQ=255)
                      cGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAg                                                                      <  2:1039216/91‑1 (MQ=255)
                      cGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAg                                                                      <  1:3528588/91‑1 (MQ=255)
                        aaaTTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTg                                                                    <  2:199660/91‑1 (MQ=255)
                          aTTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGt                                                                  >  1:1652282/1‑91 (MQ=255)
                          aTTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGt                                                                  <  1:3377132/91‑1 (MQ=255)
                            tAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTg                                                                <  1:806509/91‑1 (MQ=255)
                             aaaCTTTTTTGTATCGGCTTTGCTTTATATCCTGAGAGAGATGTACTATAATGGTTGTGAGATAATTCACAAAGTATAACGGAGTGGTTGt                                                               <  2:2187838/91‑1 (MQ=255)
                              aaCTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGta                                                              <  1:889232/91‑1 (MQ=255)
                               aCTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGtat                                                             <  2:1436488/91‑1 (MQ=255)
                                cTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGtata                                                            >  2:617889/1‑91 (MQ=255)
                                 ttttttGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTatat                                                           >  1:919022/1‑91 (MQ=255)
                                   ttttGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGt                                                         >  1:1809968/1‑91 (MQ=255)
                                   ttttGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGt                                                         <  1:1423207/91‑1 (MQ=255)
                                   ttttGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGt                                                         <  2:3628631/91‑1 (MQ=255)
                                      tGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAATTATAAAGGAGTGGTTGTATATGTTaa                                                      <  2:2454680/91‑1 (MQ=255)
                                         ataGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACta                                                   <  2:542209/91‑1 (MQ=255)
                                          taGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACtat                                                  >  1:2022854/1‑91 (MQ=255)
                                            gAATTTTCTTTATATCCTGAGAGACATGTACTATAAAGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATAc                                                <  2:1534395/91‑1 (MQ=255)
                                                   cTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGaaaa                                         >  2:3857866/1‑91 (MQ=255)
                                                   cTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGaaaa                                         >  2:505968/1‑91 (MQ=255)
                                                    tttATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGaaaaa                                        <  1:2547245/91‑1 (MQ=255)
                                                    tttATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGaaaaa                                        >  2:2675072/1‑91 (MQ=255)
                                                     ttATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAg                                       <  2:1086435/91‑1 (MQ=255)
                                                     ttATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAg                                       <  2:2763035/91‑1 (MQ=255)
                                                     ttATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAg                                       <  1:3619371/91‑1 (MQ=255)
                                                      tataTCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGa                                      >  1:1373219/1‑91 (MQ=255)
                                                       atatCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGaa                                     <  2:1063130/91‑1 (MQ=255)
                                                        tatCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGaaa                                    >  1:1706404/1‑91 (MQ=255)
                                                         atCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGaaaa                                   <  2:1853171/91‑1 (MQ=255)
                                                         atCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACATGAAAAAGaaaa                                   <  2:2292332/91‑1 (MQ=255)
                                                           ccTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATc                                 >  1:3234828/1‑91 (MQ=255)
                                                             tGAGAGACATGTACTATAATGTTTGTGGAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATGGt                               <  2:2080818/91‑1 (MQ=255)
                                                              gagagaCATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTg                              <  1:2367194/91‑1 (MQ=255)
                                                              gagagaCATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTg                              <  2:806958/91‑1 (MQ=255)
                                                               agagaCATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTgg                             >  2:573326/1‑91 (MQ=255)
                                                                gagaCATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGa                            <  2:3328000/91‑1 (MQ=255)
                                                                 agaCATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGAt                           <  2:2838278/91‑1 (MQ=255)
                                                                  gaCATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATc                          >  1:452441/1‑91 (MQ=255)
                                                                   aCATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCg                         <  2:1676374/91‑1 (MQ=255)
                                                                    cATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGa                        >  1:3420065/1‑91 (MQ=255)
                                                                    cATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGa                        <  1:2446522/91‑1 (MQ=255)
                                                                    cATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGCTAACTATACCTGAAAAAGAAAATCGTGGATCGa                        >  1:3181180/1‑91 (MQ=255)
                                                                    aaTGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGa                        <  2:2527033/90‑1 (MQ=255)
                                                                     aTGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGaa                       >  2:1264427/1‑91 (MQ=255)
                                                                      tGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGaaa                      >  1:1785576/1‑91 (MQ=255)
                                                                      tGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGaaa                      <  1:617889/91‑1 (MQ=255)
                                                                         aCTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGaa                   <  1:2504109/91‑1 (MQ=255)
                                                                          cTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAAc                  <  1:1286518/91‑1 (MQ=255)
                                                                          cTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAAc                  <  1:940376/91‑1 (MQ=255)
                                                                          cTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAAc                  <  2:1360587/91‑1 (MQ=255)
                                                                            ataATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACaa                <  2:1016835/91‑1 (MQ=255)
                                                                             taATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACaag               <  1:2876814/91‑1 (MQ=255)
                                                                              aaTGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAGAAGAAAATCGTGGATCGAAAGAACaaga              <  1:1667249/91‑1 (MQ=255)
                                                                              aaTGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACaaga              >  2:1207783/1‑91 (MQ=255)
                                                                              aaTGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACaaga              >  1:2681868/1‑91 (MQ=255)
                                                                               aTGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACaagaa             <  1:1028134/91‑1 (MQ=255)
                                                                                  tttGTGAATTATTTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTg          >  2:1142692/1‑91 (MQ=255)
                                                                                  tttGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTg          <  1:1985693/91‑1 (MQ=255)
                                                                                   ttGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTgg         >  1:260906/1‑91 (MQ=255)
                                                                                    tgtgAATTATTTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGc        >  2:2265536/1‑91 (MQ=255)
                                                                                    tgtgAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGc        <  2:1861868/91‑1 (MQ=255)
                                                                                    tgtgAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGc        <  1:2172122/91‑1 (MQ=255)
                                                                                    tgtgAAACAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACTAGAAGTGGc        <  1:1950516/91‑1 (MQ=255)
                                                                                     gtgAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCa       >  1:2774447/1‑91 (MQ=255)
                                                                                     gtgAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCa       >  1:2189421/1‑91 (MQ=255)
                                                                                      tgAATTATTTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCaa      >  1:2179439/1‑91 (MQ=255)
                                                                                      tgAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCaa      >  1:1140398/1‑91 (MQ=255)
                                                                                      tgAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCaa      >  1:2383343/1‑91 (MQ=255)
                                                                                       gAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAAt     <  2:3265387/91‑1 (MQ=255)
                                                                                       gAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAAt     <  1:1526523/91‑1 (MQ=255)
                                                                                        aaaTAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAAtt    <  1:833420/91‑1 (MQ=255)
                                                                                         aataatTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTa   >  1:722652/1‑91 (MQ=255)
                                                                                          ataatTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTAt  <  2:2388799/91‑1 (MQ=255)
                                                                                          |                                                                                          
TGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTAT  >  USA300TCH1516_ALE/170979‑171159

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