BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	147,141	A→G	G58G (GGA→GGG)	yadD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	147,141	0	A	G	86.9%	775.2 / 97.2	312	G58G (GGA→GGG)	yadD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (18/23); new base G (138/133); total (156/156)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 320 total aligned reads displayed.

AGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGA  >  NC_000913/147050‑147232

  aTACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTTGGGGAGTTTCATTGAAGAGAGCCTGAAAGGgcac                                                                                          <  1:55100/93‑1 (MQ=255)

  aTACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGgcac                                                                                          <  2:1616625/93‑1 (MQ=255)

   tACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGcaca                                                                                         <  1:4779717/93‑1 (MQ=255)

   tACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGcaca                                                                                         >  1:1889467/1‑93 (MQ=255)

   tACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGcaca                                                                                         >  1:3614831/1‑93 (MQ=255)

    aCATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAg                                                                                        >  2:719944/1‑93 (MQ=255)

     cATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGc                                                                                       <  2:2222646/93‑1 (MQ=255)

      aTTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCa                                                                                      >  1:5800883/1‑93 (MQ=255)

      aTTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCa                                                                                      <  1:4386536/93‑1 (MQ=255)

       tttGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCAc                                                                                     >  1:3232523/1‑93 (MQ=255)

       tttGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCAc                                                                                     <  1:4270556/93‑1 (MQ=255)

        ttGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACg                                                                                    >  2:5810139/1‑93 (MQ=255)

           ccAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGAc                                                                                 <  2:3326412/93‑1 (MQ=255)

            cAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACg                                                                                <  1:6576943/93‑1 (MQ=255)

             aGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGt                                                                               >  2:6655028/1‑93 (MQ=255)

             aGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGt                                                                               >  2:685830/1‑93 (MQ=255)

              gTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTg                                                                              >  2:2677841/1‑93 (MQ=255)

               tGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGc                                                                             <  2:7312767/93‑1 (MQ=255)

                ggAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGtcgc                                                                            >  1:2426397/1‑90 (MQ=255)

                ggAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGct                                                                            >  1:2930051/1‑93 (MQ=255)

                 gAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGctc                                                                           <  2:2833586/93‑1 (MQ=255)

                  aaTTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGctct                                                                          <  1:1644922/93‑1 (MQ=255)

                    ttACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTAt                                                                        <  1:5463942/93‑1 (MQ=255)

                      aCGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTc                                                                      <  1:6966629/93‑1 (MQ=255)

                       cgcgAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTcc                                                                     <  2:2090861/93‑1 (MQ=255)

                          gAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTg                                                                  >  2:4252181/1‑93 (MQ=255)

                          gAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTg                                                                  >  1:7468843/1‑93 (MQ=255)

                            aCTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCa                                                                >  2:2466155/1‑93 (MQ=255)

                            aCTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCa                                                                >  2:606522/1‑93 (MQ=255)

                                tgtgACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATg                                                            <  2:7425156/93‑1 (MQ=255)

                                tgtgACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGGGCTCTATTCCGTGCAAATg                                                            >  1:7296254/1‑93 (MQ=255)

                                 gtgACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGc                                                           <  1:5274209/93‑1 (MQ=255)

                                 gtgACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGc                                                           >  2:1867053/1‑93 (MQ=255)

                                   gACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAg                                                         <  1:624010/93‑1 (MQ=255)

                                   gACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAg                                                         <  1:1354348/93‑1 (MQ=255)

                                    aCCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAgg                                                        <  1:4651989/93‑1 (MQ=255)

                                    aCCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAgg                                                        <  1:3021798/93‑1 (MQ=255)

                                        cAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCaa                                                    <  1:5766207/93‑1 (MQ=255)

                                         aaCACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAAt                                                   >  1:988404/1‑93 (MQ=255)

                                         aaCACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAAt                                                   <  1:1483136/93‑1 (MQ=255)

                                           cacGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATcc                                                 <  2:437619/93‑1 (MQ=255)

                                             cGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCg                                               <  1:5398771/93‑1 (MQ=255)

                                              gCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCgg                                              >  2:4748116/1‑93 (MQ=255)

                                              gCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCgg                                              <  1:2968808/93‑1 (MQ=255)

                                              gCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCgg                                              >  2:2550249/1‑93 (MQ=255)

                                               cTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGt                                             >  2:5141144/1‑93 (MQ=255)

                                               cTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGt                                             <  2:5843445/93‑1 (MQ=255)

                                                ttCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGtt                                            >  1:4828918/1‑93 (MQ=255)

                                                ttCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGtt                                            >  1:758460/1‑93 (MQ=255)

                                                 tCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                                           <  1:2695752/93‑1 (MQ=255)

                                                  cATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTAt                                          >  2:733640/1‑93 (MQ=255)

                                                   aTTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATc                                         >  2:5658087/1‑93 (MQ=255)

                                                   aTTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATc                                         >  2:4584212/1‑93 (MQ=255)

                                                   aTTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATc                                         >  2:3062975/1‑93 (MQ=255)

                                                   aTTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATc                                         >  1:1817872/1‑93 (MQ=255)

                                                    tttAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCt                                        >  1:5936257/1‑93 (MQ=255)

                                                      tAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGc                                      <  2:4119917/93‑1 (MQ=255)

                                                       agagTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCa                                     >  2:5552258/1‑93 (MQ=255)

                                                        gagTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAt                                    <  2:5986094/93‑1 (MQ=255)

                                                          gTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgt                                  >  1:1610551/1‑93 (MQ=255)

                                                           tCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgtt                                 <  1:3666470/93‑1 (MQ=255)

                                                           tCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgtt                                 <  2:2260330/93‑1 (MQ=255)

                                                           tCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgtt                                 <  1:1060066/93‑1 (MQ=255)

                                                           tCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgtt                                 <  2:5471961/93‑1 (MQ=255)

                                                           tCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgtt                                 <  1:5655944/93‑1 (MQ=255)

                                                             ggggAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgttgt                               >  2:6924038/1‑93 (MQ=255)

                                                              gggAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg                              >  2:3835718/1‑93 (MQ=255)

                                                              gggAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg                              >  1:1972997/1‑93 (MQ=255)

                                                               ggAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGa                             >  2:5233865/1‑93 (MQ=255)

                                                                gAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGAt                            >  2:5058707/1‑93 (MQ=255)

                                                                gAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGAt                            >  2:2408977/1‑93 (MQ=255)

                                                                 aGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGAtt                           <  1:3829437/93‑1 (MQ=255)

                                                                 aGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGAtt                           <  1:515968/93‑1 (MQ=255)

                                                                 aGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGAtt                           >  1:2171748/1‑93 (MQ=255)

                                                                  gTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTg                          >  2:3005517/1‑93 (MQ=255)

                                                                  gTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTg                          <  2:1335861/93‑1 (MQ=255)

                                                                    ttCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGaa                        <  1:830656/93‑1 (MQ=255)

                                                                     tCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAac                       <  2:3683186/93‑1 (MQ=255)

                                                                       aTTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAacac                     >  2:6586105/1‑93 (MQ=255)

                                                                        ttGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACAcc                    <  1:6675207/93‑1 (MQ=255)

                                                                          gAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCaa                  >  1:6333766/1‑93 (MQ=255)

                                                                           aaGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCaaa                 <  1:2834516/93‑1 (MQ=255)

                                                                              agagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagca              >  1:5195262/1‑93 (MQ=255)

                                                                               gagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaa             <  2:150498/93‑1 (MQ=255)

                                                                                 gCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaagc           <  2:899047/93‑1 (MQ=255)

                                                                                 gCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaagc           >  1:3127025/1‑93 (MQ=255)

                                                                                 gCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaagc           <  1:6833396/93‑1 (MQ=255)

                                                                                  ccTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGcc          <  2:6151244/93‑1 (MQ=255)

                                                                                  ccTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGcc          >  2:3411596/1‑93 (MQ=255)

                                                                                  ccTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGcc          >  1:2482018/1‑93 (MQ=255)

                                                                                   cTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCg         <  1:2092230/93‑1 (MQ=255)

                                                                                     gAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGa       <  1:2560664/93‑1 (MQ=255)

                                                                                     gAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGa       >  2:3601465/1‑93 (MQ=255)

                                                                                     gAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGa       <  2:7004013/93‑1 (MQ=255)

                                                                                      aaaGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGAt      >  2:1179736/1‑93 (MQ=255)

                                                                                        aGGGCACGGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATaa    <  1:4113881/93‑1 (MQ=255)

                                                                                        aGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATaa    <  2:4693397/93‑1 (MQ=255)

                                                                                        aGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATaa    <  1:3358735/93‑1 (MQ=255)

                                                                                         ggACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAg   <  2:2141892/93‑1 (MQ=255)

                                                                                          gACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGa  <  2:4838790/93‑1 (MQ=255)

AGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGA  >  NC_000913/147050‑147232

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