BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	462,826	T→A	6.2%	intergenic (+85/+18)	USA300HOU_RS02185 → / ← USA300HOU_RS02190	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	462,826	0	T	A	6.2%	450.6 / 15.8	130	intergenic (+85/+18)	USA300HOU_RS02185/USA300HOU_RS02190	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (52/70); new base A (6/2); total (58/72)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

Only 100 of 134 total aligned reads displayed.

TAATTAAATAAAACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTTTGATTGAA  >  USA300TCH1516_ALE/462739‑462916

taatTAAATAAAACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTAACa                                                                                         <  1:378474/91‑1 (MQ=255)

 aatTAAATAAAACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACaa                                                                                        <  1:205943/91‑1 (MQ=255)

   ttAAATAAAACAGCTTAATGATGTAATGAAATTTGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAg                                                                                      <  1:3391467/91‑1 (MQ=255)

    taaataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGc                                                                                     >  1:2065193/1‑91 (MQ=255)

    taaataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGc                                                                                     <  1:782804/91‑1 (MQ=255)

    taaataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGc                                                                                     <  1:2131160/91‑1 (MQ=255)

    taaataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGc                                                                                     >  1:2576307/1‑91 (MQ=255)

    taaataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGc                                                                                     >  1:2576564/1‑91 (MQ=255)

    taaataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGc                                                                                     >  1:3612924/1‑91 (MQ=255)

      aataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCat                                                                                   >  1:2503460/1‑91 (MQ=255)

       atataACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCatt                                                                                  <  1:3273860/91‑1 (MQ=255)

       atagaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAATCatt                                                                                  <  1:3273798/91‑1 (MQ=255)

       ataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCatt                                                                                  <  1:1214590/91‑1 (MQ=255)

        taaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCatta                                                                                 <  1:1027928/91‑1 (MQ=255)

        taaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCatta                                                                                 <  1:3597805/91‑1 (MQ=255)

         aaaaCAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCattat                                                                                <  1:4000473/91‑1 (MQ=255)

             cAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGaa                                                                            >  1:3945634/1‑91 (MQ=255)

             cAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGaa                                                                            >  1:1218691/1‑91 (MQ=255)

             cAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGaa                                                                            <  1:2232966/91‑1 (MQ=255)

             cAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGaa                                                                            <  1:1899277/91‑1 (MQ=255)

              aGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAat                                                                           >  1:1465718/1‑91 (MQ=255)

              aGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAat                                                                           >  1:705562/1‑91 (MQ=255)

               gCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAatt                                                                          >  1:1451746/1‑91 (MQ=255)

                cTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAattt                                                                         <  1:3061443/91‑1 (MQ=255)

                cTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAattt                                                                         <  1:494280/91‑1 (MQ=255)

                cTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAattt                                                                         <  1:2231102/91‑1 (MQ=255)

                cTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAattt                                                                         <  1:2230815/91‑1 (MQ=255)

                cTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAattt                                                                         <  1:818818/91‑1 (MQ=255)

                cTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAattt                                                                         >  1:1893969/1‑91 (MQ=255)

                  tAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAatttat                                                                       <  1:856323/91‑1 (MQ=255)

                   aaTGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAatttatt                                                                      <  1:2926092/91‑1 (MQ=255)

                   aaTGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAatttatt                                                                      <  1:2926091/91‑1 (MQ=255)

                   aaTGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAatttatt                                                                      >  1:280181/1‑91 (MQ=255)

                    atgatgTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAatttattt                                                                     <  1:797097/91‑1 (MQ=255)

                     tgatgTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTAtttt                                                                    >  1:1428674/1‑91 (MQ=255)

                     tgatgTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTAtttt                                                                    >  1:1428659/1‑91 (MQ=255)

                        tgTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACg                                                                 >  1:2075016/1‑91 (MQ=255)

                         gTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACgt                                                                >  1:2710600/1‑91 (MQ=255)

                            aTGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGtt                                                             <  1:1453951/91‑1 (MQ=255)

                            aTGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTGATTTTACGTGtt                                                             <  1:3204478/91‑1 (MQ=255)

                             tGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTc                                                            <  1:1241451/91‑1 (MQ=255)

                              gAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCa                                                           <  1:3040938/91‑1 (MQ=255)

                                aaTTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCata                                                         <  1:853828/91‑1 (MQ=255)

                                  ttAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATAtt                                                       <  1:3519011/91‑1 (MQ=255)

                                   tAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATAttt                                                      >  1:1744203/1‑91 (MQ=255)

                                      tGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGa                                                   >  1:1411206/1‑91 (MQ=255)

                                       gAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGaa                                                  >  1:323707/1‑91 (MQ=255)

                                       gAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGaa                                                  <  1:387991/91‑1 (MQ=255)

                                       gAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGaa                                                  >  1:1181391/1‑91 (MQ=255)

                                       gAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGaa                                                  <  1:3751230/91‑1 (MQ=255)

                                         gTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTATGCTTTGTTACAACGCATTATTGAATTTATTTTACGTGTTCATATTTCGAAAc                                                >  1:3852245/1‑91 (MQ=255)

                                            aaTCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATc                                             >  1:3073261/1‑91 (MQ=255)

                                              tCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCaa                                           >  1:461464/1‑91 (MQ=255)

                                               cACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCaaa                                          >  1:816190/1‑91 (MQ=255)

                                                 cTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGc                                        <  1:4071472/91‑1 (MQ=255)

                                                 cTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGc                                        <  1:134791/91‑1 (MQ=255)

                                                 cTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGc                                        >  1:443055/1‑91 (MQ=255)

                                                 cTGACTTCTACGTCATTGAGCTGTGTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGc                                        >  1:2755653/1‑91 (MQ=255)

                                                  tGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGcc                                       >  1:181072/1‑91 (MQ=255)

                                                  tGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGcc                                       >  1:181010/1‑91 (MQ=255)

                                                   gACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCg                                      >  1:1237005/1‑91 (MQ=255)

                                                     cTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTc                                    <  1:1212548/91‑1 (MQ=255)

                                                        cTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTg                                 <  1:131099/91‑1 (MQ=255)

                                                         tACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGt                                <  1:528226/91‑1 (MQ=255)

                                                          aCGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGtt                               <  1:3171677/91‑1 (MQ=255)

                                                          aCGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGtt                               <  1:1528637/91‑1 (MQ=255)

                                                           cGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGttt                              <  1:1980526/91‑1 (MQ=255)

                                                           cGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGttt                              <  1:1980587/91‑1 (MQ=255)

                                                           cGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGttt                              <  1:2060575/91‑1 (MQ=255)

                                                             tCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAg                            >  1:1508566/1‑91 (MQ=255)

                                                             tCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAg                            <  1:1119624/91‑1 (MQ=255)

                                                              cATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGc                           <  1:3771890/91‑1 (MQ=255)

                                                                   aGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTttgt                      <  1:1762248/91‑1 (MQ=255)

                                                                     cTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTttgttg                    <  1:1533361/91‑1 (MQ=255)

                                                                     cTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTttgttg                    <  1:2467004/91‑1 (MQ=255)

                                                                      tGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTGTGTTGa                   <  1:3223396/91‑1 (MQ=255)

                                                                      tGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTGTGTTGa                   <  1:3223397/91‑1 (MQ=255)

                                                                        ttttttGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATa                 <  1:190945/91‑1 (MQ=255)

                                                                        ttttttGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATa                 >  1:30897/1‑91 (MQ=255)

                                                                         ttttaGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATaa                >  1:3280652/1‑91 (MQ=255)

                                                                           tttGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATg              <  1:2973783/91‑1 (MQ=255)

                                                                           tttGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATg              <  1:2626634/91‑1 (MQ=255)

                                                                            ttGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAATGCCGTCTTGTTTAGCTTTGTTGATAATGt             >  1:259591/1‑91 (MQ=255)

                                                                            ttGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGt             >  1:210149/1‑91 (MQ=255)

                                                                            ttGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGt             <  1:1817799/91‑1 (MQ=255)

                                                                            ttGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGt             >  1:3827130/1‑91 (MQ=255)

                                                                            ttGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGt             <  1:1956311/91‑1 (MQ=255)

                                                                              gtgCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCt           <  1:4066087/91‑1 (MQ=255)

                                                                              gtgCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCt           >  1:699230/1‑91 (MQ=255)

                                                                               tgCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCtt          >  1:2724138/1‑91 (MQ=255)

                                                                                gCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCttt         >  1:2415940/1‑91 (MQ=255)

                                                                                 cTTTGTAACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTttg        >  1:3798561/1‑91 (MQ=255)

                                                                                   ttGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTttgat      >  1:1393893/1‑91 (MQ=255)

                                                                                   ttGTAACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTttgat      >  1:1013554/1‑91 (MQ=255)

                                                                                   ttGTAACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTttgat      >  1:3905388/1‑91 (MQ=255)

                                                                                   ttGTAACAAAGCATTATTGAATTTATTTTACGTGTTCATATCTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTttgat      >  1:3596452/1‑91 (MQ=255)

                                                                                     gTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTttgattg    <  1:1710138/91‑1 (MQ=255)

                                                                                     gTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTttgattg    >  1:3654289/1‑91 (MQ=255)

                                                                                       tACAATGCATTATTGAATTTATTTTTCGTGTTCATATTTGGAAACATCAAAGCTGTCTTGTTTAGCTTTGTTGATAATGTCTTTGATTGaa  <  1:3856837/91‑1 (MQ=255)

                                                                                       tACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTTTGATTGaa  <  1:2791997/91‑1 (MQ=255)

TAATTAAATAAAACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTTTGATTGAA  >  USA300TCH1516_ALE/462739‑462916

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