Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 665,392 T→G L79* (TTA→TGA)  paiA_1 → Spermidine/spermine N(1)‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE665,3920TG99.0% 379.4 / ‑3.6 100L79* (TTA→TGA) paiA_1Spermidine/spermine N(1)‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (57/42);  total (57/43)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 101 total aligned reads displayed.

AGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGA  >  USA300TCH1516_ALE/665305‑665479
                                                                                       |                                                                                       
agaaaagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATgata                                                                                      <  1:1573710/91‑1 (MQ=255)
 gaaaagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGatat                                                                                     >  1:2039452/1‑91 (MQ=255)
  aaaagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATAtt                                                                                    <  1:654681/91‑1 (MQ=255)
   aaagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATAttt                                                                                   >  1:3663536/1‑91 (MQ=255)
   aaagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATAttt                                                                                   >  1:2280253/1‑91 (MQ=255)
    aagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTa                                                                                  >  1:2339995/1‑91 (MQ=255)
     agaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTAt                                                                                 >  1:631071/1‑91 (MQ=255)
      gaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATc                                                                                <  1:2098523/91‑1 (MQ=255)
       aCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCa                                                                               <  1:934902/91‑1 (MQ=255)
       aCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCa                                                                               <  1:1775256/91‑1 (MQ=255)
       aCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCa                                                                               <  1:2763727/91‑1 (MQ=255)
       aCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCa                                                                               <  1:3054880/91‑1 (MQ=255)
        cTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAg                                                                              >  1:1523648/1‑91 (MQ=255)
         ttCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGc                                                                             >  1:1009452/1‑91 (MQ=255)
         ttCAAGAACAAGTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGc                                                                             <  1:3356044/91‑1 (MQ=255)
          tCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCt                                                                            <  1:385981/91‑1 (MQ=255)
          tCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCt                                                                            >  1:2755460/1‑91 (MQ=255)
           cAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTc                                                                           >  1:1113997/1‑91 (MQ=255)
             aGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCAt                                                                         <  1:918139/91‑1 (MQ=255)
             aGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCAt                                                                         <  1:656518/91‑1 (MQ=255)
             aGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCAt                                                                         <  1:656410/91‑1 (MQ=255)
             aGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGATGAAGAATGATATTTATCAGCTCAt                                                                         >  1:3433599/1‑91 (MQ=255)
               aaCAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTa                                                                       >  1:1270405/1‑91 (MQ=255)
                ataaTTGTTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATAAt                                                                      >  1:387663/3‑91 (MQ=255)
                 taaTTGTTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATg                                                                     >  1:1137151/2‑91 (MQ=255)
                 cAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATg                                                                     >  1:1809334/1‑91 (MQ=255)
                  aaTTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGt                                                                    >  1:2481491/1‑91 (MQ=255)
                  aaTTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGt                                                                    >  1:1587321/1‑91 (MQ=255)
                    ttattCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTa                                                                  <  1:3118346/91‑1 (MQ=255)
                      attCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTaaa                                                                <  1:1543394/91‑1 (MQ=255)
                      attCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTaaa                                                                <  1:3259266/91‑1 (MQ=255)
                      attCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTaaa                                                                <  1:3536380/91‑1 (MQ=255)
                       ttCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAAc                                                               <  1:424201/91‑1 (MQ=255)
                       ttCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAAc                                                               >  1:2376583/1‑91 (MQ=255)
                       ttCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAAc                                                               <  1:65197/91‑1 (MQ=255)
                        tCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAAcc                                                              >  1:198675/1‑91 (MQ=255)
                        tCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAAcc                                                              >  1:2834376/1‑91 (MQ=255)
                        tCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAAcc                                                              <  1:1605208/91‑1 (MQ=255)
                        tCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAAcc                                                              >  1:198584/1‑91 (MQ=255)
                          ttAGTCGTTAAAGAAGGTAATGACATCTTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAg                                                            <  1:1308215/91‑1 (MQ=255)
                           tAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGa                                                           >  1:3499039/1‑91 (MQ=255)
                            aGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGaa                                                          <  1:821950/91‑1 (MQ=255)
                            aGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGaa                                                          >  1:318288/1‑91 (MQ=255)
                             gTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAAt                                                         >  1:2286461/1‑91 (MQ=255)
                             gTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAAt                                                         <  1:561027/91‑1 (MQ=255)
                                  tGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCaa                                                    <  1:563119/91‑1 (MQ=255)
                                    aagaagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACTAGAATCGCAACa                                                  <  1:2837967/91‑1 (MQ=255)
                                    aagaagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACa                                                  <  1:3144965/91‑1 (MQ=255)
                                      gaagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATa                                                >  1:1472822/1‑91 (MQ=255)
                                       aagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATAc                                               >  1:1351094/1‑91 (MQ=255)
                                       aagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATAc                                               <  1:3449425/91‑1 (MQ=255)
                                       aagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATAc                                               >  1:825616/1‑91 (MQ=255)
                                         ggTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAg                                             <  1:1049704/91‑1 (MQ=255)
                                          gTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAgg                                            <  1:3389181/91‑1 (MQ=255)
                                          gTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAgg                                            <  1:3389173/91‑1 (MQ=255)
                                           tAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGt                                           >  1:276367/1‑91 (MQ=255)
                                           tAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGt                                           >  1:2038072/1‑91 (MQ=255)
                                            aaTGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGtt                                          <  1:2083666/91‑1 (MQ=255)
                                              tGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTAt                                        >  1:1488366/1‑91 (MQ=255)
                                              tGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTAt                                        >  1:3342736/1‑91 (MQ=255)
                                              tGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTAt                                        >  1:3342735/1‑91 (MQ=255)
                                                  aTCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTa                                    <  1:3380065/91‑1 (MQ=255)
                                                    cGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACa                                  >  1:2802257/1‑91 (MQ=255)
                                                    cGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACa                                  >  1:1992445/1‑91 (MQ=255)
                                                     gTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAg                                 >  1:2822184/1‑91 (MQ=255)
                                                     gTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAg                                 >  1:918733/1‑91 (MQ=255)
                                                       tGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCa                               >  1:159053/1‑91 (MQ=255)
                                                       tGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCa                               >  1:457626/1‑91 (MQ=255)
                                                        ggCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCTTTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCAt                              >  1:2932817/1‑91 (MQ=255)
                                                        ggCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCAt                              >  1:2931428/1‑91 (MQ=255)
                                                        ggCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCAt                              >  1:1667796/1‑91 (MQ=255)
                                                        ggCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCAt                              >  1:1667866/1‑91 (MQ=255)
                                                        ggCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCAt                              >  1:1710263/1‑91 (MQ=255)
                                                           tttGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGt                           >  1:2520944/1‑91 (MQ=255)
                                                           tttGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGt                           >  1:750788/1‑91 (MQ=255)
                                                           tttGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGt                           >  1:3490875/1‑91 (MQ=255)
                                                            ttGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCGCATTATGTTAAACCAGAAGCGCAACATACAGGTTATGGTACAGCATTGtt                          <  1:1592942/91‑1 (MQ=255)
                                                              gCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTaa                        >  1:2365943/1‑91 (MQ=255)
                                                               cTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTaaa                       <  1:2301964/91‑1 (MQ=255)
                                                                  aCTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGa                    >  1:1752547/1‑91 (MQ=255)
                                                                  aCTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGa                    <  1:676456/91‑1 (MQ=255)
                                                                  aCTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGa                    <  1:1074881/91‑1 (MQ=255)
                                                                   cTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGaa                   >  1:3496130/1‑91 (MQ=255)
                                                                   cTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGaa                   >  1:1181943/1‑91 (MQ=255)
                                                                   cTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGaa                   >  1:1182056/1‑91 (MQ=255)
                                                                   cTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGaa                   <  1:1614733/91‑1 (MQ=255)
                                                                    tttatttaCGGTGAAGAATGATATTTATCAGCTCATTATGTTAATCCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAg                  <  1:2939746/91‑1 (MQ=255)
                                                                    tttatttaCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAg                  >  1:1429308/1‑91 (MQ=255)
                                                                    tttatttaCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAg                  >  1:1429315/1‑91 (MQ=255)
                                                                          taCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTAt            <  1:1575143/91‑1 (MQ=255)
                                                                           aCGGTGGAGGATGATATTTATCTGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATc           <  1:3496616/91‑1 (MQ=255)
                                                                           aCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATc           <  1:2060736/91‑1 (MQ=255)
                                                                           aCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATc           >  1:2357144/1‑91 (MQ=255)
                                                                           aCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATc           <  1:2222065/91‑1 (MQ=255)
                                                                           aCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATc           <  1:1169456/91‑1 (MQ=255)
                                                                            cGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCa          <  1:1195105/91‑1 (MQ=255)
                                                                            cGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCa          >  1:628627/1‑91 (MQ=255)
                                                                               tGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGt       >  1:1192064/1‑91 (MQ=255)
                                                                                  agaaTGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCACTGTTAAATGAAGGATTATCACGtttt    >  1:3001988/1‑91 (MQ=255)
                                                                                    aaTGATATTTATCAGCTCATTATGTTAAACCAGAATCACAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGa  >  1:3127234/1‑91 (MQ=255)
                                                                                       |                                                                                       
AGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGA  >  USA300TCH1516_ALE/665305‑665479

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: