Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 593,433 T→A 7.4% intergenic (+25/‑112) rpoC → / → rplGB DNA‑directed RNA polymerase subunit beta'/Ribosome‑associated protein L7Ae‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE593,4330TA 7.4% 738.1 / 35.2 226intergenic (+25/‑112)rpoC/rplGBDNA‑directed RNA polymerase subunit beta'/Ribosome‑associated protein L7Ae‑like protein
Reads supporting (aligned to +/- strand):  ref base T (116/93);  new base A (6/11);  total (122/104)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

Only 100 of 234 total aligned reads displayed.

GTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGC  >  USA300TCH1516_ALE/593343‑593523
                                                                                           |                                                                                          
gTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCt                                                                                            >  1:384651/1‑91 (MQ=255)
 tAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCtt                                                                                           <  2:2424735/91‑1 (MQ=255)
 tAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCtt                                                                                           <  2:3214933/91‑1 (MQ=255)
   aaaTACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTa                                                                                         >  1:2185580/1‑91 (MQ=255)
   aaaTACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTa                                                                                         >  1:2922654/1‑91 (MQ=255)
    aaTACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAg                                                                                        <  2:3194027/91‑1 (MQ=255)
      tACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGcc                                                                                      >  1:2105971/1‑91 (MQ=255)
        cGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAAAGCATTAGCctc                                                                                    <  1:1280756/91‑1 (MQ=255)
          aaaaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCtt                                                                                  <  2:404035/91‑1 (MQ=255)
           aaaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAAGAGGCTAAAGCATTAGCctctg                                                                                  <  2:1917178/91‑2 (MQ=38)
           aaaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTg                                                                                 >  1:3290340/1‑91 (MQ=255)
           aaaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTg                                                                                 <  1:1716334/91‑1 (MQ=255)
            aaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAAGAGGCTAAAGCATTAGCctctgt                                                                                 <  1:3565482/91‑3 (MQ=38)
            aaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGt                                                                                <  2:1666020/91‑1 (MQ=255)
            aaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGt                                                                                <  2:1066389/91‑1 (MQ=255)
             aaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGtt                                                                               <  2:971364/91‑1 (MQ=255)
             aaaCAGCTAAACCAGTTGCAGAAGTTGAATCGCAAACTGAAGTAACGGAATAACAAGTATATAACAAGAGGCTAAAGCATTAGCctctgtt                                                                                <  2:9355/91‑4 (MQ=25)
              aaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAAGAGGCTAAAGCATTAGCctctgtta                                                                               <  1:40688/91‑5 (MQ=37)
               aCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAAAAGAGGCTAAAGCATTAGCCTCTGTTAta                                                                            <  1:1608154/91‑2 (MQ=37)
                  gCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAAGAGGCTAAAGCATTAGCCTCTGTTAtata                                                                          <  2:762964/91‑4 (MQ=21)
                  gCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTAtttt                                                                          >  2:1701414/1‑91 (MQ=255)
                  gCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTAtttt                                                                          >  2:3399016/1‑91 (MQ=255)
                   cTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTAttttt                                                                         >  1:104275/1‑91 (MQ=255)
                   cTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTAttttt                                                                         >  2:3792404/1‑91 (MQ=255)
                    tAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTa                                                                        >  1:36444/1‑91 (MQ=255)
                    tAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTCTTa                                                                        >  2:2930492/1‑91 (MQ=255)
                      aaCCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATg                                                                      >  2:1118814/1‑91 (MQ=255)
                        ccAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTa                                                                    >  1:1875522/1‑91 (MQ=255)
                         cAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTaa                                                                   <  1:569033/91‑1 (MQ=255)
                          aGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTaaa                                                                  >  1:3747025/1‑91 (MQ=255)
                           gTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAat                                                                 <  1:1338190/91‑1 (MQ=255)
                           gTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAat                                                                 <  1:344464/91‑1 (MQ=255)
                            ttGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAatt                                                                >  1:2460597/1‑91 (MQ=255)
                            ttGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAatt                                                                >  1:2460104/1‑91 (MQ=255)
                               cAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAattatt                                                             >  1:1074348/1‑91 (MQ=255)
                               cAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAattatt                                                             <  1:3062903/91‑1 (MQ=255)
                                aGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTAttt                                                            >  1:3101726/1‑91 (MQ=255)
                                 gAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTg                                                           >  2:3255427/1‑91 (MQ=255)
                                  aaGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGa                                                          <  2:2025594/91‑1 (MQ=255)
                                    gTTGAATCTGAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGAtt                                                        <  1:3707716/91‑1 (MQ=255)
                                     ttGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGAttt                                                       <  2:2573203/91‑1 (MQ=255)
                                     ttGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGAttt                                                       >  2:1085719/1‑91 (MQ=255)
                                        aaTCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAAt                                                    >  1:1177150/1‑91 (MQ=255)
                                        aaTCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAAt                                                    <  1:2308980/91‑1 (MQ=255)
                                         aTCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATg                                                   <  2:2710180/91‑1 (MQ=255)
                                          tctcAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGt                                                  <  2:856150/91‑1 (MQ=255)
                                             cAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGa                                               >  2:218328/1‑91 (MQ=255)
                                             cAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGa                                               >  1:2551797/1‑91 (MQ=255)
                                             cAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGa                                               >  1:3064568/1‑91 (MQ=255)
                                               aaCTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACg                                             >  2:642023/1‑91 (MQ=255)
                                                aCTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGa                                            >  2:3880505/1‑91 (MQ=255)
                                                  tGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAAt                                          <  1:1324092/91‑1 (MQ=255)
                                                  tGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAAt                                          <  2:3120111/91‑1 (MQ=255)
                                                      gTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATtctc                                      >  2:2788147/1‑91 (MQ=255)
                                                      gTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATtctc                                      >  2:1399528/1‑91 (MQ=255)
                                                      gTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATtctc                                      >  2:2299408/1‑91 (MQ=255)
                                                      gTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATtctc                                      >  2:2719046/1‑91 (MQ=255)
                                                         aCGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTg                                   <  1:1807554/91‑1 (MQ=255)
                                                           ggAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGtt                                 <  1:2075597/91‑1 (MQ=255)
                                                             aaTAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCa                               <  2:1177150/91‑1 (MQ=255)
                                                             aaTAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCa                               >  2:1254437/1‑91 (MQ=255)
                                                              aTAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCaa                              <  1:1742220/91‑1 (MQ=255)
                                                               tAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAAt                             >  2:3040600/1‑91 (MQ=255)
                                                               tAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTAATTTAATGTTGACGAATTCTCTTGTTCAAt                             >  1:2702568/1‑91 (MQ=255)
                                                                aaCAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATg                            >  1:1968326/1‑91 (MQ=255)
                                                                aaCAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATg                            >  2:2519024/1‑91 (MQ=255)
                                                                aaCAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATg                            <  1:1461732/91‑1 (MQ=255)
                                                                    aGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTaa                        >  2:2225814/1‑91 (MQ=255)
                                                                      tatataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAata                      >  1:2318988/1‑91 (MQ=255)
                                                                       atataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAatat                     >  2:2797975/1‑91 (MQ=255)
                                                                       atataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAatat                     >  2:3346846/1‑91 (MQ=255)
                                                                       atataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAatata                    <  1:3164504/91‑1 (MQ=255)
                                                                        tataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAatata                    <  2:3793630/91‑1 (MQ=255)
                                                                        aataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAatata                    <  1:728846/90‑1 (MQ=255)
                                                                         ataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAatatat                   <  1:3346846/91‑1 (MQ=255)
                                                                         ataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAatatat                   <  1:789477/91‑1 (MQ=255)
                                                                         ataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAatatat                   >  2:1597866/1‑91 (MQ=255)
                                                                         ataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAatatat                   >  1:2601041/1‑91 (MQ=255)
                                                                          taACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATAtt                  >  2:3512871/1‑91 (MQ=255)
                                                                           aacaagagGCTAAAGCATTAGCCTCTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTa                 >  2:3658268/5‑91 (MQ=37)
                                                                           aaCAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTa                 >  1:3002822/1‑91 (MQ=255)
                                                                           aaCAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTa                 <  1:873229/91‑1 (MQ=255)
                                                                           aaCAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTa                 <  1:873082/91‑1 (MQ=255)
                                                                           aaCAGAGGCTAATGCTTTAGCCTCTTGTTAATTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTa                 >  1:1701210/1‑91 (MQ=255)
                                                                           aaCAGAGGCTAATGCTATAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTa                 <  2:2979687/91‑1 (MQ=255)
                                                                            aCAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTaa                >  2:2621397/1‑91 (MQ=255)
                                                                            aCAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTaa                >  2:1854596/1‑91 (MQ=255)
                                                                             cAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTaaa               <  2:3345749/91‑1 (MQ=255)
                                                                              aagagGCTAAAGCATTAGCCTCTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAg              >  2:291297/2‑91 (MQ=38)
                                                                                gagGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAgg             <  2:571422/91‑1 (MQ=255)
                                                                                 agGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGt            <  2:1417825/91‑1 (MQ=255)
                                                                                   gCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTg          >  1:869213/1‑91 (MQ=255)
                                                                                   gCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTg          <  1:1647216/91‑1 (MQ=255)
                                                                                     tAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGt            <  1:2795189/87‑1 (MQ=255)
                                                                                     tAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGAt        >  2:2341101/1‑91 (MQ=255)
                                                                                      aaTGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATg       >  2:2515492/1‑91 (MQ=255)
                                                                                        tGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCa     >  2:790443/1‑91 (MQ=255)
                                                                                        tGCTTTAGCCTCCTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCa     >  1:2968163/1‑91 (MQ=255)
                                                                                           tttAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGc  <  2:3719934/91‑1 (MQ=255)
                                                                                           tttAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGc  <  1:919702/91‑1 (MQ=255)
                                                                                           |                                                                                          
GTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGC  >  USA300TCH1516_ALE/593343‑593523

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