BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BX571857	1,167,297	A→G	5.2%	K71R (AAA→AGA)	SAS1093 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BX571857	1,167,297	0	A	G	5.2%	696.1 / 14.3	186	K71R (AAA→AGA)	SAS1093	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (95/79); new base G (2/9); total (97/89)
Fisher's exact test for biased strand distribution p-value = 2.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

GTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATACTTT  >  BX571857/1167200‑1167394

gtggggCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCcaacaa                                                                                                   <  1:920107/98‑1 (MQ=255)

 tggggCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCaacaaa                                                                                                  >  1:137672/1‑98 (MQ=255)

  ggggCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCaacaaac                                                                                                 >  2:653335/1‑98 (MQ=255)

   gggCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCaacaaaca                                                                                                >  2:1699287/1‑98 (MQ=255)

   gggCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCaacaaaca                                                                                                <  2:54603/98‑1 (MQ=255)

     acccAACACAGAAGCTGGCGAAAAGACAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAAt                                                                                              <  1:2787815/97‑1 (MQ=255)

         aaCCCAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCCACAGAGAGAAATAGGATCACCAATTCCCACAAACAATGCaa                                                                                          >  2:3572501/1‑98 (MQ=255)

          acacAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAg                                                                                         <  2:2401259/98‑1 (MQ=255)

           cacaGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGt                                                                                        >  2:2290948/1‑98 (MQ=255)

           cacaGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGt                                                                                        >  1:1065900/1‑98 (MQ=255)

            acaGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGtt                                                                                       >  2:3233201/1‑98 (MQ=255)

             caGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTg                                                                                      <  1:1817337/98‑1 (MQ=255)

             caGAAACTGGCGAAAAGCCAGCATACAAGAATGTGCAAGTTGGCGGGGCCCGAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTg                                                                                      <  2:2798068/98‑1 (MQ=21)

              aGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTgg                                                                                     >  2:2198762/1‑98 (MQ=255)

              aGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTgg                                                                                     <  1:1714747/98‑1 (MQ=255)

                  gCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCggg                                                                                 >  1:375972/1‑98 (MQ=255)

                    tGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGc                                                                               <  1:1323434/98‑1 (MQ=255)

                    tGGCGAAAAGCCAGCAAACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGc                                                                               <  2:2895502/98‑1 (MQ=255)

                     ggCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGcc                                                                              >  2:2856923/1‑98 (MQ=255)

                       cGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGcccc                                                                            >  1:3505357/1‑98 (MQ=255)

                        gAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCa                                                                           >  2:1112076/1‑98 (MQ=255)

                          aaaGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAAc                                                                         >  1:1115186/1‑98 (MQ=255)

                          aaaGCCAGCAAACAAAAATGTGAAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAAc                                                                         <  2:292527/98‑1 (MQ=255)

                            aGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACAt                                                                       >  1:356802/1‑98 (MQ=255)

                             gTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATa                                                                      >  2:2069873/1‑98 (MQ=255)

                             gTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATa                                                                      >  1:2988538/1‑98 (MQ=255)

                              tCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATag                                                                     <  1:2607031/98‑1 (MQ=255)

                              tCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATag                                                                     <  1:2862709/98‑1 (MQ=255)

                              tCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATag                                                                     <  1:3539995/98‑1 (MQ=255)

                               cAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATaga                                                                    >  1:1146413/1‑98 (MQ=255)

                                aGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATagag                                                                   >  1:1670687/1‑98 (MQ=255)

                                 gCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAgaga                                                                  >  2:1598058/1‑98 (MQ=255)

                                  cATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGaa                                                                 >  2:3491889/1‑98 (MQ=255)

                                   aaaCAGAAATGTGAAAGCTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGaaa                                                                <  2:3625478/96‑1 (MQ=17)

                                   aTACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGaaa                                                                >  1:2185400/1‑98 (MQ=255)

                                   aTACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGaaa                                                                >  2:407218/1‑98 (MQ=255)

                                     aCAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAAtt                                                              >  1:2909270/1‑98 (MQ=255)

                                     aCAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAAtt                                                              >  2:1608103/1‑98 (MQ=255)

                                     aCAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAAtt                                                              >  2:1160236/1‑98 (MQ=255)

                                       aaaaaTGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTgg                                                            >  1:2249827/1‑98 (MQ=255)

                                       aaaaaTGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTgg                                                            >  2:1531415/1‑98 (MQ=255)

                                         aaaTGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGt                                                          >  1:3593883/1‑98 (MQ=255)

                                          caTGGGCCAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAGCAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTc                                                         <  2:1906498/97‑1 (MQ=17)

                                          aaTGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTc                                                         <  1:881793/98‑1 (MQ=255)

                                               ccAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCaa                                                    <  1:1937720/97‑1 (MQ=255)

                                                  aGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAAttt                                                 >  1:405715/1‑98 (MQ=255)

                                                   gTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTc                                                >  1:2699890/1‑98 (MQ=255)

                                                    ttGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCa                                               >  2:873937/1‑98 (MQ=255)

                                                    ttGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCa                                               >  2:2343221/1‑98 (MQ=255)

                                                           ggCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACa                                        <  2:2654658/98‑1 (MQ=255)

                                                           ggCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACa                                        <  2:2145045/98‑1 (MQ=255)

                                                           ggCCACAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACa                                        <  1:1389640/98‑1 (MQ=21)

                                                           ggCCACAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACa                                        <  1:2500459/98‑1 (MQ=17)

                                                            gCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACaa                                       <  2:3389085/98‑1 (MQ=255)

                                                            gCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACaa                                       >  1:1762331/1‑98 (MQ=255)

                                                             ccccAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAAt                                      <  2:122867/98‑1 (MQ=255)

                                                                 aaCAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCaa                                  >  2:1434270/1‑98 (MQ=38)

                                                                  aCAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAg                                 <  1:1195559/98‑1 (MQ=39)

                                                                  aCAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAg                                 >  1:1253596/1‑98 (MQ=39)

                                                                  aCAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAg                                 >  1:2061110/1‑98 (MQ=39)

                                                                  aCAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAg                                 <  1:3309399/98‑1 (MQ=39)

                                                                     gagagaAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTg                              >  1:813028/1‑98 (MQ=38)

                                                                     gagagaAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTg                              <  1:3205495/98‑1 (MQ=39)

                                                                      agagaAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTgg                             <  1:2070214/98‑1 (MQ=38)

                                                                       gagaAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACa                                        <  1:1898486/86‑1 (MQ=39)

                                                                       gagaAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTggg                            <  2:529875/98‑1 (MQ=38)

                                                                        agaAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTgggg                           >  2:1426762/1‑98 (MQ=39)

                                                                        agaAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTgggg                           >  2:2670109/1‑98 (MQ=38)

                                                                           aaTAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAg                        >  1:2166951/1‑98 (MQ=255)

                                                                            aTAGGATCACCAATTCCAACATACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAgg                       <  2:1252721/98‑1 (MQ=37)

                                                                              aGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGAc                     <  1:149167/98‑1 (MQ=255)

                                                                              aGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGAc                     <  1:537491/98‑1 (MQ=255)

                                                                               ggATCCCCAATTCCAACAAACGATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCACATTAAACAGACAACGGCAGTTGGGGGAGGAcc                    >  2:1575937/1‑97 (MQ=14)

                                                                                 aTCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATc                  <  2:1763350/98‑1 (MQ=255)

                                                                                  tCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCg                 <  2:1205579/98‑1 (MQ=255)

                                                                                   cACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGa                <  2:3280535/98‑1 (MQ=255)

                                                                                   cACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGa                >  1:2461235/1‑98 (MQ=255)

                                                                                    aCCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGAt               <  1:692622/98‑1 (MQ=255)

                                                                                    aCCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGAt               <  2:3268647/98‑1 (MQ=255)

                                                                                    aCCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGAt               <  1:64347/98‑1 (MQ=255)

                                                                                     ccAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATa              <  2:1378445/98‑1 (MQ=255)

                                                                                      cAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaa             <  2:1568226/98‑1 (MQ=255)

                                                                                      cAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaa             >  1:2674884/1‑98 (MQ=255)

                                                                                       aaTTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaaa            >  2:1843072/1‑98 (MQ=255)

                                                                                       aaTTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaaa            >  1:1653622/1‑98 (MQ=255)

                                                                                       aaTTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaaa            >  1:1427160/1‑98 (MQ=255)

                                                                                        aTTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaaaa           <  1:948740/98‑1 (MQ=255)

                                                                                         ttCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaaaaa          >  2:1801093/1‑98 (MQ=255)

                                                                                          tCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaaaaaa         <  1:1847266/98‑1 (MQ=255)

                                                                                          tCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaaaaaa         <  1:3166061/98‑1 (MQ=255)

                                                                                            caacaaACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAAt       >  2:1723546/1‑98 (MQ=255)

                                                                                             aacaaacaATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATa      <  1:645100/98‑1 (MQ=255)

                                                                                              acaaacaaTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATAc     <  1:2963889/98‑1 (MQ=255)

                                                                                              acaaacaaTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATAc     >  1:3355058/1‑98 (MQ=255)

                                                                                              acaaacaaTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATAc     <  1:1697010/98‑1 (MQ=255)

                                                                                               caaacaaTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATACt    <  1:3235337/98‑1 (MQ=255)

                                                                                                aaacaaTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATACtt   <  2:484099/98‑1 (MQ=255)

                                                                                                aaacaaTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATACtt   <  1:1795658/98‑1 (MQ=255)

                                                                                                aaacaaTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATACtt   <  2:2812320/98‑1 (MQ=255)

                                                                                                 aacaaTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATACttt  <  2:1385734/98‑1 (MQ=255)

GTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATACTTT  >  BX571857/1167200‑1167394

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