BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BX571857	431,505	T→C	10.2%	S77P (TCA→CCA)	SAS0387 →	exotoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BX571857	431,505	0	T	C	10.2%	446.7 / 8.2	137	S77P (TCA→CCA)	SAS0387	exotoxin
Reads supporting (aligned to +/- strand): ref base T (47/70); new base C (8/8); total (56/81)
Fisher's exact test for biased strand distribution p-value = 5.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 139 total aligned reads displayed.

AAATGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGACACAT  >  BX571857/431411‑431601

aaaTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCa                                                                                                <  2:204760/97‑1 (MQ=255)

aaaTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCa                                                                                                <  2:1519182/97‑1 (MQ=255)

aaaTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAACCa                                                                                                >  2:1456484/1‑97 (MQ=255)

aaaTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCAACTAAAGTAGAAGCACCGCAACCa                                                                                                >  2:1456413/1‑97 (MQ=255)

 aaTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCAACTAAAGTAGAAGCACCGCAACCaa                                                                                               >  2:1394239/1‑97 (MQ=255)

 aaTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCAACTAAAGTAGAAACACCGCAACCaa                                                                                               >  2:1116837/1‑97 (MQ=255)

  aTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCaaa                                                                                              <  2:244626/97‑1 (MQ=255)

   tGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCaaaa                                                                                             >  1:12715/1‑97 (MQ=255)

   tGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCaaaa                                                                                             >  1:224624/1‑97 (MQ=255)

    gCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCCAACGCGACAACACCATCTTCAACTAATGTAGAAGCACCGCAATCCAAAc                                                                                            >  1:1834379/1‑97 (MQ=255)

      aacaacGCCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCa                                                                                          <  2:488434/97‑1 (MQ=255)

      aacaacACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAAACAAAACCa                                                                                          <  2:1157043/97‑1 (MQ=255)

        caacaCCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAaccaaaaccaaa                                                                                        >  2:882353/1‑97 (MQ=255)

        caacaCCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCaaa                                                                                        <  2:2135337/97‑1 (MQ=255)

        caacaCCATCTGCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCaaa                                                                                        <  2:638535/97‑1 (MQ=255)

         aacaCCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAAc                                                                                       >  1:1048874/1‑97 (MQ=255)

         aacaCCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAAc                                                                                       <  2:1812566/97‑1 (MQ=255)

           cacCATCTTCAACTAAAGTGGAAGCCCCGCAACAAGCAGCAAACGCGACAACACCATCTTCACCTAAAGTAGAAGCACCGCAATCAAAACCAAAcgc                                                                                     <  2:1535984/97‑1 (MQ=255)

           cacCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAAcgc                                                                                     <  1:1497780/97‑1 (MQ=255)

            acCATCTTCAACTAAAGTGGACGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGACTAAAGTAGACGCACCGCAAGCAAAACCACAcgcc                                                                                    >  2:293700/1‑96 (MQ=14)

            acCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCCATCAAAACCAAAcacc                                                                                    >  2:8168/1‑94 (MQ=255)

              cATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGac                                                                                  <  2:1876140/97‑1 (MQ=255)

                tCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGacaa                                                                                <  1:1579111/97‑1 (MQ=255)

                tCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGacaa                                                                                <  2:2149953/97‑1 (MQ=255)

                 cTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAGCGCGACAACACCAGCTTCAACGAAAGTAGAAGCACCGCAATCAAAACCAAACGCGacaac                                                                               <  1:662190/97‑1 (MQ=255)

                  ttCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGacaaca                                                                              >  1:2151011/1‑97 (MQ=255)

                  ttCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGacaaca                                                                              >  1:90764/1‑97 (MQ=255)

                  tgcAACTAAAGTGGAAGAACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACAGCAACCAAAACCAAACGCGacaaaa                                                                              >  2:2147500/3‑95 (MQ=18)

                   tCAACTAAAGTGAAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAacac                                                                             <  2:571585/97‑1 (MQ=255)

                    cAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACAcc                                                                            <  2:58684/97‑1 (MQ=255)

                    cAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACAcc                                                                            <  1:1973255/97‑1 (MQ=255)

                    cAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACAcc                                                                            >  1:785062/1‑97 (MQ=255)

                      aCTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCAt                                                                          <  1:440606/97‑1 (MQ=255)

                       cTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATc                                                                         >  1:2049977/1‑97 (MQ=255)

                        tAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCt                                                                        >  1:1767909/1‑97 (MQ=255)

                         aaaGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAGGCACCACAATCAAAACCAAAAGCGACAACACAATCtt                                                                       >  2:434828/1‑97 (MQ=21)

                          aaGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTc                                                                      >  1:1995641/1‑97 (MQ=255)

                                aaGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTaa                                                                >  1:1425157/1‑97 (MQ=255)

                                aaGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTaa                                                                <  1:868913/97‑1 (MQ=255)

                                 aGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCGACTAAAGTAGAAGCACCGCCACCACAACCAAACGCGACAACACCACCTGCAGCTaaa                                                               >  2:2180625/1‑97 (MQ=18)

                                 aGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAACCAAAACCAAACGCGACAACACCATCTTCAACTaaa                                                               >  2:1570552/1‑97 (MQ=255)

                                  gCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAg                                                              <  2:1362810/97‑1 (MQ=255)

                                   cACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGt                                                             >  2:1152669/1‑97 (MQ=255)

                                   cACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGt                                                             >  1:611435/1‑97 (MQ=255)

                                    aCCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTg                                                            <  2:1864492/97‑1 (MQ=255)

                                     ccGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTgg                                                           <  1:1365239/97‑1 (MQ=255)

                                       gCAACGCGCAGCAACCGCGACAACACCCTCCTCAACTAAAGTCGACGCACCGCAATCAAAACCCAACGCGACAACACCATCTTCAACTAAAGTGGaa                                                         <  2:623174/97‑1 (MQ=18)

                                       gCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGaa                                                         >  2:1722389/1‑97 (MQ=255)

                                        caacaaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAg                                                        >  1:1122919/1‑97 (MQ=255)

                                        caacaaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAg                                                        <  1:1351074/97‑1 (MQ=255)

                                         aacaaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGc                                                       <  1:686200/97‑1 (MQ=255)

                                         aacaaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGc                                                       <  1:556351/97‑1 (MQ=255)

                                         aacaaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGc                                                       <  1:1748718/97‑1 (MQ=255)

                                          acaaGCAGCAAGCGCGACAACACCAGCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCa                                                      <  2:402707/97‑1 (MQ=255)

                                          acaaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCa                                                      >  2:1302366/1‑97 (MQ=255)

                                           caaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCAc                                                     >  2:32840/1‑97 (MQ=255)

                                              gcagcaAACGCGACAACACCCTCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGc                                                  <  2:125538/97‑1 (MQ=255)

                                                agcaAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGcaa                                                >  1:1203530/1‑97 (MQ=255)

                                                agcaAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGcaa                                                >  1:1244096/1‑97 (MQ=255)

                                                         gACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAgcagca                                       <  1:764016/97‑1 (MQ=255)

                                                            aacaCCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAAc                                    >  1:1384133/1‑97 (MQ=255)

                                                            aacaCCACCTTCACCTAAGGTGGAAGCACCGCACCCAGAGCCAAACGCGCCAACACCAGCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAAc                                    <  2:1553342/97‑1 (MQ=14)

                                                               acCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAAcgcg                                 >  1:1102147/1‑97 (MQ=255)

                                                               acCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAAcgcg                                 >  1:203779/1‑97 (MQ=255)

                                                                ccATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGa                                >  1:640009/1‑97 (MQ=255)

                                                                ccATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGa                                >  1:640010/1‑97 (MQ=255)

                                                                 cATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGac                               <  1:333064/97‑1 (MQ=255)

                                                                 cATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGac                               <  1:381498/97‑1 (MQ=255)

                                                                 cATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGac                               <  1:136488/97‑1 (MQ=255)

                                                                   tCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGacaa                             <  1:829125/97‑1 (MQ=255)

                                                                    cTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGacaac                            >  1:312861/1‑97 (MQ=255)

                                                                    cTTCAACTAAAGTAGAAGCACCGCAACCAAAAGCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGacaac                            <  1:1357326/97‑1 (MQ=255)

                                                                      tCAACTAAAGTAGAAGCACCGCAACCAGAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAacac                          <  2:811861/97‑1 (MQ=255)

                                                                        aaCTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAAcacca                        <  2:536871/97‑1 (MQ=255)

                                                                          cTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACaccacc                      <  1:1992214/97‑1 (MQ=255)

                                                                          cTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACaccacc                      <  1:1992112/97‑1 (MQ=255)

                                                                          cTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACaccacc                      >  2:2165407/1‑97 (MQ=255)

                                                                           tAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCt                     >  1:967331/1‑97 (MQ=255)

                                                                            aaaGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCtt                    <  2:1717734/97‑1 (MQ=255)

                                                                             aaGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTc                   >  2:890302/1‑97 (MQ=255)

                                                                             aaGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTc                   >  1:1838716/1‑97 (MQ=255)

                                                                              aGTGGAAGCACCGCAACCAACACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCg                  <  2:1101699/97‑1 (MQ=255)

                                                                              aGTCGAAGCGCCTCAACCAAAACCAAACGCCACAACACCATCTTCAACTAAAGTGGAAGCACCGCATCAAGCAGCAAACGCGACATCACCACCTGCg                  <  2:576489/97‑1 (MQ=17)

                                                                              aGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCg                  <  1:1630667/97‑1 (MQ=255)

                                                                               gTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGt                 >  1:90581/1‑97 (MQ=255)

                                                                               gTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGt                 >  1:171086/1‑97 (MQ=255)

                                                                                 ggAAGCACCGCAATCAACACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCt               <  2:2040101/96‑1 (MQ=255)

                                                                                 aGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCt               >  1:794538/1‑97 (MQ=255)

                                                                                  gAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTa              <  2:1794175/97‑1 (MQ=255)

                                                                                  gAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTa              >  2:1695671/1‑97 (MQ=255)

                                                                                    aGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAAt            <  1:968470/97‑1 (MQ=255)

                                                                                    aGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAAt            >  1:232283/1‑97 (MQ=255)

                                                                                      cACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGt          >  2:1041752/1‑97 (MQ=255)

                                                                                      cACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGt          >  1:700284/1‑97 (MQ=255)

                                                                                       cccGCAATCAACACCAGACGCGACAACGCCATCTTCAACTCAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTa         <  2:1667362/96‑1 (MQ=18)

                                                                                       cccGCAACCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTa         <  2:1881975/96‑1 (MQ=255)

                                                                                       cccGCAACAAAGACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTa         <  2:796374/96‑1 (MQ=17)

                                                                                            aaTCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGacac    <  1:863995/97‑1 (MQ=255)

                                                                                             aTCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGAcaca   <  1:427344/97‑1 (MQ=255)

                                                                                              tCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGACACAt  >  2:45467/1‑97 (MQ=255)

AAATGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGACACAT  >  BX571857/431411‑431601

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