Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,918,555 G→A 5.3% intergenic (‑273/+225) USA300HOU_RS07235 ← / ← USA300HOU_RS09510 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,918,5550GA 5.3% 638.8 / 15.3 189intergenic (‑273/+225)USA300HOU_RS07235/USA300HOU_RS09510hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (98/81);  new base A (2/8);  total (100/89)
Fisher's exact test for biased strand distribution p-value = 4.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 201 total aligned reads displayed.

GTAGAGATGCATTTTTTGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTACAACAATGTACATTTAAGGTGA  >  USA300TCH1516_ALE/1918465‑1918645
                                                                                           |                                                                                          
gTAGAGATGCATTTTTTGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCg                                                                                            >  2:2796436/1‑91 (MQ=255)
  agagATGCATTTTTTGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCgtt                                                                                          <  1:331258/91‑1 (MQ=255)
   gagaTGCATTTTTTGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTGTTACAATTAAAGAGCCAACCgttg                                                                                         >  2:3601272/1‑91 (MQ=255)
     gaTGCATTTTTTGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCgttgtt                                                                                       <  1:496083/91‑1 (MQ=255)
     gaTGCATTTTTTGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCgttgtt                                                                                       <  1:258132/91‑1 (MQ=255)
     gaTGCATTTTTTGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCgttgtt                                                                                       >  1:1782068/1‑91 (MQ=255)
      aTGCATTTTTTGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTta                                                                                      >  1:1885676/1‑91 (MQ=255)
       tGCATTTTTTGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTtat                                                                                     <  2:819262/91‑1 (MQ=255)
         cATTTTTTGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAg                                                                                   <  2:3205451/91‑1 (MQ=255)
         cATTTTTTGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAg                                                                                   >  1:2214790/1‑91 (MQ=255)
           ttttttGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTc                                                                                 >  1:3041715/1‑91 (MQ=255)
           ttttttGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTc                                                                                 >  1:3572726/1‑91 (MQ=255)
            tttttGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTGCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCt                                                                                <  2:1222784/91‑1 (MQ=255)
             ttttGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTa                                                                               >  2:3555242/1‑91 (MQ=255)
             ttttGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTa                                                                               >  2:3559973/1‑91 (MQ=255)
              tttGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTaa                                                                              <  1:1629617/91‑1 (MQ=255)
               ttGTGACTATGGTAGATCCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAATGAGACAACCGTTGTTATAGGCTAAc                                                                             <  2:2573539/91‑1 (MQ=18)
               ttGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTCTGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAAc                                                                             <  1:3868118/91‑1 (MQ=255)
               ttGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAAc                                                                             >  2:118452/1‑91 (MQ=255)
               ttGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAAc                                                                             >  1:2768557/1‑91 (MQ=255)
                tgtgACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACGATTAAAGAGCCAACCGTTGTTATAGGCTAACa                                                                            <  2:787522/91‑1 (MQ=255)
                  ttataaTGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAAt                                                                          >  2:3234142/5‑91 (MQ=255)
                  tgACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAGCAAt                                                                          <  2:3802961/91‑1 (MQ=255)
                  tgACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAAt                                                                          <  2:3086444/91‑1 (MQ=255)
                  tgACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAAt                                                                          >  1:2470904/1‑91 (MQ=255)
                      aaTGGTAGAACCTTTTCTGAACCATACACATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCATTGTTAGACTATAACAACGGtt                                                                      <  1:3724458/91‑1 (MQ=17)
                       aTGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCATTGTTAGACTATAACAACGGTTg                                                                     <  2:3311524/91‑1 (MQ=17)
                       aTGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCATTGTTAGACTATAACAACGGTTg                                                                     <  1:1461618/91‑1 (MQ=17)
                          gTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCATTGTTAGACTATAACAACGGTTGGct                                                                  <  1:2657276/91‑1 (MQ=18)
                           tAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCATTGTTAGACTATAACAACGGTTGGctc                                                                 <  1:3528568/91‑1 (MQ=18)
                             gAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGctcct                                                               >  2:1480174/1‑91 (MQ=255)
                             gAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGctcct                                                               <  2:3104641/91‑1 (MQ=255)
                             gAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGctcct                                                               <  1:784831/91‑1 (MQ=255)
                               aCCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCctct                                                             >  2:1775577/1‑91 (MQ=255)
                               aCCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCctct                                                             <  1:3457096/91‑1 (MQ=255)
                                      cTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTa                                                      >  2:236313/1‑91 (MQ=255)
                                      cTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCATTGTTAGACTATAACAACGGTTGGctc                                                                 >  2:2298368/1‑80 (MQ=18)
                                       tGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTAt                                                     >  2:2620107/1‑91 (MQ=255)
                                        gAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTAtg                                                    <  1:3203361/91‑1 (MQ=255)
                                        aaaCCATACGCTATGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTAtg                                                    >  1:1328775/2‑91 (MQ=39)
                                         aaCCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTAtgt                                                   >  1:547377/1‑91 (MQ=255)
                                         aaCCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTAtgt                                                   >  1:1718716/1‑91 (MQ=255)
                                          aCCATACGCTATGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTCAtgtg                                                  >  2:929472/1‑91 (MQ=255)
                                            cATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCt                                                <  2:2375397/91‑1 (MQ=255)
                                             aTACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTa                                               <  2:1782068/91‑1 (MQ=255)
                                             aTACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTa                                               <  2:391806/91‑1 (MQ=255)
                                               aCGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTaaa                                             >  1:2266564/1‑91 (MQ=255)
                                                  cATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAAt                                          >  2:2156627/1‑91 (MQ=255)
                                                  cATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAAt                                          <  2:2490321/91‑1 (MQ=255)
                                                   ctatgCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAAt                                          >  1:1622304/5‑91 (MQ=255)
                                                   ctatgCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAAt                                          >  2:1985160/5‑91 (MQ=255)
                                                   ctatgCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAAt                                          >  2:2185044/5‑91 (MQ=255)
                                                   ctatgCGTATGGCTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAAt                                          >  1:2982153/5‑91 (MQ=255)
                                                    tatgCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAAtt                                         >  1:1512919/4‑91 (MQ=255)
                                                    tatgCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAAtt                                         >  2:2017544/4‑91 (MQ=255)
                                                     atgCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAAttt                                        >  1:3056925/3‑91 (MQ=255)
                                                     aGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTACAAATTta                                       >  2:2741171/1‑91 (MQ=255)
                                                     aGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATCta                                       >  2:2741393/1‑91 (MQ=255)
                                                       gCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTtat                                      >  1:395860/1‑91 (MQ=255)
                                                          tATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAgg                                   <  2:3388483/91‑1 (MQ=255)
                                                           aTGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGc                                  >  1:1841279/1‑91 (MQ=255)
                                                           aTGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGc                                  >  2:2241660/1‑91 (MQ=255)
                                                           aTGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGACTCCTCTTATTTTATGTGCTAAAAATTTATAGGc                                  >  2:1335529/1‑91 (MQ=255)
                                                            tGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCa                                 <  2:1259178/91‑1 (MQ=255)
                                                             ggTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCaa                                >  2:2527895/1‑91 (MQ=255)
                                                             ggTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCaa                                >  1:1117142/1‑91 (MQ=255)
                                                             ggTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCaa                                <  1:3066516/91‑1 (MQ=255)
                                                               ttttCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAAtt                              >  1:1011396/1‑91 (MQ=255)
                                                                 ttCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAAtttt                            <  1:95646/91‑1 (MQ=255)
                                                                   cTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTttat                          <  2:1567076/91‑1 (MQ=255)
                                                                    tttttACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTttatt                         <  1:1167709/91‑1 (MQ=255)
                                                                    tttttACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTttatt                         <  2:1088412/91‑1 (MQ=255)
                                                                    tttttACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTttatt                         <  2:2698648/91‑1 (MQ=255)
                                                                    tttttACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTttatt                         <  2:2699120/91‑1 (MQ=255)
                                                                     ttttACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTttatta                        <  1:2063701/91‑1 (MQ=255)
                                                                     ttttACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTttatta                        <  2:3458466/91‑1 (MQ=255)
                                                                       ttACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTa                                               >  1:2304806/1‑66 (MQ=255)
                                                                       ttACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTaca                      <  1:1272306/91‑1 (MQ=255)
                                                                       ttACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTaca                      <  1:2084190/91‑1 (MQ=255)
                                                                        tACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTacaa                     >  1:670950/1‑91 (MQ=255)
                                                                         aCAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTacaac                    >  1:3123469/1‑91 (MQ=255)
                                                                          cAATTAAAGATCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTacaaca                   <  2:1242689/91‑1 (MQ=255)
                                                                          cAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTacaaca                   >  1:994462/1‑91 (MQ=255)
                                                                           aaTTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTAcaacaa                  <  2:1454750/91‑1 (MQ=255)
                                                                            aTTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTACAACAAt                 >  1:3633038/1‑91 (MQ=255)
                                                                             ttAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTACAACAATg                >  2:1796566/1‑91 (MQ=255)
                                                                             ttAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTACAACAATg                <  1:3859045/91‑1 (MQ=255)
                                                                              tAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTACAACAATGt               >  1:2609481/1‑91 (MQ=255)
                                                                              tAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTACAACAATGt               <  1:1771575/91‑1 (MQ=255)
                                                                               aaaGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAGAAATTTATAGGCAAATTTATTACAACAATGTa              >  1:3323341/1‑91 (MQ=255)
                                                                               aaaGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAAGTTTATTACAACAATGTa              >  1:3322455/1‑91 (MQ=255)
                                                                                    gCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTACAACAATGTACAttt         >  1:2844682/1‑91 (MQ=255)
                                                                                     ccAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTACAACAATGTACATTTa        <  2:1670200/91‑1 (MQ=255)
                                                                                       aaCCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTACAACAATGTACATTTAAg      <  2:855942/91‑1 (MQ=255)
                                                                                        aCCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTACAACAATGTACATTTAAgg     >  2:1727155/1‑91 (MQ=255)
                                                                                        aCCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTACAACAATGTACATTTAAgg     >  1:1240339/1‑91 (MQ=255)
                                                                                        aCCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTACAACAATGTACATTTAAgg     <  2:2834938/91‑1 (MQ=255)
                                                                                          cGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTACAACAATGTACATTTAAGGTg   <  1:2107036/91‑1 (MQ=255)
                                                                                           gttgttATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTACAACAATGTACATTTAAGGTGa  <  2:71426/91‑1 (MQ=255)
                                                                                           gttgttATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTACAACAATGTACATTTAAGGTGa  <  1:2000349/91‑1 (MQ=255)
                                                                                           |                                                                                          
GTAGAGATGCATTTTTTGTGACAATGGTAGAACCTTTTCTGAACCATACGCATAGCGTATGGTTTTCTTTTTACAATTAAAGAGCCAACCGTTGTTATAGTCTAACAATGGTTGGCTCCTCTTATTTTATGTGCTAAAAATTTATAGGCAATTTTATTACAACAATGTACATTTAAGGTGA  >  USA300TCH1516_ALE/1918465‑1918645

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