BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,513,966	A→G	5.1%	L57S (TTG→TCG)	USA300HOU_RS12940 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,513,966	0	A	G	5.1%	826.4 / 23.4	257	L57S (TTG→TCG)	USA300HOU_RS12940	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (130/113); new base G (10/3); total (141/116)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 260 total aligned reads displayed.

AGCTAAACCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCAA  >  USA300TCH1516_ALE/2513881‑2514056

aGCTAAACCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTa                                                                                       >  1:3736541/1‑91 (MQ=255)

aGCTAAACCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTa                                                                                       >  1:90880/1‑91 (MQ=255)

 gCTAAACCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTaa                                                                                      >  2:2921589/1‑91 (MQ=255)

  cTAAACCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTTCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTaaa                                                                                     >  1:2529828/1‑91 (MQ=255)

  cTAAACCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTaaa                                                                                     >  1:1669780/1‑91 (MQ=255)

  cTAAACCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACACTTCCAGCTAATCCACGACCTTTTCTTAAACTAATCAATTTaaa                                                                                     >  1:3876917/1‑91 (MQ=255)

   tAAACCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAAt                                                                                    >  1:3492963/1‑91 (MQ=255)

   tAAACCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAAt                                                                                    <  1:1905348/91‑1 (MQ=255)

    aaaCCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATc                                                                                   <  2:1155168/91‑1 (MQ=255)

    aaaCCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATc                                                                                   <  2:3239169/91‑1 (MQ=255)

     aaCCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCg                                                                                  >  2:1001590/1‑91 (MQ=255)

     aaCCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCg                                                                                  <  2:3097320/91‑1 (MQ=255)

       ccAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCGATTTAAATTGat                                                                                <  2:1213637/91‑1 (MQ=255)

       ccAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGat                                                                                >  1:1340755/1‑91 (MQ=255)

        cAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGatt                                                                               >  2:669368/1‑91 (MQ=255)

          aCATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGattat                                                                             >  2:2338045/1‑91 (MQ=255)

          aCATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGattat                                                                             >  2:112798/1‑91 (MQ=255)

          aCATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGattat                                                                             >  1:2424101/1‑91 (MQ=255)

              tGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTAttttt                                                                         >  1:1589867/1‑91 (MQ=255)

               ggCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTtta                                                                        >  2:3402897/1‑91 (MQ=255)

                gCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTttat                                                                       >  2:3379662/1‑91 (MQ=255)

                    ttACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTAcc                                                                   >  2:2626201/1‑91 (MQ=255)

                      aCCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTg                                                                 <  2:2286920/91‑1 (MQ=255)

                         aTACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGaaa                                                              <  1:808860/91‑1 (MQ=255)

                         aTACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGaaa                                                              <  2:3288285/91‑1 (MQ=255)

                         aTACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGaaa                                                              >  1:3176803/1‑91 (MQ=255)

                          tACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAAc                                                             <  2:1762129/91‑1 (MQ=255)

                          tACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAAc                                                             <  1:2324904/91‑1 (MQ=255)

                           aCGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACa                                                            >  2:2008187/1‑91 (MQ=255)

                            cGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACAt                                                           <  1:950665/91‑1 (MQ=255)

                               tttCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGt                                                        <  1:403054/91‑1 (MQ=255)

                                ttCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTg                                                       <  1:1299271/91‑1 (MQ=255)

                                 tCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGc                                                      >  1:2599328/1‑91 (MQ=255)

                                 tCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGc                                                      <  1:848604/91‑1 (MQ=255)

                                 tCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGc                                                      <  2:46747/91‑1 (MQ=255)

                                 tCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGc                                                      <  1:1125297/91‑1 (MQ=255)

                                  ccTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGcc                                                     >  1:1051601/1‑91 (MQ=255)

                                  ccTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGcc                                                     >  2:1682952/1‑91 (MQ=255)

                                  ccTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGcc                                                     <  2:1948011/91‑1 (MQ=255)

                                  ccTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGcc                                                     <  1:3771636/91‑1 (MQ=255)

                                  ccTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGcc                                                     <  1:1017613/91‑1 (MQ=255)

                                   cTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCa                                                    <  2:861136/91‑1 (MQ=255)

                                   cTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCa                                                    >  1:2773486/1‑91 (MQ=255)

                                        ttCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTa                                               >  2:2072785/1‑91 (MQ=255)

                                          cATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTAtt                                             >  2:2609351/1‑91 (MQ=255)

                                            tAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTgg                                           <  1:1482098/91‑1 (MQ=255)

                                             aaCAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGga                                          >  1:537062/1‑91 (MQ=255)

                                             aaCAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGga                                          >  2:1409581/1‑91 (MQ=255)

                                             aaCAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGga                                          >  1:1885716/1‑91 (MQ=255)

                                             aaCAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGga                                          <  2:1099762/91‑1 (MQ=255)

                                               cAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGgaga                                        <  2:1819863/91‑1 (MQ=255)

                                               cAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGgaga                                        <  1:2964443/91‑1 (MQ=255)

                                                aGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgat                                       <  2:135190/91‑1 (MQ=255)

                                                 gTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgatg                                      <  2:266367/91‑1 (MQ=255)

                                                  ttCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgatga                                     >  1:3027909/1‑91 (MQ=255)

                                                  ttCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgatga                                     <  2:1471182/91‑1 (MQ=255)

                                                   tCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgatgat                                    <  2:3650864/91‑1 (MQ=255)

                                                   tCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgatgat                                    >  1:2551551/1‑91 (MQ=255)

                                                    ccAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGatgatg                                   >  1:2108054/1‑91 (MQ=255)

                                                    ccAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGatgatg                                   >  1:2107729/1‑91 (MQ=255)

                                                    ccAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGatgatg                                   >  1:3055310/1‑91 (MQ=255)

                                                    ccAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGatgatg                                   >  2:528018/1‑91 (MQ=255)

                                                    ccAGCTAATCCACGACCCTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACTTGAAACATAGTGCCATTTTATTGGAGatgatg                                   >  1:337791/1‑91 (MQ=255)

                                                     cAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATgg                                  >  2:2833833/1‑91 (MQ=255)

                                                     cAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATgg                                  <  2:2270659/91‑1 (MQ=255)

                                                      aGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGt                                 <  1:3256887/91‑1 (MQ=255)

                                                       gCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGtt                                <  2:2168631/91‑1 (MQ=255)

                                                          aaTCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGt                             <  2:934589/91‑1 (MQ=255)

                                                           aTCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGtt                            >  2:2568696/1‑91 (MQ=255)

                                                           aTCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGtt                            >  2:1094147/1‑91 (MQ=255)

                                                           aTCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGtt                            <  2:3877708/91‑1 (MQ=255)

                                                             ccACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAg                          >  2:3513952/1‑91 (MQ=255)

                                                             ccACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAg                          >  2:1631700/1‑91 (MQ=255)

                                                              cACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGa                         >  1:3342611/1‑91 (MQ=255)

                                                                cGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGAtt                       <  2:1917822/91‑1 (MQ=255)

                                                                cGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCCGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGAtt                       >  2:81835/1‑91 (MQ=255)

                                                                 gACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTc                      >  2:3478258/1‑91 (MQ=255)

                                                                  aCCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCa                     <  1:671810/91‑1 (MQ=255)

                                                                   ccTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCAt                    <  1:668111/91‑1 (MQ=255)

                                                                     ttttCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAg                  >  2:241865/1‑91 (MQ=255)

                                                                       ttCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGaa                >  2:3517632/1‑91 (MQ=255)

                                                                        tCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGaaa               <  2:1441708/91‑1 (MQ=255)

                                                                        tCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGaaa               <  2:1219607/91‑1 (MQ=255)

                                                                        tCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGaaa               <  2:2012140/91‑1 (MQ=255)

                                                                        tCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGaaa               <  2:2158729/91‑1 (MQ=255)

                                                                         cTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAg              >  1:3105064/1‑91 (MQ=255)

                                                                          ttAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGc             <  1:230737/91‑1 (MQ=255)

                                                                           tAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCg            <  2:3312644/91‑1 (MQ=255)

                                                                            aaaaTAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGa           <  2:2152177/91‑1 (MQ=255)

                                                                             aaaTAATCGATTTAAATTGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGAt          >  2:1297586/1‑91 (MQ=255)

                                                                             aaaTAATCGATTTAAATTGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGAt          >  2:428566/1‑91 (MQ=255)

                                                                             aaaTAATCGATTTAAATTGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGAt          >  1:2221633/1‑91 (MQ=255)

                                                                             aaaTAATCGATTTAAATTGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGAt          >  2:697105/1‑91 (MQ=255)

                                                                               ataatCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTg        <  1:1258429/91‑1 (MQ=255)

                                                                               ataatCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTg        >  1:3509747/1‑91 (MQ=255)

                                                                                taatCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGc       >  1:3051698/1‑91 (MQ=255)

                                                                                 aatCAATTTAAATCGTTTATTTTTAGCACCTGAAACATAGGGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGcc      <  2:3357217/91‑1 (MQ=255)

                                                                                   tCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGc    >  1:2464326/1‑91 (MQ=255)

                                                                                    cAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCa   >  1:308880/1‑91 (MQ=255)

                                                                                    cAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCa   <  1:2934965/91‑1 (MQ=255)

AGCTAAACCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCAA  >  USA300TCH1516_ALE/2513881‑2514056

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: `‑`