BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,318,996	G→A	D312N (GAC→AAC)	atoS →	AtoS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,318,996	0	G	A	99.6%	844.5 / ‑5.7	233	D312N (GAC→AAC)	atoS	AtoS
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (112/120); total (112/121)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.94e-01

Only 100 of 239 total aligned reads displayed.

GGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCGTGAC  >  NC_000913/2318906‑2319085

 gTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACa                                                                                        <  1:5164650/93‑1 (MQ=255)

  tGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTcaacaa                                                                                       <  1:6489453/93‑1 (MQ=255)

  tGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTcaacaa                                                                                       <  1:6348390/93‑1 (MQ=255)

    aTGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATa                                                                                     >  2:6257310/1‑93 (MQ=255)

    aTGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATa                                                                                     >  2:3976344/1‑93 (MQ=255)

    aTGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATa                                                                                     <  1:1629381/93‑1 (MQ=255)

     tGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATAc                                                                                    >  1:7393824/1‑93 (MQ=255)

     tGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATAc                                                                                    >  1:3378936/1‑93 (MQ=255)

       tAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTc                                                                                  <  2:4713590/93‑1 (MQ=255)

       tAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTc                                                                                  >  1:2896473/1‑93 (MQ=255)

        aaCCACAATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCa                                                                                 <  1:2917124/93‑1 (MQ=255)

         accaccATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAg                                                                                <  2:2541343/93‑1 (MQ=255)

         accaccATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAg                                                                                >  2:1250915/1‑93 (MQ=255)

          ccaccaTGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGt                                                                               >  1:4101446/1‑93 (MQ=255)

           caccaTGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGtt                                                                              <  1:1940637/93‑1 (MQ=255)

           caccaTGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGtt                                                                              >  1:7917765/1‑93 (MQ=255)

             ccaTGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCt                                                                            <  2:269273/93‑1 (MQ=255)

              caTGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTa                                                                           >  2:3632890/1‑93 (MQ=255)

                tGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACa                                                                         <  1:6772630/93‑1 (MQ=255)

                tGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACa                                                                         <  2:6787510/93‑1 (MQ=255)

                tGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACa                                                                         <  2:4342372/93‑1 (MQ=255)

                   aCCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTc                                                                      >  2:1999634/1‑93 (MQ=255)

                    cccAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTcc                                                                     >  2:6203423/1‑93 (MQ=255)

                       agcagcagAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGt                                                                  <  1:6630212/93‑1 (MQ=255)

                        gcagcagAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTa                                                                 >  2:3323868/1‑93 (MQ=255)

                         cagcagAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTAc                                                                <  2:7159167/93‑1 (MQ=255)

                          agcagAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACt                                                               >  2:6073328/1‑93 (MQ=255)

                          agcagAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACt                                                               >  2:5319690/1‑93 (MQ=255)

                           gcagAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTg                                                              <  1:5987104/93‑1 (MQ=255)

                           gcagAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTg                                                              <  2:810200/93‑1 (MQ=255)

                             agAAGTTATCACTGGCTATCAACGCCATGAGCTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGa                                                            <  2:4217772/93‑1 (MQ=255)

                             agAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGa                                                            <  2:7480863/93‑1 (MQ=255)

                             agAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGa                                                            >  2:6338117/1‑93 (MQ=255)

                             agAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGa                                                            <  1:5004490/93‑1 (MQ=255)

                              gAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGAt                                                           <  2:2995320/93‑1 (MQ=255)

                              gAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGAt                                                           >  1:3563379/1‑93 (MQ=255)

                               aaGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATa                                                          <  2:4612504/93‑1 (MQ=255)

                                 gTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACg                                                        >  1:2682428/1‑93 (MQ=255)

                                  ttATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTAGTGGCTaacc                                                       >  1:4286223/1‑90 (MQ=255)

                                  ttATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGc                                                       >  1:3459487/1‑93 (MQ=255)

                                  ttATCACTGGCTATCAACGCCATGAACTGGAAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGc                                                       >  2:1703887/1‑93 (MQ=255)

                                   tATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCt                                                      >  2:1369716/1‑93 (MQ=255)

                                   tATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCt                                                      >  2:2754619/1‑93 (MQ=255)

                                    aTCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTg                                                     >  1:6737514/1‑93 (MQ=255)

                                       aCTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATAGTCAGTTCTACAGTCCAGTACTGGATACGCTGGaa                                                  >  2:6575359/1‑93 (MQ=255)

                                        cTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAAc                                                 >  1:2659214/1‑93 (MQ=255)

                                        cTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAAc                                                 <  2:389373/93‑1 (MQ=255)

                                        cTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAAc                                                 >  1:5349435/1‑93 (MQ=255)

                                        cTGGCTATCAACGCCATGAACTGCTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAAc                                                 <  2:5653281/93‑1 (MQ=255)

                                         tGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACa                                                >  1:3802682/1‑93 (MQ=255)

                                         tGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACa                                                <  1:5917441/93‑1 (MQ=255)

                                           gCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATg                                              <  2:5951919/93‑1 (MQ=255)

                                           gCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATg                                              <  2:6898336/93‑1 (MQ=255)

                                               tCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCAc                                          >  1:2201750/1‑93 (MQ=255)

                                               tCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCAc                                          <  1:4460328/93‑1 (MQ=255)

                                                     ccATGATCTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACa                                    <  1:5552988/93‑1 (MQ=255)

                                                     ccATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACa                                    <  2:6468936/93‑1 (MQ=255)

                                                      taTGATCTGGTAGGGCAGCTTTACTCCATGTTGTTCAACAATACTCATTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACAt                                   >  1:6998330/2‑93 (MQ=255)

                                                      cATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACAt                                   <  2:5963627/93‑1 (MQ=255)

                                                      cATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACAt                                   <  1:461074/93‑1 (MQ=255)

                                                          aaCTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTgg                               <  2:7494784/93‑1 (MQ=255)

                                                          aaCTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTgg                               >  1:1437827/1‑93 (MQ=255)

                                                           aCTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGgc                              >  1:1373965/1‑93 (MQ=255)

                                                           aCTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGgc                              <  1:2541460/93‑1 (MQ=255)

                                                           aCTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGgc                              <  1:2398279/93‑1 (MQ=255)

                                                            cTGGTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGCCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGgcc                             >  2:5115355/1‑92 (MQ=255)

                                                              ggTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCt                           >  2:3202520/1‑93 (MQ=255)

                                                              ggTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGGTCTACAGCCCAGTACTGGATACGCTGGACCATGGCACCGACCATGTGGCGCt                           >  2:5816787/1‑93 (MQ=255)

                                                              ggTAGGGCAGCCTTACTCCATGTTGTTAAACAATAGTCAGTTCTACAGTCCAGAACTGGATACGCGGGAACATGGCACGCAACATGTGgcgcg                           >  2:3826541/1‑92 (MQ=255)

                                                               gTAGGGCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTg                          >  2:2841535/1‑93 (MQ=255)

                                                                 aGGGCAGCCTTACTCCATGGTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACAGGGCACCGAACATGTGGCGGTGga                        >  1:5387586/1‑93 (MQ=255)

                                                                  gggCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGgag                       >  1:5100855/1‑93 (MQ=255)

                                                                   ggCAGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGgaga                      >  2:5156769/1‑93 (MQ=255)

                                                                      aGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCa                   <  1:6426270/93‑1 (MQ=255)

                                                                      aGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCa                   <  1:2367413/93‑1 (MQ=255)

                                                                      aGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCa                   <  1:493442/93‑1 (MQ=255)

                                                                      aGCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCa                   <  2:2575440/93‑1 (MQ=255)

                                                                       gCCTTACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAg                  <  2:4465873/93‑1 (MQ=255)

                                                                       gCCTTAATCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGTTCat                  >  2:7787306/1‑92 (MQ=255)

                                                                          ttACTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGttt               >  1:7449549/1‑93 (MQ=255)

                                                                            aCTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTc             >  2:4742877/1‑93 (MQ=255)

                                                                            aCTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTc             <  2:879647/93‑1 (MQ=255)

                                                                             cTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTcc            >  1:4370992/1‑93 (MQ=255)

                                                                             cTCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTcc            >  1:1203310/1‑93 (MQ=255)

                                                                              tCCCTGTTGTTCAACAATAATCTGTTCTACAGTACAGTACTGGATACGCTGGAACATGGCACCGTACATGTGGCGCTGGAGGTCAGTTTTCCa           <  1:7024486/93‑1 (MQ=255)

                                                                              tCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCa           <  2:6958426/93‑1 (MQ=255)

                                                                              tCCATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCa           >  1:109507/1‑93 (MQ=255)

                                                                               ccATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAg          >  1:5851041/1‑93 (MQ=255)

                                                                                cATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAgg         >  1:1099168/1‑93 (MQ=255)

                                                                                cATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAgg         >  2:5426581/1‑93 (MQ=255)

                                                                                cATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAgg         <  2:2315433/93‑1 (MQ=255)

                                                                                cATGTTGTTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAgg         >  1:5795435/1‑93 (MQ=255)

                                                                                  tgttgtTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTc       >  2:6559996/1‑93 (MQ=255)

                                                                                  tgttgtTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTc       <  2:2469523/93‑1 (MQ=255)

                                                                                  tgttgtTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTc       <  2:2188631/93‑1 (MQ=255)

                                                                                  tgttgtTCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTc       >  2:1473576/1‑93 (MQ=255)

                                                                                   gttgttCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGCACATGTGGCGCTGGAGATCAGTTTTCCAGGTCg      >  1:7636469/1‑93 (MQ=255)

                                                                                      gttCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCGGGTCGTGa   >  1:4880419/1‑93 (MQ=255)

                                                                                      gttCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCGTGa   <  2:2246459/93‑1 (MQ=255)

                                                                                       ttCAACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCGTGAc  >  2:7222517/1‑93 (MQ=255)

GGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCGTGAC  >  NC_000913/2318906‑2319085

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