Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,386,589 A→T 5.2% intergenic (‑130/‑32) USA300HOU_RS12250 ← / → USA300HOU_RS12255 hypothetical protein/putative HTH‑type transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,386,5890AT 5.2% 706.1 / 16.1 208intergenic (‑130/‑32)USA300HOU_RS12250/USA300HOU_RS12255hypothetical protein/putative HTH‑type transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (90/107);  new base T (9/2);  total (99/109)
Fisher's exact test for biased strand distribution p-value = 2.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 248 total aligned reads displayed.

AAAGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTATTTAAAATT  >  USA300TCH1516_ALE/2386500‑2386679
                                                                                         |                                                                                          
aaaGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTagctaa                                                                                           <  2:2845345/91‑1 (MQ=255)
  aGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGttagcttagc                                                                                         <  2:2643653/91‑1 (MQ=255)
  aGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGttagcttagc                                                                                         <  1:3043847/91‑1 (MQ=255)
   gTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTtagcttagct                                                                                        <  1:1798016/91‑1 (MQ=255)
   gTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTagctaagct                                                                                        <  1:3138636/91‑1 (MQ=255)
    tataATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTagcttagcta                                                                                       <  1:2441345/91‑1 (MQ=255)
    tataATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTagcttagcta                                                                                       <  2:3662888/91‑1 (MQ=255)
     ataATTATTTTTGCTCTTAATAATCATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAgcttagctaa                                                                                      <  2:3602206/91‑1 (MQ=255)
      taATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAAc                                                                                     <  2:1449928/91‑1 (MQ=255)
         ttattTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTAt                                                                                  >  2:3349397/1‑91 (MQ=255)
          tattTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTAtt                                                                                 <  2:605808/91‑1 (MQ=255)
          tattTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTAtt                                                                                 <  2:3746702/91‑1 (MQ=255)
          tattTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTAtt                                                                                 <  2:2109520/91‑1 (MQ=255)
           attTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTAttt                                                                                <  1:1051139/91‑1 (MQ=255)
           attTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTAttt                                                                                >  1:2605236/1‑91 (MQ=255)
            tttttGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACt                                                                                    <  2:223293/86‑1 (MQ=255)
             ttttGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTg                                                                              <  2:2340381/91‑1 (MQ=255)
               ttGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTc                                                                            <  2:2140469/91‑1 (MQ=255)
                 gCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCtt                                                                          <  1:2644387/91‑1 (MQ=255)
                 gCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCtt                                                                          <  2:591819/91‑1 (MQ=255)
                  ctctTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTa                                                                         >  2:1347647/1‑91 (MQ=255)
                  ctctTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTa                                                                         >  2:3528579/1‑91 (MQ=255)
                   tctTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAg                                                                        <  1:2344103/91‑1 (MQ=255)
                   tctTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAg                                                                        >  1:1350220/1‑91 (MQ=255)
                   tctTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAg                                                                        >  1:1186869/1‑91 (MQ=255)
                       aataatTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTGTCTTAGGAGAt                                                                   <  1:1538816/91‑1 (MQ=255)
                        ataatTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGAt                                                                   <  1:631250/91‑1 (MQ=255)
                          aatTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGt                                                                 <  1:3487021/91‑1 (MQ=255)
                          aatTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGt                                                                 <  1:2898441/91‑1 (MQ=255)
                          aatTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGt                                                                 <  1:2170573/91‑1 (MQ=255)
                            ttATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCa                                                               <  2:692320/91‑1 (MQ=255)
                            ttATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCa                                                               <  2:2010795/91‑1 (MQ=255)
                                atatTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGt                                                                 >  2:342881/1‑85 (MQ=255)
                                atatTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTAt                                                           <  2:2826500/91‑1 (MQ=255)
                                atatTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTAt                                                           >  2:573430/1‑91 (MQ=255)
                                 tatTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATg                                                          <  1:76894/91‑1 (MQ=255)
                                        ttCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATcac                                                   <  2:1639802/91‑1 (MQ=255)
                                         tCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATcaca                                                  <  1:301210/91‑1 (MQ=255)
                                          cTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACaa                                                 <  1:2256594/91‑1 (MQ=255)
                                           tAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAg                                                <  1:1146464/91‑1 (MQ=255)
                                           tAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAg                                                >  1:2308079/1‑91 (MQ=255)
                                            aaTTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGa                                               <  2:3741458/91‑1 (MQ=255)
                                             aTTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGaa                                              >  1:3215507/1‑91 (MQ=255)
                                             aTTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGaa                                              >  2:803415/1‑91 (MQ=255)
                                              ttGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAAt                                             <  2:3265227/91‑1 (MQ=255)
                                               tGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAAtt                                            <  1:718665/91‑1 (MQ=255)
                                               tGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAAtt                                            >  1:333753/1‑91 (MQ=255)
                                                gCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTagctaagct                                                                                        >  1:1917049/1‑46 (MQ=255)
                                                gCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAAttt                                           <  1:124266/91‑1 (MQ=255)
                                                 cTCCCGTTTTATAATTCCTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAAtttt                                          >  2:174570/1‑91 (MQ=255)
                                                 cTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAAtttt                                          <  1:2499066/91‑1 (MQ=255)
                                                 cTCCCGTTTTATAATTACTATTGTTGTAAAAAGCTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAAtttt                                          >  1:2724738/1‑91 (MQ=255)
                                                   cccGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTc                                        <  2:3215507/91‑1 (MQ=255)
                                                     cGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCaa                                      >  1:3425808/1‑91 (MQ=255)
                                                      gTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAAt                                     >  1:2205576/1‑91 (MQ=255)
                                                         ttATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGt                                  >  2:190104/1‑91 (MQ=255)
                                                          tataATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGtt                                 <  1:1691228/91‑1 (MQ=255)
                                                          tataATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGtt                                 >  2:226404/1‑91 (MQ=255)
                                                          tataATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGtt                                 >  1:1117182/1‑91 (MQ=255)
                                                           ataATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGttt                                <  1:1572633/91‑1 (MQ=255)
                                                             aaTTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTta                              >  2:3841712/1‑91 (MQ=255)
                                                             aaTTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTta                              <  1:40059/91‑1 (MQ=255)
                                                              aTTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTtat                             <  1:1398705/91‑1 (MQ=255)
                                                                tACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATaa                           <  1:2544000/91‑1 (MQ=255)
                                                                 aCTATTGTTGTAAAAAGCTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAAc                          <  2:2724738/91‑1 (MQ=255)
                                                                   tATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACaa                        <  2:3084335/91‑1 (MQ=255)
                                                                   tATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACaa                        <  2:3084100/91‑1 (MQ=255)
                                                                    aTTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAat                       >  2:157969/1‑91 (MQ=255)
                                                                    aTTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAat                       <  2:1413185/91‑1 (MQ=255)
                                                                    aTTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAat                       >  1:679356/1‑91 (MQ=255)
                                                                    aTTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAat                       <  1:2355261/91‑1 (MQ=255)
                                                                     ttgttgTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAata                      >  2:3109030/1‑91 (MQ=255)
                                                                        ttgtAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAatatat                   <  2:507852/91‑1 (MQ=255)
                                                                         tgtAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATc                  <  2:1080952/91‑1 (MQ=255)
                                                                         tgtAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATc                  >  2:412605/1‑91 (MQ=255)
                                                                         tgtAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATc                  <  1:516446/91‑1 (MQ=255)
                                                                          gtAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCg                 <  1:890104/91‑1 (MQ=255)
                                                                          gtAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCg                 >  1:3853117/1‑91 (MQ=255)
                                                                              aatagTTAGCTTAGCTAACCTTTTTGTCTTAGGAGATGTCACTATGTTATAACAAGACTACTCCCATAGTTTTATAACAATATATCGcccc             >  2:1741625/5‑91 (MQ=19)
                                                                              aaaGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGccca             >  2:3851213/1‑90 (MQ=255)
                                                                              aaaGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTACGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGcccc             <  1:3090090/91‑1 (MQ=255)
                                                                               atagTTAGCTTAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCt            >  1:1122688/4‑91 (MQ=255)
                                                                               aaGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCt            <  1:1214899/91‑1 (MQ=255)
                                                                               aaGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCt            >  1:3738521/1‑91 (MQ=255)
                                                                                tagTTAGCTTAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTa           >  1:3565073/3‑91 (MQ=255)
                                                                                 ggTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTAt          >  2:2425812/1‑91 (MQ=255)
                                                                                 ggTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTAt          >  2:76908/1‑91 (MQ=255)
                                                                                 ggTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTAt          >  2:1909497/1‑91 (MQ=255)
                                                                                 agTTAGCTTAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTAt          >  2:85764/2‑91 (MQ=255)
                                                                                  gTTAGCTTAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTAtt         >  2:1025882/1‑91 (MQ=255)
                                                                                  gTTAGCTTAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTAtt         >  2:1974297/1‑91 (MQ=255)
                                                                                   ttAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTAttt        >  2:810916/1‑91 (MQ=255)
                                                                                      gctaagctaaCTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTATTTaaa     >  2:3627909/1‑91 (MQ=255)
                                                                                      gctaagctaaCTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTATTTaaa     >  2:111343/1‑91 (MQ=255)
                                                                                       ctaagctaaCTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTATTTaaaa    >  1:3298010/1‑91 (MQ=255)
                                                                                       ctaagctaaCTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTATTTaaaa    >  1:2953916/1‑91 (MQ=255)
                                                                                        taagctaaCTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTATTTAAAAt   <  1:360109/91‑1 (MQ=255)
                                                                                        taagctaaCTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTATTTAAAAt   <  1:2068075/91‑1 (MQ=255)
                                                                                         aagctaaCTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTATTTAAAAtt  <  2:2308079/91‑1 (MQ=255)
                                                                                         aagctaaCTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTATTTAAAAtt  >  2:2516597/1‑91 (MQ=255)
                                                                                         |                                                                                          
AAAGTATAATTATTTTTGCTCTTAATAATTATATATTGATTTCTAATTGCTCCCGTTTTATAATTACTATTGTTGTAAAAAGGTTAGCTAAGCTAACTATTTTGTCTTAGGAGATGTCACTATGCTATCACAAGAATTTTTCAATAGTTTTATAACAATATATCGCCCCTATTTAAAATT  >  USA300TCH1516_ALE/2386500‑2386679

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