Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,218,730 G→T 7.1% intergenic (+221/+42) USA300HOU_RS06065 → / ← USA300HOU_RS06070 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,218,7300GT 7.1% 795.2 / 30.4 232intergenic (+221/+42)USA300HOU_RS06065/USA300HOU_RS06070hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (116/98);  new base T (7/10);  total (123/109)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

TAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATAC  >  USA300TCH1516_ALE/1218640‑1218820
                                                                                          |                                                                                           
tAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAgg                                                                                             >  2:2595824/1‑91 (MQ=255)
tAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAgg                                                                                             <  1:377420/91‑1 (MQ=255)
tAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAgg                                                                                             <  1:2301818/91‑1 (MQ=255)
  aCTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGac                                                                                           <  2:514073/91‑1 (MQ=255)
  aCTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGac                                                                                           <  1:2572307/91‑1 (MQ=255)
   cTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTATCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGact                                                                                          >  2:3198227/1‑91 (MQ=255)
     tAAATGATGAAACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGactac                                                                                        <  2:2674063/91‑1 (MQ=255)
      aaaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGactact                                                                                       <  1:2098722/91‑1 (MQ=255)
       aaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActt                                                                                      <  1:1060237/91‑1 (MQ=255)
        atgatgATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActtt                                                                                     <  1:691039/91‑1 (MQ=255)
         tgatgaTACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttc                                                                                    <  1:581553/91‑1 (MQ=255)
         tgatgaTACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttc                                                                                    >  1:2974633/1‑91 (MQ=255)
          gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttct                                                                                   >  2:1459732/1‑91 (MQ=255)
           atgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttctt                                                                                  >  2:941556/1‑91 (MQ=255)
            tgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTACCTAATAAATATACACCTCACCATAAGAAAGGACTActttcttt                                                                                <  2:116995/91‑1 (MQ=255)
            tgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttcttt                                                                                >  1:3235590/1‑91 (MQ=255)
            tgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttcttt                                                                                <  2:647186/91‑1 (MQ=255)
             gatACGTTATTTAGAGACACGAGACTTCGTGTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAATGTAGTCCTTTCTTAtg                                                                              <  2:2803184/91‑1 (MQ=12)
             gatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCtttt                                                                               <  2:603255/91‑1 (MQ=255)
              atAAGTTATTTAAAAAAACGATACTTCGTTTGAATAAACGCATTAAATAATAAATAAACCCCTCACAAAAGAAAGTAGTCCTTTCTTAtgg                                                                             <  2:2111485/91‑2 (MQ=12)
               tACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTtgt                                                                             >  2:2606153/1‑91 (MQ=255)
                aCGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTtgtg                                                                            <  2:1979975/91‑1 (MQ=255)
                aCGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTtgtg                                                                            >  1:1688173/1‑91 (MQ=255)
                aCGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTtgtg                                                                            >  1:2196706/1‑91 (MQ=255)
                aCGTTATTTAAAAACACGATACTTCGTTTCAAAGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTtgtg                                                                            >  2:3801020/1‑91 (MQ=255)
                 cGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGa                                                                           <  1:195358/91‑1 (MQ=255)
                  gTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAg                                                                          >  1:1755352/1‑91 (MQ=255)
                   ttATTTAAAAGCGCTATACTTCGTGTCAATGAACGCATTAAATCATAAATGAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAgg                                                                         <  2:3603414/91‑1 (MQ=255)
                   ttATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAgg                                                                         >  1:3361492/1‑91 (MQ=255)
                   ttATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCATTTGTGAgg                                                                         >  2:3502686/1‑91 (MQ=255)
                     aTTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtg                                                                       >  1:656775/1‑91 (MQ=255)
                      tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGactact                                                                                       <  1:480691/75‑1 (MQ=255)
                      tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                                      <  1:2672189/91‑1 (MQ=255)
                      tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                                      <  2:2183721/91‑1 (MQ=255)
                         aaaaaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTa                                                                   >  2:3103986/1‑91 (MQ=255)
                           aaaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAtt                                                                 >  2:1435128/1‑91 (MQ=255)
                           aaaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAtt                                                                 >  1:3302540/1‑91 (MQ=255)
                           aaaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAtt                                                                 >  1:2874496/1‑91 (MQ=255)
                           aaaCACGATACTTCGTGTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGTAAGGACTACTTTCTTTTGTGAGGTGTTTAtt                                                                 <  2:2970463/91‑1 (MQ=255)
                                cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                            <  1:3198602/91‑1 (MQ=255)
                                 gATACTTCGTTTCAATGAACGCATTAAATAATAAATGAACACCTCACCATGCGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGa                                                           <  2:1701818/91‑1 (MQ=255)
                                 gATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGa                                                           >  1:2249437/1‑91 (MQ=255)
                                  aTACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGAc                                                          >  1:1354463/1‑91 (MQ=255)
                                   tACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACa                                                         >  2:2035189/1‑91 (MQ=255)
                                   tACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACa                                                         >  2:2224930/1‑91 (MQ=255)
                                    ccTTCGTTTCAATGAACACATTAAATAATAAATAAACACCTCACCATAAGAAAGGCCTACTTTCTTTTGTGAGGTGTTTCTTTGTCGACaa                                                        >  2:759542/2‑91 (MQ=255)
                                    aCTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACaa                                                        >  1:751503/1‑91 (MQ=255)
                                       tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                                     >  1:1002745/1‑91 (MQ=255)
                                        cGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTg                                                    >  2:1495819/1‑91 (MQ=255)
                                          tttCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGAt                                                  >  2:325533/1‑91 (MQ=255)
                                          tttCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGAt                                                  >  1:3137684/1‑91 (MQ=255)
                                           ttCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGAtt                                                 <  2:1231463/91‑1 (MQ=255)
                                            tCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTa                                                <  2:2272389/91‑1 (MQ=255)
                                            tCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTATCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTa                                                <  1:3283508/91‑1 (MQ=255)
                                             cAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTAt                                               <  2:3176564/91‑1 (MQ=255)
                                             cAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTAt                                               >  1:1617134/1‑91 (MQ=255)
                                               aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                             >  1:2947183/1‑91 (MQ=255)
                                                  aaCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTg                                          >  1:576737/1‑91 (MQ=255)
                                                   aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                         <  1:982641/91‑1 (MQ=255)
                                                   aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                         >  2:2857027/1‑91 (MQ=255)
                                                   aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                         <  2:1354463/91‑1 (MQ=255)
                                                    cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                        <  1:1110520/91‑1 (MQ=255)
                                                    cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                        <  1:2512403/91‑1 (MQ=255)
                                                    cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                        <  2:2973184/91‑1 (MQ=255)
                                                    cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                        <  1:1458130/91‑1 (MQ=255)
                                                      cATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTc                                      <  2:3637281/91‑1 (MQ=255)
                                                       aTTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCATCTGTGAGGAGTTTATTAGTCGACAACTTGATTATGATTGTTTCa                                     >  2:1811735/1‑91 (MQ=255)
                                                         taaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtt                                   >  1:599926/1‑91 (MQ=255)
                                                           aataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtttt                                 >  1:1615408/1‑91 (MQ=255)
                                                           aataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtttt                                 <  1:1415530/91‑1 (MQ=255)
                                                           aataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtttt                                 <  1:719202/91‑1 (MQ=255)
                                                            ataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttttt                                >  2:2682939/1‑91 (MQ=255)
                                                              aataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTg                              <  1:2883479/91‑1 (MQ=255)
                                                               atagataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                             <  1:3494475/91‑1 (MQ=255)
                                                               ataaataaACACCCCACCAGAAGAAAGGACTACGTTCTTTTTTGAGGGGTTTATTTGGCGACAACTGGATTATGATTGTTTCATTTTGTGa                             <  2:173196/91‑1 (MQ=255)
                                                                taaataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGaa                            >  1:859416/1‑91 (MQ=255)
                                                                  aataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATa                          >  1:2642351/1‑91 (MQ=255)
                                                                  aataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATa                          >  2:885273/1‑91 (MQ=255)
                                                                    taaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAg                        >  1:2576072/1‑91 (MQ=255)
                                                                      aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt                      >  2:2243494/1‑91 (MQ=255)
                                                                       acacCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTc                     >  1:2762401/1‑91 (MQ=255)
                                                                        cacCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCa                    >  1:2134192/1‑91 (MQ=17)
                                                                         acCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCAt                   >  2:1131647/1‑91 (MQ=255)
                                                                         acCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCAt                   >  2:478306/1‑91 (MQ=255)
                                                                          ccTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATa                  <  1:3836432/91‑1 (MQ=255)
                                                                          ccTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATa                  >  1:3734071/1‑91 (MQ=255)
                                                                          ccTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATa                  >  2:2388750/1‑91 (MQ=255)
                                                                           cTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATaa                 >  2:3166741/1‑91 (MQ=255)
                                                                              aCCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCa              >  2:3876511/1‑91 (MQ=255)
                                                                              aCCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCa              <  1:956315/91‑1 (MQ=255)
                                                                               ccATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGaa                            <  2:255837/76‑1 (MQ=255)
                                                                               ccATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg             >  2:1525245/1‑91 (MQ=255)
                                                                                caaaaGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg             >  2:3729774/4‑91 (MQ=21)
                                                                                 aTAAGAGAGGACTACTTTCTTTTGTGAGGTGTTTATCTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgt           <  2:3835192/91‑1 (MQ=255)
                                                                                  tAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtg          <  2:576737/91‑1 (MQ=255)
                                                                                  tAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtg          >  1:1241385/1‑91 (MQ=255)
                                                                                    aGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTa        <  2:3750531/91‑1 (MQ=255)
                                                                                         ggACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATa   >  1:266788/1‑91 (MQ=255)
                                                                                         ggACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATa   <  1:2820557/91‑1 (MQ=255)
                                                                                          gACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATAc  >  2:3693712/1‑91 (MQ=255)
                                                                                          |                                                                                           
TAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATAC  >  USA300TCH1516_ALE/1218640‑1218820

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