Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,797,772 G→C 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,7720GC97.9% 708.6 / 3.7 195R7G (CGC→GGC) phoRAlkaline phosphatase synthesis sensor protein PhoR
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base C (78/112);  total (81/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 199 total aligned reads displayed.

GATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGG  >  USA300TCH1516_ALE/1797684‑1797858
                                                                                        |                                                                                      
   aaaaCATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCcgtggt                                                                                   <  1:1863675/91‑1 (MQ=255)
    aaaCATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCgtggtg                                                                                  >  2:2186238/1‑91 (MQ=255)
    aaaCATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCgtggtg                                                                                  >  2:1153232/1‑91 (MQ=255)
     aaCATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGa                                                                                 >  2:2715159/1‑91 (MQ=255)
     aaCATATGGGCTGGAAGCGTGAAGAAGAATGCCTAATGTAACCAAGCTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGa                                                                                 <  1:2996525/91‑1 (MQ=255)
      aCATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGaa                                                                                <  2:899175/91‑1 (MQ=255)
      aCATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGaa                                                                                <  1:100588/91‑1 (MQ=255)
        atatGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAAc                                                                              <  1:255242/91‑1 (MQ=255)
        atatGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAAc                                                                              <  1:3052364/91‑1 (MQ=255)
         tatGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACt                                                                             >  2:409384/1‑91 (MQ=255)
          atGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGa                                                                                 >  1:1891051/1‑86 (MQ=255)
          atGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACtt                                                                            >  1:2677356/1‑91 (MQ=255)
          atGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACtt                                                                            >  2:1418757/1‑91 (MQ=255)
          atGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACtt                                                                            >  1:3300666/1‑91 (MQ=255)
          atGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACtt                                                                            >  2:669492/1‑91 (MQ=255)
          atGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACtt                                                                            >  1:1335794/1‑91 (MQ=255)
             tgATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTtcat                                                                         >  2:185235/1‑91 (MQ=255)
              gATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTtcatc                                                                        >  2:1943126/1‑91 (MQ=255)
              gATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTtcatc                                                                        >  2:169614/1‑91 (MQ=255)
               aTGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTtcatca                                                                       <  1:1699585/91‑1 (MQ=255)
               aTGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTtcatca                                                                       <  1:1701064/91‑1 (MQ=255)
                 ggAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCAtt                                                                     <  2:2318820/91‑1 (MQ=255)
                  gAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGTCGTGGTGAAACTTCATCATTg                                                                    <  1:1785022/91‑1 (MQ=255)
                   aaGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGt                                                                   <  1:3418352/91‑1 (MQ=255)
                      cGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCt                                                                >  2:1264339/1‑91 (MQ=255)
                      cGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCt                                                                >  1:2988461/1‑91 (MQ=255)
                      cGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCt                                                                <  2:733687/91‑1 (MQ=255)
                        tGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCttt                                                              >  1:1206958/1‑91 (MQ=255)
                        tGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCttt                                                              <  1:1567901/91‑1 (MQ=255)
                            taataatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGt                                                          <  2:588863/91‑1 (MQ=255)
                            taataatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGt                                                          >  2:888498/1‑91 (MQ=255)
                            taataatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGt                                                          <  1:1008009/91‑1 (MQ=255)
                             aataatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTc                                                         <  1:1178201/91‑1 (MQ=255)
                               taatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCtt                                                       <  2:1236164/91‑1 (MQ=255)
                                 atGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTGGGTCTTTc                                                     <  2:691105/91‑1 (MQ=255)
                                 atGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTc                                                     <  1:232960/91‑1 (MQ=255)
                                 atGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTc                                                     <  2:1544503/91‑1 (MQ=255)
                                  tGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTcc                                                    >  2:3547294/1‑91 (MQ=255)
                                     cTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAAt                                                 <  2:3532262/91‑1 (MQ=255)
                                     cTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAAt                                                 <  2:3060222/91‑1 (MQ=255)
                                      tAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAAtt                                                >  2:2606514/1‑91 (MQ=255)
                                      tAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAAtt                                                <  1:2829319/91‑1 (MQ=255)
                                       aaTGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAAttt                                               <  1:1259457/91‑1 (MQ=255)
                                       aaTGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAAttt                                               <  2:2958304/91‑1 (MQ=255)
                                       aaTGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAAttt                                               <  2:2957019/91‑1 (MQ=255)
                                       aaTGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAAttt                                               >  2:1912417/1‑91 (MQ=255)
                                         tGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTtat                                             <  1:2020549/91‑1 (MQ=255)
                                           tAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAg                                           <  1:1710615/91‑1 (MQ=255)
                                           tAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAg                                           >  1:2004628/1‑91 (MQ=255)
                                            aaCAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGc                                          <  1:3552841/91‑1 (MQ=255)
                                            aaCAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGc                                          <  1:685720/91‑1 (MQ=255)
                                             aCAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGcc                                         >  2:2439938/1‑91 (MQ=255)
                                             aCAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGcc                                         >  2:3394258/1‑91 (MQ=255)
                                             aCAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGcc                                         >  1:2064456/1‑91 (MQ=255)
                                               aaaaCTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTa                                       <  1:1390995/91‑1 (MQ=255)
                                               aaaaCTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTa                                       >  2:2693535/1‑91 (MQ=255)
                                                  aCTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGc                                    <  1:1044622/91‑1 (MQ=255)
                                                  aCTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGc                                    <  2:1195352/91‑1 (MQ=255)
                                                  aCTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGc                                    >  1:249270/1‑91 (MQ=255)
                                                   cTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGcc                                   >  2:3027973/1‑91 (MQ=255)
                                                    tGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCa                                  <  1:1505524/91‑1 (MQ=255)
                                                     gATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCAc                                 <  1:130568/91‑1 (MQ=255)
                                                     gATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCAc                                 <  1:638484/91‑1 (MQ=255)
                                                      aTAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACg                                <  1:1643301/91‑1 (MQ=255)
                                                      aTAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACg                                <  2:1437759/91‑1 (MQ=255)
                                                      aTAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACg                                >  2:1871431/1‑91 (MQ=255)
                                                       taattaatATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGc                               <  2:105316/91‑1 (MQ=255)
                                                       taattaatATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGc                               >  2:2268373/1‑91 (MQ=255)
                                                        aattaatATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGca                              <  2:3292344/91‑1 (MQ=255)
                                                        aattaatATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGca                              <  1:2122527/91‑1 (MQ=255)
                                                           taatATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAAATTATAGCCTAAGCCACGCACAg                           >  2:2931226/1‑91 (MQ=255)
                                                             atatACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGtt                         <  2:2030659/91‑1 (MQ=255)
                                                               ataCTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGtttt                       <  1:1116292/91‑1 (MQ=255)
                                                                taCTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTta                      <  1:1183082/91‑1 (MQ=255)
                                                                taCTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTta                      <  2:776432/91‑1 (MQ=255)
                                                                taCTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTta                      <  2:2064630/91‑1 (MQ=255)
                                                                  cTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTGGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTtaat                    >  2:3194329/1‑91 (MQ=255)
                                                                  cTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTtaat                    >  2:200995/1‑91 (MQ=255)
                                                                   tGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTtaata                   >  2:1211979/1‑91 (MQ=255)
                                                                      ttATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGt                >  2:2076849/1‑91 (MQ=255)
                                                                      ttATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGt                <  1:3544014/91‑1 (MQ=255)
                                                                       tATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGtt               <  1:888498/91‑1 (MQ=255)
                                                                        aTCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTtg              <  1:33142/91‑1 (MQ=255)
                                                                         tCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTtgt             <  2:451833/91‑1 (MQ=255)
                                                                         tCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTtgt             >  2:2682934/1‑91 (MQ=255)
                                                                         tCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTtgt             <  1:1420515/91‑1 (MQ=255)
                                                                         tCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTtgt             <  1:3148609/91‑1 (MQ=255)
                                                                         tCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTtgt             >  1:773315/1‑91 (MQ=255)
                                                                            aTAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGc          >  1:1815711/1‑91 (MQ=255)
                                                                            aTAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGc          >  1:1813235/1‑91 (MQ=255)
                                                                                cATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGcttct      <  1:1210819/91‑1 (MQ=255)
                                                                                  ttAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTa    <  1:1682399/91‑1 (MQ=255)
                                                                                  ttAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTa    >  2:1954218/1‑91 (MQ=255)
                                                                                   tAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAg   <  2:2739935/91‑1 (MQ=255)
                                                                                    aaGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACTGTTTTAATAAGTTGTGGCTTCTTAgg  <  2:2064397/91‑1 (MQ=255)
                                                                                    aaGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAgg  >  2:2578947/1‑91 (MQ=255)
                                                                                    aaGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAgg  <  1:2890405/91‑1 (MQ=255)
                                                                                    aaGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAgg  <  2:3348092/91‑1 (MQ=255)
                                                                                    aaGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAgg  >  2:2645067/1‑91 (MQ=255)
                                                                                    aaGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAgg  >  1:827722/1‑91 (MQ=255)
                                                                                        |                                                                                      
GATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGG  >  USA300TCH1516_ALE/1797684‑1797858

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: