BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	553,102	C→T	intergenic (‑2/+116)	purE ← / ← lpxH	N(5)‑carboxyaminoimidazole ribonucleotide mutase/UDP‑2,3‑diacylglucosamine diphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	553,102	0	C	T	100.0%	1175.5 / NA	274	intergenic (‑2/+116)	purE/lpxH	N(5)‑carboxyaminoimidazole ribonucleotide mutase/UDP‑2,3‑diacylglucosamine diphosphatase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (131/143); total (131/143)

Only 100 of 278 total aligned reads displayed.

GGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATC  >  NC_000913/553004‑553199

  aTTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCt                                                                                               >  1:940756/1‑101 (MQ=255)

    ttCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTa                                                                                             >  2:940758/1‑101 (MQ=255)

      cGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTaaa                                                                                           <  1:940461/101‑1 (MQ=255)

       gAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTaaaa                                                                                          <  1:940243/101‑1 (MQ=255)

       gAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCCTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTaaaa                                                                                          <  2:940557/101‑1 (MQ=255)

        aaGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAAc                                                                                         >  1:940764/1‑101 (MQ=255)

         aGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACt                                                                                        >  2:940766/1‑101 (MQ=255)

              tCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCtgt                                                                                   >  2:940771/1‑101 (MQ=255)

                ggcggcGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACt                                                                                        <  1:940214/94‑1 (MQ=255)

                ggcggcGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTgtgg                                                                                 >  1:940774/1‑101 (MQ=255)

                  cggcgAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTgtggtg                                                                               <  1:939871/101‑1 (MQ=255)

                  cggcgAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTgtggtg                                                                               <  1:940341/101‑1 (MQ=255)

                     cgAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGcac                                                                            <  2:940213/101‑1 (MQ=255)

                      gAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGcaca                                                                           >  2:940782/1‑101 (MQ=255)

                      gAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGcaca                                                                           >  2:940779/1‑101 (MQ=255)

                       aaCTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACaa                                                                          >  1:940787/1‑101 (MQ=255)

                       aaCTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACaa                                                                          >  1:940784/1‑101 (MQ=255)

                        aCTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAAc                                                                         >  2:940788/1‑101 (MQ=255)

                          tGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAActc                                                                       <  1:940353/101‑1 (MQ=255)

                           gCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAActct                                                                      <  1:940333/101‑1 (MQ=255)

                           gCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAActct                                                                      >  2:940792/1‑101 (MQ=255)

                           gCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAActct                                                                      <  1:940309/101‑1 (MQ=255)

                            cATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACtctc                                                                     <  2:940484/101‑1 (MQ=255)

                             aTGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCg                                                                    >  1:940796/1‑101 (MQ=255)

                              tGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCgg                                                                   <  1:940320/101‑1 (MQ=255)

                              tGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCgg                                                                   <  2:940342/101‑1 (MQ=255)

                              tGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCgg                                                                   <  1:940367/101‑1 (MQ=255)

                                gTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCt                                                                 >  2:940800/1‑101 (MQ=255)

                                gTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCt                                                                 <  2:940365/101‑1 (MQ=255)

                                 tAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCtt                                                                >  1:940801/1‑101 (MQ=255)

                                 tAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCtt                                                                <  2:940445/101‑1 (MQ=255)

                                  aGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCttt                                                               <  2:940453/101‑1 (MQ=255)

                                   gCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTa                                                              <  1:940303/101‑1 (MQ=255)

                                     ccAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTaga                                                            <  1:940266/101‑1 (MQ=255)

                                      cAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTagag                                                           >  2:940805/1‑101 (MQ=255)

                                      cAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTagag                                                           <  2:940433/101‑1 (MQ=255)

                                      cAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTagag                                                           <  2:940385/101‑1 (MQ=255)

                                       aGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAAc                                                                         <  1:940282/86‑1 (MQ=255)

                                       aGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAgg                                                          >  2:940808/1‑101 (MQ=255)

                                       aGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAgg                                                          <  1:940542/101‑1 (MQ=255)

                                       aGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAgg                                                          >  1:940809/1‑101 (MQ=255)

                                          cGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGca                                                       >  2:940812/1‑101 (MQ=255)

                                           gCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGcac                                                      <  1:940203/101‑1 (MQ=255)

                                            cTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGcaca                                                     <  2:940504/101‑1 (MQ=255)

                                            cTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGcaca                                                     >  1:940816/1‑101 (MQ=255)

                                             ttttGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACag                                                    <  2:940399/101‑1 (MQ=255)

                                             ttttGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACag                                                    <  2:940240/101‑1 (MQ=255)

                                              tttGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACaga                                                   >  2:940818/1‑101 (MQ=255)

                                               ttGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACagag                                                  >  2:940822/1‑101 (MQ=255)

                                                tGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAgaga                                                 >  2:940825/1‑101 (MQ=255)

                                                 ggACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGaa                                                >  1:940827/1‑101 (MQ=255)

                                                  gACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAAt                                               >  2:940829/1‑101 (MQ=255)

                                                     cccATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGc                                            >  2:940834/1‑101 (MQ=255)

                                                     cccATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGc                                            <  2:940292/101‑1 (MQ=255)

                                                      ccATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCa                                           >  1:940836/1‑101 (MQ=255)

                                                        aTCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACg                                         <  2:940269/101‑1 (MQ=255)

                                                          cACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGa                                       <  2:940336/101‑1 (MQ=255)

                                                           aCGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGaa                                      >  1:940840/1‑101 (MQ=255)

                                                            cGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGaa                                      <  2:940155/100‑1 (MQ=255)

                                                            cGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGCGCACAGAGAATAGCACGGaaa                                     >  2:940841/1‑101 (MQ=255)

                                                              aTGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAaga                                   >  1:940844/1‑101 (MQ=255)

                                                              aTGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAaga                                   <  2:940526/101‑1 (MQ=255)

                                                               tGGCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAagag                                  >  2:940846/1‑101 (MQ=255)

                                                                ggCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAagaga                                 <  2:940476/101‑1 (MQ=255)

                                                                 gCGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAagagag                                >  2:940849/1‑101 (MQ=255)

                                                                  cGACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGNAAGAGAGc                               <  1:940441/101‑1 (MQ=255)

                                                                   gACACGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCa                              <  1:940499/101‑1 (MQ=255)

                                                                    acacGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCaa                             >  2:940853/1‑101 (MQ=255)

                                                                     cacGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAg                                             <  1:940409/84‑1 (MQ=255)

                                                                     cacGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAg                            <  1:940184/101‑1 (MQ=255)

                                                                     cacGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAg                            <  1:940391/101‑1 (MQ=255)

                                                                      acGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAgg                           >  2:940857/1‑101 (MQ=255)

                                                                      acGCGCCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAgg                           >  1:940856/1‑101 (MQ=255)

                                                                        gcgcCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGGGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAagag                                  <  1:940099/92‑1 (MQ=255)

                                                                         cgcCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGaaa                        <  1:940744/101‑1 (MQ=255)

                                                                          gcCGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGaaaa                       <  1:940407/101‑1 (MQ=255)

                                                                           ccGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAAc                      <  2:940497/101‑1 (MQ=255)

                                                                           ccGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAAc                      <  1:940574/101‑1 (MQ=255)

                                                                            cGGATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACg                     <  2:940308/101‑1 (MQ=255)

                                                                             ggATTATTGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGaaa                        <  1:940375/97‑1 (MQ=255)

                                                                               attattGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGtt                  >  2:940863/1‑101 (MQ=255)

                                                                               attattGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGtt                  >  1:940862/1‑101 (MQ=255)

                                                                                ttattGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTg                 >  1:940866/1‑101 (MQ=255)

                                                                                ttattGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTg                 <  2:940420/101‑1 (MQ=255)

                                                                                 tattGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGc                <  1:939953/101‑1 (MQ=255)

                                                                                 tattGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGc                <  2:940472/101‑1 (MQ=255)

                                                                                 tattGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGc                >  1:940869/1‑101 (MQ=255)

                                                                                  attGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCg               >  1:940874/1‑101 (MQ=255)

                                                                                  attGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCg               >  2:940872/1‑101 (MQ=255)

                                                                                   ttGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGt              >  1:940876/1‑101 (MQ=255)

                                                                                   ttGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGt              <  1:940480/101‑1 (MQ=255)

                                                                                    tGCGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCg               <  1:940729/99‑1 (MQ=255)

                                                                                      cGGGAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGc           >  1:940881/1‑101 (MQ=255)

                                                                                       gggAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCt          >  2:940883/1‑101 (MQ=255)

                                                                                        ggAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCtg         >  1:940886/1‑101 (MQ=255)

                                                                                         gAAGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCtgt        <  2:940516/101‑1 (MQ=255)

                                                                                           aGACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGa      >  2:940889/1‑101 (MQ=255)

                                                                                            gACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGaa     <  2:940614/101‑1 (MQ=255)

                                                                                            gACATGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGaa     >  2:940892/1‑101 (MQ=255)

                                                                                               aTGTGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATc  <  2:940322/101‑1 (MQ=255)

GGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATC  >  NC_000913/553004‑553199

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