BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,596,823	T→C	L21L (TTA→TTG)	nlpB ←	lipoprotein‑34

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,596,823	0	T	C	99.4%	497.5 / ‑4.8	157	L21L (TTA→TTG)	nlpB	lipoprotein‑34
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (77/79); total (77/80)
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.31e-01

Only 100 of 168 total aligned reads displayed.

ATGAAGCTCCGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGTAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAACTTTACAGCAAAC  >  NC_000913/2596733‑2596912

aTGAAGCTCCGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGcaa                                                                                         >  1:3155888/1‑93 (MQ=255)

aTGAAGCTCCGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGcaa                                                                                         >  1:2195920/1‑93 (MQ=255)

 tGAAGCTCCGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGcaat                                                                                        >  2:1295926/1‑93 (MQ=255)

 tGAAGCTCCGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGcaat                                                                                        >  2:2802244/1‑93 (MQ=255)

 tGAAGCTCCGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGcaat                                                                                        <  1:3647459/93‑1 (MQ=255)

  gAAGCTCCGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGcaata                                                                                       >  2:924635/1‑93 (MQ=255)

  gAAGCTCCGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGcaata                                                                                       <  2:147868/93‑1 (MQ=255)

  gAAGCTCCGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGcaata                                                                                       >  1:1292639/1‑93 (MQ=255)

  gAAGCGCCGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGcaata                                                                                       <  1:278250/93‑1 (MQ=255)

    aGCTCCGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATaaa                                                                                     <  2:4409503/93‑1 (MQ=255)

     gCTCCGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATaaaa                                                                                    <  1:69567/93‑1 (MQ=255)

     gCTCCGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATaaaa                                                                                    >  1:4083795/1‑93 (MQ=255)

       tCCGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACa                                                                                  >  2:3507161/1‑93 (MQ=255)

       tCCGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACa                                                                                  <  2:759195/93‑1 (MQ=255)

        ccGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACaa                                                                                 >  1:3964035/1‑93 (MQ=255)

         cGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAg                                                                                <  1:3265212/93‑1 (MQ=255)

          gCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGc                                                                               >  2:4207000/1‑93 (MQ=255)

             aGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGaa                                                                            >  2:3700465/1‑93 (MQ=255)

              gCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGaaa                                                                           >  2:631712/1‑93 (MQ=255)

              gCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGaaa                                                                           <  1:567363/93‑1 (MQ=255)

               cGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGCCCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAac                                                                          <  1:1247770/93‑1 (MQ=255)

                ggTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAaca                                                                         <  2:4128677/93‑1 (MQ=255)

                ggTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAaca                                                                         <  2:2198468/93‑1 (MQ=255)

                   gCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACAccc                                                                      <  2:4531130/93‑1 (MQ=255)

                   gCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACAccc                                                                      <  1:1757511/93‑1 (MQ=255)

                       cTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCaa                                                                  <  1:4440837/93‑1 (MQ=255)

                       cTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCaa                                                                  >  1:32507/1‑93 (MQ=255)

                        ttCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAAc                                                                 >  2:1882932/1‑93 (MQ=255)

                        ttCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAAc                                                                 >  2:4419481/1‑93 (MQ=255)

                          ccAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCt                                                               <  2:3409424/93‑1 (MQ=255)

                           cAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCtt                                                              >  1:1405384/1‑93 (MQ=255)

                              gTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGc                                                           >  2:723192/1‑93 (MQ=255)

                               tAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGcc                                                          <  1:2890583/93‑1 (MQ=255)

                                aGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCa                                                         <  1:3855434/93‑1 (MQ=255)

                                aGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCa                                                         <  1:757148/93‑1 (MQ=255)

                                 ggCTTCATCACCACTGACCTGACGCTTATAGCGTGCGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAg                                                        >  2:799403/1‑93 (MQ=255)

                                 ggCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAg                                                        <  1:3724577/93‑1 (MQ=255)

                                  gCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAgg                                                       >  1:487273/1‑93 (MQ=255)

                                   cTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGgc                                                      >  2:2233063/1‑93 (MQ=255)

                                   cTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGgc                                                      <  2:3561721/93‑1 (MQ=255)

                                    ttCATCACCACTGACCTGGCGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGgcc                                                     >  2:4122064/1‑92 (MQ=255)

                                    ttCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGgcg                                                     <  2:2903688/93‑1 (MQ=255)

                                       atcaCCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGa                                                  >  1:3655009/1‑93 (MQ=255)

                                          accacTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACtt                                               >  2:2232432/1‑93 (MQ=255)

                                          accacTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACtt                                               <  1:687081/93‑1 (MQ=255)

                                             acTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTg                                            >  2:1492784/1‑93 (MQ=255)

                                              ctgacctgacGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACCTTTGa                                           >  2:966006/1‑93 (MQ=255)

                                               tgacctgacGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGaa                                          <  1:2074304/93‑1 (MQ=255)

                                               tgacctgacGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAACGAAACACCCGCAACCTTTGCCAGGCGCGACTTGTGaa                                          >  2:1264989/1‑93 (MQ=255)

                                                gacctgacGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAAc                                         <  2:49065/93‑1 (MQ=255)

                                                 acctgacGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACa                                        <  1:1097070/93‑1 (MQ=255)

                                                 acctgacGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACa                                        >  2:1557884/1‑93 (MQ=255)

                                                 acctgacGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACa                                        >  1:2010602/1‑93 (MQ=255)

                                                 acctgacGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCAAGGCGCGACTTTTGAACa                                        <  2:1359919/93‑1 (MQ=255)

                                                  cctgacGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACag                                       <  2:2596191/93‑1 (MQ=255)

                                                    tgacGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACagag                                     >  1:3024531/1‑93 (MQ=255)

                                                    tgacGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACagag                                     >  2:158328/1‑93 (MQ=255)

                                                      acGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTa                                   <  1:1232573/93‑1 (MQ=255)

                                                       cGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTaa                                  >  1:2319435/1‑93 (MQ=255)

                                                       cGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTaa                                  <  2:3698971/93‑1 (MQ=255)

                                                       cGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTaa                                  >  1:1468501/1‑93 (MQ=255)

                                                       cGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTaa                                  <  2:889142/93‑1 (MQ=255)

                                                       cGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTaa                                  <  1:668466/93‑1 (MQ=255)

                                                        gCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAg                                 <  2:1878084/93‑1 (MQ=255)

                                                        gCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAg                                 <  1:264512/93‑1 (MQ=255)

                                                         cTTATAGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGc                                <  2:1602081/93‑1 (MQ=255)

                                                          ttATAGCGTGAGTCAGAACCACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGcc                               >  2:4580614/1‑93 (MQ=255)

                                                           tataGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCa                              >  2:210364/1‑93 (MQ=255)

                                                             taGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATc                            <  1:4067987/93‑1 (MQ=255)

                                                             taGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATc                            <  1:4207000/93‑1 (MQ=255)

                                                              aGCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCa                           <  1:4317411/93‑1 (MQ=255)

                                                               gCGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCaa                          >  1:987198/1‑93 (MQ=255)

                                                                cGTGAGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCaaa                         <  2:2131467/93‑1 (MQ=255)

                                                                    aGTCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAAtctc                     <  2:1921837/93‑1 (MQ=255)

                                                                      tCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTccc                   <  1:4114292/93‑1 (MQ=255)

                                                                      tCAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTccc                   <  2:1693971/93‑1 (MQ=255)

                                                                       cAGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCt                  <  1:3166978/93‑1 (MQ=255)

                                                                        aGAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTa                 <  1:3031723/93‑1 (MQ=255)

                                                                         gAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTaa                >  1:1562106/1‑93 (MQ=255)

                                                                         gAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTaa                >  1:908163/1‑93 (MQ=255)

                                                                         gAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTaa                <  1:3880059/93‑1 (MQ=255)

                                                                         gAACTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTaa                >  2:2356743/1‑93 (MQ=255)

                                                                          aaCTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTaaa               <  2:2630778/93‑1 (MQ=255)

                                                                           aCTACAGGCAGCGAGTAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAAc              <  1:2841593/93‑1 (MQ=255)

                                                                           aCTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAAc              >  2:4490967/1‑93 (MQ=255)

                                                                           aCTACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAAc              <  1:1130875/93‑1 (MQ=255)

                                                                             tACAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCTACTTTTGAACAGAGTAAGCCCTCAAATCTCCCTCCACtt            >  1:2120856/1‑93 (MQ=255)

                                                                               cAGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAACTTTa          >  1:2530496/1‑93 (MQ=255)

                                                                                aGGCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAACTTTAc         >  1:3399880/1‑93 (MQ=255)

                                                                                 ggCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAACTTTACa        >  1:2955181/1‑93 (MQ=255)

                                                                                 ggCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAACTTTACa        >  2:3099792/1‑93 (MQ=255)

                                                                                  gCAGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAACTTTACAg       >  2:3017983/1‑93 (MQ=255)

                                                                                    aGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAACTTTACAGCa     >  2:881400/1‑93 (MQ=255)

                                                                                    aGCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAACTTTACAGCa     <  1:952640/93‑1 (MQ=255)

                                                                                     gCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAACTTTACAGCaa    >  2:536967/1‑93 (MQ=255)

                                                                                     gCGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAACTTTACAGCaa    >  2:1191937/1‑93 (MQ=255)

                                                                                      cGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAACTTTACAGCaaa   >  2:3791495/1‑93 (MQ=255)

                                                                                      cGAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAACTTTACAGCaaa   >  2:556981/1‑93 (MQ=255)

                                                                                       gAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAACTTTACAGCAAAc  >  1:2272997/1‑93 (MQ=255)

                                                                                       gAGCAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAACTTTACAGCAAAc  <  2:3161126/93‑1 (MQ=255)

ATGAAGCTCCGCAAGCGGTGCCGCTTCCAGGTAGGCTTCATCACCACTGACCTGACGCTTATAGCGTGAGTCAGAACTACAGGCAGCGAGTAATAAAACAAGCGAAACACCCGCAACCTTTGCCAGGCGCGACTTTTGAACAGAGTAAGCCATCAAATCTCCCTAAACTTTACAGCAAAC  >  NC_000913/2596733‑2596912

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