Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 171,072 A→T 6.4% intergenic (+421/‑29) 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,0720AT 6.4% 631.2 / 22.3 188intergenic (+421/‑29)USA300HOU_RS07235/adhEhypothetical protein/Aldehyde‑alcohol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (67/109);  new base T (8/4);  total (75/113)
Fisher's exact test for biased strand distribution p-value = 6.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 210 total aligned reads displayed.

TATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTATGAT  >  USA300TCH1516_ALE/170982‑171162
                                                                                          |                                                                                          
tatTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAaataa                                                                                            >  2:1028134/1‑91 (MQ=255)
tatTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAaataa                                                                                            <  2:162503/91‑1 (MQ=255)
tatTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAaataa                                                                                            <  2:2076406/91‑1 (MQ=255)
  tttGTTCAAATAATATTCGAAATTAAACTTTTTTTTATTTAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAAtt                                                                                          >  2:1832353/1‑91 (MQ=255)
  tttGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAAtt                                                                                          <  2:1855844/91‑1 (MQ=255)
   ttGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTc                                                                                         <  1:266868/91‑1 (MQ=255)
    tGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAATTATTTca                                                                                        <  2:1042444/91‑1 (MQ=255)
    tGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTca                                                                                        >  1:1086435/1‑91 (MQ=255)
    tGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTca                                                                                        >  1:2556048/1‑91 (MQ=255)
      ttCAAATAATACTCGAAATTAAAGTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGATATGTACTATAATGTTTGTGAAATAATTcaca                                                                                      <  2:2961322/91‑1 (MQ=255)
       tCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACaa                                                                                     >  2:2042650/1‑91 (MQ=255)
       tCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACaa                                                                                     <  1:676737/91‑1 (MQ=255)
       tCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACaa                                                                                     <  1:2995359/91‑1 (MQ=255)
       tCAAACAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAATTATTTCACaa                                                                                     <  2:3148150/91‑1 (MQ=255)
         aaaTAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAg                                                                                   <  2:3564996/91‑1 (MQ=255)
          aataatATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTcac                                                                                       <  2:2425895/86‑1 (MQ=255)
          aataatATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTcac                                                                                       >  1:2425895/1‑86 (MQ=255)
           ataataTTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGta                                                                                 <  2:2735191/91‑1 (MQ=255)
             aataTTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGtata                                                                               <  2:788540/91‑1 (MQ=255)
               tatTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATaaa                                                                             <  2:962108/91‑1 (MQ=255)
               tatTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATaaa                                                                             <  1:1118074/91‑1 (MQ=255)
                atTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAg                                                                            <  2:2333075/91‑1 (MQ=255)
                atTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAg                                                                            <  2:91824/91‑1 (MQ=255)
                atTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAg                                                                            <  1:2494122/91‑1 (MQ=255)
                 ttCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAgg                                                                           <  2:1337945/91‑1 (MQ=255)
                 ttCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAgg                                                                           >  1:642487/1‑91 (MQ=255)
                  tCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGa                                                                          >  1:2226588/1‑91 (MQ=255)
                  tCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGa                                                                          <  2:2614366/91‑1 (MQ=255)
                   cGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAg                                                                         <  1:1684287/91‑1 (MQ=255)
                   cGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAg                                                                         <  1:3528588/91‑1 (MQ=255)
                   cGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAg                                                                         <  1:1875071/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)
                          aaaCTTTTTTTTTTATTTTTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGt                                                                  >  2:1250533/1‑91 (MQ=255)
                          aaaCTTTTTTGTATCGGCTTTGCTTTATATCCTGAGAGAGATGTACTATAATGGTTGTGAGATAATTCACAAAGTATAACGGAGTGGTTGt                                                                  <  2:2187838/91‑1 (MQ=255)
                           aaCTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGta                                                                 <  1:1832530/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:3812639/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)
                                       taGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACtat                                                     >  1:2022854/1‑91 (MQ=255)
                                         gAATTTTCTTTATATCCTGAGAGACATGTACTATAAAGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATAc                                                   <  2:1534395/91‑1 (MQ=255)
                                           aTTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCt                                                 <  1:2508579/91‑1 (MQ=255)
                                                 tttATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGaaaaa                                           <  1:2547084/91‑1 (MQ=255)
                                                  ttATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAg                                          <  1:461472/91‑1 (MQ=255)
                                                   tataTCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGa                                         >  1:1373219/1‑91 (MQ=255)
                                                   tataTCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGa                                         <  1:2379532/91‑1 (MQ=255)
                                                    atatCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGaa                                        <  2:1063674/91‑1 (MQ=255)
                                                     tatCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGaaa                                       >  1:1706486/1‑91 (MQ=255)
                                                      atCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACATGAAAAAGaaaa                                      <  2:2292332/91‑1 (MQ=255)
                                                        ccTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATc                                    >  1:3234828/1‑91 (MQ=255)
                                                          tGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGt                                  <  2:2824400/91‑1 (MQ=255)
                                                           gagagaCATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTg                                 <  2:806958/91‑1 (MQ=255)
                                                            agagaCATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTgg                                >  2:573366/1‑91 (MQ=255)
                                                            agagaCATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTgg                                >  1:976876/1‑91 (MQ=255)
                                                            agagaCATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTgg                                >  2:573318/1‑91 (MQ=255)
                                                             gagaCATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGa                               <  2:3328000/91‑1 (MQ=255)
                                                                aCATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCg                            <  2:1676374/91‑1 (MQ=255)
                                                                 cATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGa                           >  1:3420065/1‑91 (MQ=255)
                                                                 cATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGCTAACTATACCTGAAAAAGAAAATCGTGGATCGa                           >  1:3181180/1‑91 (MQ=255)
                                                                 aaTGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGa                           <  2:2527033/90‑1 (MQ=255)
                                                                  aTGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGaa                          >  2:1264551/1‑91 (MQ=255)
                                                                   tGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGaaa                         <  1:617889/91‑1 (MQ=255)
                                                                       cTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAAc                     <  2:1360587/91‑1 (MQ=255)
                                                                       cTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAAc                     <  1:260117/91‑1 (MQ=255)
                                                                       cTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAAc                     <  1:1286518/91‑1 (MQ=255)
                                                                       cTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAAc                     <  1:939238/91‑1 (MQ=255)
                                                                          taATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACaag                  <  1:2876814/91‑1 (MQ=255)
                                                                          taATGGTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACaag                  <  2:2877244/91‑1 (MQ=255)
                                                                           aaTGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACaaga                 >  1:2681305/1‑91 (MQ=255)
                                                                           aaTGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACaaga                 >  2:1207783/1‑91 (MQ=255)
                                                                           aaTGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACaaga                 <  2:3005913/91‑1 (MQ=255)
                                                                            aTGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACaagaa                <  1:1028134/91‑1 (MQ=255)
                                                                             tGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACaagaag               <  1:2232358/91‑1 (MQ=255)
                                                                               tttGTGAATTATTTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTg             >  2:1142586/1‑91 (MQ=255)
                                                                               tttGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTg             <  1:1985693/91‑1 (MQ=255)
                                                                                ttGTGAATTATTTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAACTCGTGGATCGAACGAACAAGAAGTgg            >  2:188275/1‑91 (MQ=255)
                                                                                ttGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTgg            >  2:532730/1‑91 (MQ=255)
                                                                                 tgtgAATTATTTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGc           >  2:2265536/1‑91 (MQ=255)
                                                                                 tgtgAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTg             <  2:485850/89‑1 (MQ=255)
                                                                                 tgtgAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGc           <  1:2172122/91‑1 (MQ=255)
                                                                                 tgtgAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGc           <  1:2646749/91‑1 (MQ=255)
                                                                                 tgtgAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATAGTGGATCGAAAGAACAAGAAGTGGc           >  1:3628446/1‑91 (MQ=255)
                                                                                 tgtgAAACAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACTAGAAGTGGc           <  1:1950881/91‑1 (MQ=255)
                                                                                   tgAATTATTTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCaa         >  1:2179439/1‑91 (MQ=255)
                                                                                   tgAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCaa         >  1:1140398/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)
                                                                                     aaaTAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAAtt       >  1:2687961/1‑91 (MQ=255)
                                                                                      aataatTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTa      >  1:722652/1‑91 (MQ=255)
                                                                                      aataatTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTa      <  1:1795292/91‑1 (MQ=255)
                                                                                       ataatTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTAt     <  2:2388799/91‑1 (MQ=255)
                                                                                        taatTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTAtg    >  2:931630/1‑91 (MQ=255)
                                                                                         aatTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTAtga   <  2:2204265/91‑1 (MQ=255)
                                                                                          atTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTAtgat  >  1:1134748/1‑91 (MQ=255)
                                                                                          |                                                                                          
TATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTATGAT  >  USA300TCH1516_ALE/170982‑171162

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: