Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,219,579 A→T 6.2% V297V (GTT→GTA USA300HOU_RS06070 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,219,5790AT 6.2% 696.0 / 23.4 209V297V (GTT→GTAUSA300HOU_RS06070hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (103/93);  new base T (7/6);  total (110/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 214 total aligned reads displayed.

ATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTAATGAA  >  USA300TCH1516_ALE/1219489‑1219669
                                                                                          |                                                                                          
atatAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGa                                                                                            >  2:1614319/1‑91 (MQ=255)
 tataACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGaa                                                                                           <  2:1247087/91‑1 (MQ=255)
 tataACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGaa                                                                                           <  1:1874616/91‑1 (MQ=255)
  ataACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAAc                                                                                          <  1:3565396/91‑1 (MQ=255)
  ataACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAAc                                                                                          <  1:1276517/91‑1 (MQ=255)
  ataACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAAc                                                                                          <  2:3799783/91‑1 (MQ=255)
   taACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACa                                                                                         <  1:3484712/91‑1 (MQ=255)
   taACTGTTCAGTATCTTTATTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTGTACaa                                                                                        <  2:3014097/91‑1 (MQ=37)
    aaCTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACaa                                                                                        >  1:2509596/1‑91 (MQ=255)
     aCTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTGTACaaa                                                                                       <  1:306551/91‑1 (MQ=38)
     aCTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTGTACaaa                                                                                       <  1:307002/91‑1 (MQ=38)
     aCTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACaaa                                                                                       >  1:2716200/1‑91 (MQ=255)
      cTGTTCAGTATCTTTATTTTTTTTTTTTTTATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAAt                                                                                      >  2:2157067/1‑91 (MQ=255)
         ttCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGa                                                                                   >  1:1162286/1‑91 (MQ=255)
         ttCAGGATCTGTATTTTGAGACTGCTCATATTCTTCAAGCAACTTCGCTAATTTATTTTGATATGTGTACAACGGCTGTTGAACAAATCGa                                                                                   <  1:3718486/91‑1 (MQ=255)
          tCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGaa                                                                                  <  2:1579223/91‑1 (MQ=255)
            aGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACt                                                                                <  1:1368493/91‑1 (MQ=255)
             gTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTa                                                                               <  2:2417007/91‑1 (MQ=255)
                tCTTTATTTCTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAAt                                                                            >  1:1946357/1‑91 (MQ=255)
                 cTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATc                                                                           <  1:1095161/91‑1 (MQ=255)
                   ttATTTTTAGACTGCTCATTTTCTTCAATCAACTTCGCTAGTTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCAt                                                                         <  2:663657/91‑1 (MQ=255)
                     aTTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTGCAACTGCTGTTGAACAAATCGAACTAAATCATTc                                                                       >  1:529543/1‑91 (MQ=255)
                     aTTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTc                                                                       >  2:2858432/1‑91 (MQ=255)
                      tttttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCg                                                                      <  2:2061402/91‑1 (MQ=255)
                      tttttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCg                                                                      <  2:842500/91‑1 (MQ=255)
                       ttttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGc                                                                     >  2:1453356/1‑91 (MQ=255)
                       ttttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGc                                                                     >  1:1390064/1‑91 (MQ=255)
                        tttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCa                                                                    >  1:3797933/1‑91 (MQ=255)
                        tttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCa                                                                    <  1:3466227/91‑1 (MQ=255)
                          tAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACg                                                                  <  2:1661031/91‑1 (MQ=255)
                          tAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACg                                                                  <  2:2313676/91‑1 (MQ=255)
                           aGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgt                                                                 >  2:12908/1‑91 (MQ=255)
                           aGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgt                                                                 >  2:1039290/1‑91 (MQ=255)
                           aGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgt                                                                 >  1:3345694/1‑91 (MQ=255)
                            gACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgtt                                                                <  1:1159561/91‑1 (MQ=255)
                               tGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgttgtt                                                             <  1:3085886/91‑1 (MQ=255)
                                gCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTTAACAAATCGAACTAAATCATTCGCACGTTGttt                                                            >  2:1338072/1‑91 (MQ=255)
                                 cTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTa                                                           >  2:3277879/1‑91 (MQ=255)
                                      atTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAAcgcg                                                      >  1:3207935/1‑91 (MQ=255)
                                      atTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAAcgcg                                                      <  2:927852/91‑1 (MQ=255)
                                      atTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAAcgcg                                                      >  1:3858709/1‑91 (MQ=255)
                                       ttcttcAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGt                                                     <  2:2669343/91‑1 (MQ=255)
                                       ttcttcAATCAACTTCGCTAATTTATCTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGt                                                     >  2:2211878/1‑91 (MQ=255)
                                        tcttcAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGtt                                                    >  1:447437/1‑91 (MQ=255)
                                        tcttcAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGtt                                                    >  1:257354/1‑91 (MQ=255)
                                         cttcAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTc                                                   <  1:952231/91‑1 (MQ=255)
                                          ttcAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCa                                                  >  2:2748567/1‑91 (MQ=255)
                                             aatcaaCTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGt                                               >  2:1952122/1‑91 (MQ=255)
                                             aatcaaCTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGt                                               >  2:137715/1‑91 (MQ=255)
                                               tcaaCTTCGCTAATTTATTTTGATATTTGTGCAGCTGCTGTTGAACAAATCGAACTAAATCATTCGCATGTTGTTTAACGCGTTCACGTTc                                             <  2:3425040/91‑1 (MQ=255)
                                                  aCTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGcc                                          <  1:1330222/91‑1 (MQ=255)
                                                   cTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGccc                                         >  1:875797/1‑90 (MQ=255)
                                                   cTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCa                                         <  2:897581/91‑1 (MQ=255)
                                                   cTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCa                                         <  2:369259/91‑1 (MQ=255)
                                                   cTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGAGTTCACGTTCGTCa                                         >  1:77718/1‑91 (MQ=255)
                                                     tCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcg                                       >  2:257857/1‑91 (MQ=255)
                                                      cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgc                                      >  2:1098357/1‑91 (MQ=255)
                                                      cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgc                                      >  1:1524520/1‑91 (MQ=255)
                                                      cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgc                                      >  1:2595480/1‑91 (MQ=255)
                                                      cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgc                                      <  1:2139045/91‑1 (MQ=255)
                                                       gCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgcg                                     >  2:448965/1‑91 (MQ=255)
                                                        cTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACgcgc                                    >  2:3480844/1‑91 (MQ=255)
                                                        cTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACgcgc                                    >  2:3558455/1‑91 (MQ=255)
                                                         tAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGca                                   >  2:3622855/1‑91 (MQ=255)
                                                          aaTTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcat                                  <  1:2090486/91‑1 (MQ=255)
                                                           atttatttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcatc                                 >  2:1430066/1‑91 (MQ=255)
                                                           atttatttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcatc                                 >  1:3145926/1‑91 (MQ=255)
                                                                ttttGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaaa                            >  2:3531033/1‑91 (MQ=255)
                                                                ttttGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaaa                            <  1:2587173/91‑1 (MQ=255)
                                                                ttttGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaaa                            <  1:1614319/91‑1 (MQ=255)
                                                                 tttGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaaaa                           >  2:402569/1‑91 (MQ=255)
                                                                  ttGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaaaaa                          >  2:3746888/1‑91 (MQ=255)
                                                                    gATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCCTCATAAAAACg                        >  2:172881/1‑91 (MQ=255)
                                                                        tttGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCa                    <  2:3170058/91‑1 (MQ=255)
                                                                         ttGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCaa                   <  1:225879/91‑1 (MQ=255)
                                                                         ttGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCaa                   >  2:304365/1‑91 (MQ=255)
                                                                          tgttcAACAGCAGTTGTACAAATCGAACTAAATCATTCGCACGGTGTTTACCGCGTTCACGTTCGCCACGCGCATCCTAAAAACGCTCAAg                  >  2:1596980/5‑91 (MQ=255)
                                                                          tGTTCAACAGCAGTTGTACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAg                  >  2:714499/1‑91 (MQ=255)
                                                                           gttcAACAGCAGTTGTACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGc                 >  1:486446/4‑91 (MQ=255)
                                                                           gTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGc                 >  1:2809030/1‑91 (MQ=255)
                                                                           gTACAACTGCTGTTGAACAAATCGAACTAAACCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCATGCGCATCATAATAACGATCAAGc                 >  2:2966118/1‑91 (MQ=255)
                                                                            tACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCa                <  2:2127670/91‑1 (MQ=255)
                                                                             tcAACAGCAGTTGTACAAATCGAACTAAATCATTCGCCCGTTGTTTAACGCGTTCACGTTCGCCACGCGCCTCATAAAAACGCCCAAGCaa               >  2:54808/2‑91 (MQ=255)
                                                                              cAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCaaa              <  1:2191554/91‑1 (MQ=255)
                                                                              cAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCaaa              <  2:3590057/91‑1 (MQ=255)
                                                                                aCTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATc            >  2:1743893/1‑91 (MQ=255)
                                                                                aCTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATc            >  1:2642532/1‑91 (MQ=255)
                                                                                 ctgctgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCa           >  2:2778558/1‑91 (MQ=255)
                                                                                 ctgctgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCa           >  1:2405420/1‑91 (MQ=255)
                                                                                    cagTTGTACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCAttt        >  2:153583/3‑91 (MQ=255)
                                                                                     tgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTa       <  1:1502007/91‑1 (MQ=255)
                                                                                     tgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTa       <  1:2281175/91‑1 (MQ=255)
                                                                                     tgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTa       <  1:768377/91‑1 (MQ=255)
                                                                                     tgTTGAACAAATCGAACTAAATCATACGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTa       <  1:3607050/91‑1 (MQ=255)
                                                                                      gTTGTACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTaa      >  1:1833994/1‑91 (MQ=255)
                                                                                      gTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTaa      <  1:1450484/91‑1 (MQ=255)
                                                                                      gTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTaa      <  1:3248466/91‑1 (MQ=255)
                                                                                        tGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTAATg    >  2:2898429/1‑91 (MQ=255)
                                                                                         gAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTAATGa   >  1:1003828/1‑91 (MQ=255)
                                                                                          aaCAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTAATGaa  >  1:2348912/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTAATGAA  >  USA300TCH1516_ALE/1219489‑1219669

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