BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	569,160	0	G	T	22.0%	446.1 / 92.8	191	D87Y (GAC→TAC)	ybcK	DLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand): ref base G (71/78); new base T (20/22); total (91/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 198 total aligned reads displayed.

AGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTC  >  NC_000913/569063‑569258

aGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTATATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGa                                                                                                   <  2:2725504/99‑1 (MQ=255)

aGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGa                                                                                                   <  1:2341222/99‑1 (MQ=255)

  gAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACa                                                                                                 >  1:534662/1‑99 (MQ=255)

  gAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATCGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACa                                                                                                 <  1:638122/99‑1 (MQ=255)

        aTGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACttt                                                                                           <  2:390381/99‑1 (MQ=255)

        aTGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACttt                                                                                           >  1:310876/1‑99 (MQ=255)

        aTGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACttt                                                                                           <  1:1628558/99‑1 (MQ=255)

         tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTCTAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTc                                                                                          <  2:3116040/99‑1 (MQ=255)

         tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTc                                                                                          <  2:76866/99‑1 (MQ=255)

         tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTc                                                                                          <  2:63139/99‑1 (MQ=255)

         tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTc                                                                                          <  1:425394/99‑1 (MQ=255)

         agCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTc                                                                                          <  2:3270958/98‑1 (MQ=255)

          gCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                         >  2:2473853/1‑99 (MQ=255)

           cacaATTAGGAGCATTTTCGGAATTTTTAGATGCTATAGAGCATGGTTAGATATTGCCAGGAACTACATTGATAGTTGAAAGTCAGTCCAGACTTTCaa                                                                                        <  2:1629201/99‑1 (MQ=255)

           cacaATCAGTAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTTAAAGTCTGTCCAGACTTTCaa                                                                                        <  1:2613762/99‑1 (MQ=255)

           cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaa                                                                                        <  2:2626336/99‑1 (MQ=255)

           cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaa                                                                                        <  2:134840/99‑1 (MQ=255)

           cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaa                                                                                        <  1:1935700/99‑1 (MQ=255)

           cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaa                                                                                        <  1:533773/99‑1 (MQ=255)

           cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaa                                                                                        <  2:702348/99‑1 (MQ=255)

           cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaa                                                                                        <  1:1211303/99‑1 (MQ=255)

           cacaAACAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATCTATTGCCCGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaa                                                                                        <  1:3299688/99‑1 (MQ=255)

            acaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                       <  1:293271/99‑2 (MQ=255)

              aaTCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAagag                                                                                     <  1:1799666/99‑1 (MQ=255)

                tCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGaa                                                                                   >  1:1603340/1‑99 (MQ=255)

                    gAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAg                                                                               >  1:1567563/1‑99 (MQ=255)

                       cTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCg                                                                            <  1:1244695/99‑1 (MQ=255)

                        tttttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCgg                                                                           <  1:23142/99‑1 (MQ=255)

                        tttttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCgg                                                                           <  1:597572/99‑1 (MQ=255)

                         ttttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGt                                                                          <  2:1074472/99‑1 (MQ=255)

                           ttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGa                                                                        <  1:1569907/99‑1 (MQ=255)

                             cGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAg                                                                      <  1:1288705/99‑1 (MQ=255)

                              ggAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGc                                                                     >  1:2968259/1‑99 (MQ=255)

                              ggAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGc                                                                     >  1:1507005/1‑99 (MQ=255)

                              ggAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGc                                                                     >  1:50380/1‑99 (MQ=255)

                                aaTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGa                                                                   >  1:441823/1‑99 (MQ=255)

                                 aTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAt                                                                  <  2:1281486/99‑1 (MQ=255)

                                 aTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAt                                                                  <  2:15270/99‑1 (MQ=255)

                                  tttttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAtt                                                                 <  2:1122877/99‑1 (MQ=255)

                                  tttttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAtt                                                                 <  2:3057173/99‑1 (MQ=255)

                                  gttttAGATGTTATAGAGCATGGTTATATATTGCCTGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAtt                                                                 <  2:2800555/98‑1 (MQ=255)

                                     ttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGACAAAGTCGGTGAAGCGATTGag                                                              >  2:2584191/1‑98 (MQ=255)

                                     ttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGaa                                                              >  2:1564831/1‑99 (MQ=255)

                                          tGCTATAGAGCATGGTTATATATTGCCATGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCt                                                         <  1:1120822/99‑1 (MQ=255)

                                             tataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGaa                                                      >  2:1357153/1‑99 (MQ=255)

                                              ataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGaaa                                                     >  1:2574844/1‑99 (MQ=255)

                                                  agCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGa                                                 <  2:1577986/99‑1 (MQ=255)

                                                   gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                <  2:889260/99‑1 (MQ=255)

                                                   gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                <  1:1012726/99‑1 (MQ=255)

                                                    cATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAtt                                               >  2:2132879/1‑99 (MQ=255)

                                                     aTGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAttt                                              <  2:237183/99‑1 (MQ=255)

                                                     aTGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAttt                                              <  1:2520970/99‑1 (MQ=255)

                                                     aTGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAttt                                              <  1:1171087/99‑1 (MQ=255)

                                                        gTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGa                                           >  1:2772195/1‑99 (MQ=255)

                                                         ttATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGaa                                          >  1:2736726/1‑99 (MQ=255)

                                                         atatataTTGCAAGGAACTACAATGTTAGTAGAAAATCTAGAAAAAAATTCACGAGCAAAAGAAGGTGACGAGAATGAACGAATGAAATTGATTTTGaa                                          >  2:1184423/2‑99 (MQ=255)

                                                          tatattTTGTCATGAACTACATGGTTAGTTGAAAGTCGGTACAGACTTTCAAGAGAAAAAGGTGGGTAAGCTATTGAACGTCTGAAATTGATTTTGAAt                                         <  2:2902695/99‑1 (MQ=255)

                                                          tatataTTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTTAAt                                         >  2:826668/1‑99 (MQ=255)

                                                            tatatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCa                                       <  1:3244406/99‑1 (MQ=255)

                                                             atatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCAc                                      >  1:2628840/1‑99 (MQ=255)

                                                               atTGCCAGGAACTACATTGTTAGGTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACgg                                    <  2:3180343/99‑1 (MQ=255)

                                                               atTGCCAGGAACTACATTCTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACgg                                    >  1:2460202/1‑99 (MQ=255)

                                                                  gCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTAt                                 >  1:1778044/1‑99 (MQ=255)

                                                                   ccAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTAtt                                <  1:2118648/99‑1 (MQ=255)

                                                                         aCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGtt                          <  1:1002055/99‑1 (MQ=255)

                                                                           tACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtat                        >  1:1069239/1‑99 (MQ=255)

                                                                           tACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtat                        <  2:1549850/99‑1 (MQ=255)

                                                                           tACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtat                        >  2:1466845/1‑99 (MQ=255)

                                                                           tACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtat                        <  2:342539/99‑1 (MQ=255)

                                                                             cATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATaa                      <  1:3214551/99‑1 (MQ=255)

                                                                                 gTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAActct                  <  2:1323392/99‑1 (MQ=255)

                                                                                 gTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAActct                  >  1:3411995/1‑99 (MQ=255)

                                                                                  ttAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                 >  1:2496540/1‑99 (MQ=255)

                                                                                   ctcttGAAAGTCTGTCCAGACTTTCCAGAGAAAAAGTAGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCttt                >  1:269766/4‑99 (MQ=38)

                                                                                   ctcttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCttt                >  1:983585/4‑99 (MQ=39)

                                                                                    tcttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg               >  2:3153295/3‑99 (MQ=39)

                                                                                    aGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg               <  2:1981632/99‑1 (MQ=255)

                                                                                     gTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGc              >  2:962953/1‑99 (MQ=255)

                                                                                     gTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAAATCTTTGc              >  2:1869294/1‑99 (MQ=255)

                                                                                      ttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg             >  1:2691929/1‑99 (MQ=255)

                                                                                      ttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg             >  1:3061865/1‑99 (MQ=255)

                                                                                      ttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg             >  1:1685665/1‑99 (MQ=255)

                                                                                      ttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg             >  2:505499/1‑99 (MQ=255)

                                                                                      ttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg             >  2:2820166/1‑99 (MQ=255)

                                                                                      ttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCg             >  2:179027/1‑99 (MQ=255)

                                                                                      ttGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTgtg             >  2:2874458/1‑97 (MQ=255)

                                                                                       tGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGa            >  1:2654360/1‑99 (MQ=255)

                                                                                        gAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc           >  1:2195203/1‑99 (MQ=255)

                                                                                        gAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc           >  2:1690795/1‑99 (MQ=255)

                                                                                        gAAAGTCTGGACAGACNTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc           >  1:289411/1‑99 (MQ=255)

                                                                                         aaaGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACa          >  1:2798764/1‑99 (MQ=255)

                                                                                            gTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATa       <  1:2534263/99‑1 (MQ=255)

                                                                                            gTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATa       <  1:975833/99‑1 (MQ=255)

                                                                                             tCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATAc      >  1:2804299/1‑99 (MQ=255)

                                                                                             tCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATAc      >  1:1713781/1‑99 (MQ=255)

                                                                                                ggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATACCTCTTTGCGACAATACAGt   <  2:1014628/99‑1 (MQ=255)

                                                                                                ggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGt   >  1:2058448/1‑99 (MQ=255)

                                                                                                 gACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc  >  1:1901268/1‑99 (MQ=255)

                                                                                                 gACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc  >  2:2763335/1‑99 (MQ=255)

                                                                                                 gACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc  >  1:706873/1‑99 (MQ=255)

AGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTC  >  NC_000913/569063‑569258

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