Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,797,772 G→C 11.1% R7G (CGC→GGC)  phoR ← Alkaline phosphatase synthesis sensor protein PhoR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,797,7720GC11.1% 759.7 / 66.1 235R7G (CGC→GGC) phoRAlkaline phosphatase synthesis sensor protein PhoR
Reads supporting (aligned to +/- strand):  ref base G (103/104);  new base C (16/10);  total (119/116)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 257 total aligned reads displayed.

CTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATT  >  USA300TCH1516_ALE/1797682‑1797861
                                                                                          |                                                                                         
cTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCg                                                                                           <  2:3960592/91‑1 (MQ=255)
cTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCg                                                                                           >  2:206147/1‑91 (MQ=255)
 tGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCgg                                                                                          >  1:1176076/1‑91 (MQ=255)
 tGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCgg                                                                                          <  2:783742/91‑1 (MQ=255)
 tGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCgg                                                                                          >  1:2913277/1‑91 (MQ=255)
 tGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCgg                                                                                          >  2:1987444/1‑91 (MQ=255)
   aTAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCggtg                                                                                        <  1:3539157/91‑1 (MQ=255)
   aTAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCggtg                                                                                        >  2:1437457/1‑91 (MQ=255)
   aTAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCggtg                                                                                        <  2:2386212/91‑1 (MQ=255)
     aaaaCATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCggtggt                                                                                      >  1:1132792/1‑91 (MQ=255)
      aaaCATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGgtggtg                                                                                     <  1:3414631/91‑1 (MQ=255)
      aaaCATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGgtggtg                                                                                     >  2:1292667/1‑91 (MQ=255)
       aaCATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGa                                                                                    <  1:2319254/91‑1 (MQ=255)
          gtatGTGATGGAAGCGTGAATAGTAATGCCTGATGTGACAAAACTGATAATTAATATACTGCGAAACAATAACATTAAGAGGTGGTGAAAc                                                                                 <  2:3306380/90‑1 (MQ=255)
          atatGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAAc                                                                                 <  1:3530014/91‑1 (MQ=255)
           tatGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACt                                                                                <  1:892645/91‑1 (MQ=255)
           tatGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACt                                                                                <  1:2864827/91‑1 (MQ=255)
           tatGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACt                                                                                >  1:1381610/1‑91 (MQ=255)
           tatGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTATCAATAACATTAAGCGGTGGTGAAACtt                                                                               <  2:897739/91‑1 (MQ=255)
            atGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACtt                                                                               >  2:2828320/1‑91 (MQ=255)
            atGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACtt                                                                               >  1:1692689/1‑91 (MQ=255)
             tgtgATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTtc                                                                              <  1:2632741/91‑1 (MQ=255)
             tgtgATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTtc                                                                              >  2:429775/1‑91 (MQ=255)
             tgtgATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTtc                                                                              >  2:2090603/1‑91 (MQ=255)
                gATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTtcatc                                                                           <  1:2710001/91‑1 (MQ=255)
                gATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTtcatc                                                                           <  2:1623951/91‑1 (MQ=255)
                 aTGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTtcatca                                                                          <  1:2406198/91‑1 (MQ=255)
                  tGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTGATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTcatcat                                                                         <  2:3549147/91‑1 (MQ=255)
                   ggcaGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAACAACATTAAGCGGTGGTGACACTTCATCAtt                                                                        >  2:2294018/4‑91 (MQ=255)
                   ggAAGCGTGAATAATAATGCCAAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCAtt                                                                        <  2:3341326/91‑1 (MQ=255)
                    gAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTg                                                                       <  1:2012991/91‑1 (MQ=255)
                     aaGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTCAGCCGTGGTGAAACTTCATCATTGt                                                                      >  1:3668086/1‑91 (MQ=255)
                     aaGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGt                                                                      >  1:56542/1‑91 (MQ=255)
                         gTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCtt                                                                  <  2:1506004/91‑1 (MQ=255)
                         gTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCtt                                                                  <  1:3746860/91‑1 (MQ=255)
                          tGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCttt                                                                 <  1:1935976/91‑1 (MQ=255)
                          tGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCttt                                                                 <  2:3707954/91‑1 (MQ=255)
                          tGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCttt                                                                 >  1:261712/1‑91 (MQ=255)
                           gAATCATCATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAATCGGTGGTGAAACTTCATCATTGTTCTCTa                                                                >  2:1646396/1‑91 (MQ=255)
                           gAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTa                                                                >  2:3181004/1‑91 (MQ=255)
                             ataataatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAgg                                                              <  2:174331/91‑1 (MQ=255)
                             ataataatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAgg                                                              >  2:3125871/1‑91 (MQ=255)
                             ataataatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAgg                                                              >  2:833483/1‑91 (MQ=255)
                               aataatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTc                                                            <  2:1068570/91‑1 (MQ=255)
                               aataatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTCTAGGTc                                                            <  1:3074382/91‑1 (MQ=255)
                                ataatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCt                                                           >  1:458462/1‑91 (MQ=255)
                                ataatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCt                                                           <  2:118027/91‑1 (MQ=255)
                                  aatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCttt                                                         >  2:2217045/1‑91 (MQ=255)
                                   atGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTc                                                        >  1:655711/1‑91 (MQ=255)
                                   atGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTc                                                        >  2:11499/1‑91 (MQ=255)
                                    tGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTcc                                                       <  2:1445041/91‑1 (MQ=255)
                                     gCCTAATGTAACAAAACTGATAATTAATATAGTGCTTATCAATAGCATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCa                                                      <  1:3158204/91‑1 (MQ=255)
                                     gCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACATCATCATTGTTCTTTAGGTCTTTCCa                                                      <  1:1806841/91‑1 (MQ=255)
                                      ccTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCaa                                                     <  2:3034643/91‑1 (MQ=255)
                                       cTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAAt                                                    <  2:1172924/91‑1 (MQ=255)
                                          aTGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTta                                                 <  1:534462/91‑1 (MQ=255)
                                          aTGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTta                                                 <  2:955453/91‑1 (MQ=255)
                                          aTGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTta                                                 <  1:2062532/91‑1 (MQ=255)
                                          aTGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTta                                                 <  1:1676837/91‑1 (MQ=255)
                                           tGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTtat                                                <  2:539673/91‑1 (MQ=255)
                                            gTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTtata                                               >  2:2968747/1‑91 (MQ=255)
                                             tAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAg                                              >  1:2331947/1‑91 (MQ=255)
                                              aaCAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGc                                             <  1:929287/91‑1 (MQ=255)
                                              aaCAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGc                                             <  1:3489592/91‑1 (MQ=255)
                                               aCAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGcc                                            <  1:1010585/91‑1 (MQ=255)
                                               aCAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGcc                                            <  1:1329866/91‑1 (MQ=255)
                                               aCAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGcc                                            >  1:2317849/1‑91 (MQ=255)
                                                cAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCt                                           >  1:2294176/1‑91 (MQ=255)
                                                   aaCTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTATTTAGGTCTTTCCAATTTATAGCCTAAg                                        >  1:2250874/1‑91 (MQ=255)
                                                   aaCTGATAATTAATATACTGCTTATCAATAACATTAAGCGCTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTatg                                        >  2:3247395/1‑89 (MQ=255)
                                                       gATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCAc                                    <  2:1941427/91‑1 (MQ=255)
                                                       gATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCAc                                    <  2:1603669/91‑1 (MQ=255)
                                                       gATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCAc                                    <  1:1298134/91‑1 (MQ=255)
                                                       gATAATAAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCAc                                    >  1:3077103/1‑91 (MQ=255)
                                                        aTACTTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACg                                   <  1:3194134/91‑1 (MQ=255)
                                                        aTAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACg                                   >  2:2880669/1‑91 (MQ=255)
                                                        aTAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACg                                   <  1:3328467/91‑1 (MQ=255)
                                                         taattaatATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGc                                  <  1:1454282/91‑1 (MQ=255)
                                                          aattaatATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGca                                 <  1:1067186/91‑1 (MQ=255)
                                                          aattaatATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGca                                 <  1:2130859/91‑1 (MQ=255)
                                                               atatACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGtt                            >  1:1479657/1‑91 (MQ=255)
                                                                tataCTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGttt                           >  2:2424847/1‑91 (MQ=255)
                                                                 ataCTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGtttt                          >  1:178110/1‑91 (MQ=255)
                                                                  taCTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTta                         <  2:1779761/91‑1 (MQ=255)
                                                                     tGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTtaata                      >  1:2771304/1‑91 (MQ=255)
                                                                       cTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAg                    <  1:3257928/91‑1 (MQ=255)
                                                                        ttATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGt                   <  1:1844527/91‑1 (MQ=255)
                                                                        gtATCAATGACATGAAGCGGTGGTGAAACTGCGTCGTTGTGCTTTAGGTCTGTCGGATTTGTAGCCTGAGCCACGCACAGTTTTAATAAGt                   <  2:3724574/90‑1 (MQ=255)
                                                                          aTCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTtg                 <  2:3433360/91‑1 (MQ=255)
                                                                          aTCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTtg                 <  1:750163/91‑1 (MQ=255)
                                                                           tCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTtgt                <  2:613291/91‑1 (MQ=255)
                                                                               tAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGct            <  2:3387551/91‑1 (MQ=255)
                                                                               tAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGct            >  2:3598781/1‑91 (MQ=255)
                                                                                aaCATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGctt           >  2:1670498/1‑91 (MQ=255)
                                                                                 aCATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGcttc          >  2:496002/1‑91 (MQ=255)
                                                                                  cATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGcttct         <  2:385531/91‑1 (MQ=255)
                                                                                   aTTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGcttctt        <  1:2395346/91‑1 (MQ=255)
                                                                                   aTTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGcttctt        >  2:892319/1‑91 (MQ=255)
                                                                                      aaGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAgg     <  1:988832/91‑1 (MQ=255)
                                                                                      aaGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAgg     >  2:2208758/1‑91 (MQ=255)
                                                                                          |                                                                                         
CTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATT  >  USA300TCH1516_ALE/1797682‑1797861

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: