BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BX571857	1,167,297	A→G	6.1%	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	6.1%	438.4 / 8.3	120	K71R (AAA→AGA)	SAS1093	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (59/52); new base G (3/6); total (62/58)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

GGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATACTT  >  BX571857/1167202‑1167393

ggggCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCaacaaa                                                                                                 >  1:248640/1‑97 (MQ=255)

ggggCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCaacaaa                                                                                                 >  2:1799703/1‑97 (MQ=255)

  ggCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCaacaaaca                                                                                               <  1:726142/97‑1 (MQ=255)

   gCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAacaaacaa                                                                                              <  1:2292039/97‑1 (MQ=255)

    cccAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAAt                                                                                             >  2:813282/1‑97 (MQ=255)

     ccAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATg                                                                                            >  2:82558/1‑97 (MQ=255)

      cAACACAGAAGCTGGCGAAAAGTCAGCATACAACAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGc                                                                                           <  2:1105723/97‑1 (MQ=255)

         cacaGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAg                                                                                        >  1:1419658/1‑97 (MQ=255)

          acaGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGt                                                                                       >  2:1323726/1‑97 (MQ=255)

           caGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGtt                                                                                      >  1:510760/1‑97 (MQ=255)

             gAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTgg                                                                                    >  1:1448560/1‑97 (MQ=255)

             gAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTgg                                                                                    >  1:1448576/1‑97 (MQ=255)

             gAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTgg                                                                                    >  1:2534972/1‑97 (MQ=255)

              aaGCTGGAGAAAAGCCAGCATACAAGAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGc                                                                                   <  2:1296410/97‑1 (MQ=255)

                 cTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCggg                                                                                >  2:1951502/1‑97 (MQ=255)

                 cTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCggg                                                                                >  2:2178413/1‑97 (MQ=255)

                 cTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCggg                                                                                >  1:98223/1‑97 (MQ=255)

                      gAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAANTAGNATCACCAATTCCAACNAACAATGCAAGTTGGNGGGGcccc                                                                           >  2:2541708/1‑97 (MQ=255)

                      gAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGcccc                                                                           >  1:21189/1‑97 (MQ=255)

                       ataaGTCTACATACAATTATTTGCAAGCTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCa                                                                          <  2:1641616/95‑1 (MQ=17)

                        aaaGTCAACAAACAGTACTGTCCATGTTCTCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCaa                                                                         <  2:1949307/97‑1 (MQ=14)

                          agccAGCATACAGAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACa                                                                       <  2:247180/94‑1 (MQ=255)

                           gTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACAt                                                                      >  1:2737276/1‑97 (MQ=255)

                             cAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATag                                                                    >  1:1127325/1‑97 (MQ=255)

                              aGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATaga                                                                   >  1:1510223/1‑97 (MQ=255)

                                 aTACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGaa                                                                >  1:2863854/1‑97 (MQ=255)

                                 aTACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGaa                                                                <  1:1507036/97‑1 (MQ=255)

                                  tACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGaaa                                                               >  1:1905722/1‑97 (MQ=255)

                                   aCAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAAt                                                              >  1:1902822/1‑97 (MQ=255)

                                   aCAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAAt                                                              >  2:518206/1‑97 (MQ=255)

                                   aCAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAAt                                                              >  1:676067/1‑97 (MQ=255)

                                    cAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAAtt                                                             >  2:2062720/1‑97 (MQ=255)

                                    cAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAAtt                                                             >  2:2495856/1‑97 (MQ=255)

                                    cAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAAtt                                                             >  2:407480/1‑97 (MQ=255)

                                    aagaaaTGTGCAAACTCGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAAtt                                                             <  2:2534160/94‑1 (MQ=14)

                                        aaTGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGt                                                         >  2:2551546/1‑97 (MQ=255)

                                         aTGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTc                                                        >  1:2865459/1‑97 (MQ=255)

                                          tgtgCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCa                                                       >  2:2764617/1‑97 (MQ=255)

                                                   tGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCa                                              >  2:242716/1‑97 (MQ=255)

                                                    ggCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCaa                                             >  1:2010519/1‑97 (MQ=255)

                                                    ggCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCaa                                             <  1:2679620/97‑1 (MQ=255)

                                                    ggCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCaa                                             >  2:2099443/1‑97 (MQ=255)

                                                    ggCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCaa                                             <  1:1269373/97‑1 (MQ=255)

                                                         ggCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGAc                                        <  2:1635136/97‑1 (MQ=255)

                                                         ggCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGAc                                        <  1:1907994/97‑1 (MQ=255)

                                                          gCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACa                                       >  2:550205/1‑97 (MQ=255)

                                                          gCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGACCCAACATAGAGAAATTGGGGCACCAATTTAAACAGaaa                                       >  2:144714/1‑95 (MQ=14)

                                                              cAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGc                                   <  1:886445/97‑1 (MQ=38)

                                                               aaCAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCa                                  <  1:1590620/97‑1 (MQ=39)

                                                                  agagagAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGt                               <  1:2850987/97‑1 (MQ=38)

                                                                  agagagAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGt                               <  1:2684835/97‑1 (MQ=38)

                                                                  agagagAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGt                               >  1:1833851/1‑97 (MQ=39)

                                                                   gcgagcAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGtt                              <  1:2117681/95‑1 (MQ=255)

                                                                      agaAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTggg                           <  1:884627/97‑1 (MQ=39)

                                                                      agaAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTggg                           <  2:1191253/97‑1 (MQ=38)

                                                                        aaaTAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGt                         <  2:2059164/97‑1 (MQ=37)

                                                                        aaaTAGGATCACCAATTCCAACAAACAAGGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGt                         <  2:740053/97‑1 (MQ=37)

                                                                         aaTAGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTa                        >  2:52006/1‑97 (MQ=1)

                                                                         aaTAGGATCACCAATTCCAACAAACAATGCAAGTTGGAGGGGCCACAACATAGAGAAATTGGGGCAACAACTTCAACAGACCATGCAAGTGGGGGTa                        >  2:279931/1‑97 (MQ=14)

                                                                            aGGATCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGa                     <  2:2386870/97‑1 (MQ=25)

                                                                            aGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGa                     <  1:1044990/97‑1 (MQ=255)

                                                                            aGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGa                     <  1:386102/97‑1 (MQ=255)

                                                                            aGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGa                     >  1:794262/1‑97 (MQ=255)

                                                                            aGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGa                     <  1:1044981/97‑1 (MQ=255)

                                                                              gATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACa                   <  1:1529432/97‑1 (MQ=255)

                                                                              gATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACa                   <  1:276820/97‑1 (MQ=255)

                                                                              gATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACa                   <  2:1296654/97‑1 (MQ=255)

                                                                               aTCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACAt                  >  2:315730/1‑97 (MQ=255)

                                                                               aTCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACAt                  >  2:1090412/1‑97 (MQ=255)

                                                                                tCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATc                 <  1:108758/97‑1 (MQ=255)

                                                                                tCACCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATc                 <  1:1330965/97‑1 (MQ=255)

                                                                                tCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATc                 <  1:1558775/97‑1 (MQ=255)

                                                                                tCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATc                 <  2:1723560/97‑1 (MQ=255)

                                                                                  aCCCATGCCAACAGACCATGCCAGTTGGCGCGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAGGTTGGGGTAGGACATCGa               <  2:2790318/97‑1 (MQ=14)

                                                                                  aCCAATTCCAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGa               <  2:2977968/97‑1 (MQ=255)

                                                                                  aCCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGa               >  2:1605904/1‑97 (MQ=255)

                                                                                   ccAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGAt              >  1:112918/1‑97 (MQ=255)

                                                                                   ccAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGAt              <  1:1216283/97‑1 (MQ=255)

                                                                                     aaTTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGAc                                        <  1:2459846/69‑1 (MQ=40)

                                                                                     aaTTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaa            <  1:1336710/97‑1 (MQ=255)

                                                                                     aaTTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaa            >  1:474618/1‑97 (MQ=255)

                                                                                     aaTTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaa            >  1:414485/1‑97 (MQ=255)

                                                                                      aTTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaaa           <  2:1715976/97‑1 (MQ=255)

                                                                                       ttCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaaaa          <  1:815968/97‑1 (MQ=255)

                                                                                        tCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaaaaa         >  2:1471003/1‑97 (MQ=255)

                                                                                        tCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaaaaa         >  1:1166339/1‑97 (MQ=255)

                                                                                        cccAACAGACAAAGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaaaaa         <  2:1863145/96‑1 (MQ=255)

                                                                                         ccAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaaaaaa        >  2:668428/1‑97 (MQ=255)

                                                                                         ccAACAGACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACACCGCTaaaaaa        >  2:1802362/1‑97 (MQ=11)

                                                                                         ccAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaaaaaa        <  2:595494/97‑1 (MQ=255)

                                                                                          caacaaACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaaaaaaa       <  2:1151698/97‑1 (MQ=255)

                                                                                          caacaaACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATaaaaaaa       <  2:1151699/97‑1 (MQ=255)

                                                                                           aacagacaATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAAt      <  1:1113489/97‑1 (MQ=255)

                                                                                           aacaaacaATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAAt      >  2:586984/1‑97 (MQ=255)

                                                                                           aacaaacaATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAAt      <  1:2769441/97‑1 (MQ=255)

                                                                                            acaaacaaTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATa     >  1:2335309/1‑97 (MQ=255)

                                                                                            acaaacaaTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATa     <  2:2131771/97‑1 (MQ=255)

                                                                                             cgaacaaTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACACTGCAAGTTGGGGTAGGACATCGATAAAAAAATAc    <  1:2156659/95‑1 (MQ=255)

                                                                                             caaacaaTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATAc    >  2:1660482/1‑97 (MQ=255)

                                                                                             caaacaaTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATAc    <  2:2313523/97‑1 (MQ=255)

GGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGGTCACCAATTTCAACAGACAATGCAAGTTGGGGTAGGACATCGATAAAAAAATACTT  >  BX571857/1167202‑1167393

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