BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,579,865	T→C	S252G (AGC→GGC)	hsdM ←	DNA methylase M

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,579,865	0	T	C	99.4%	491.7 / ‑4.8	162	S252G (AGC→GGC)	hsdM	DNA methylase M
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (81/80); total (82/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 = 3.82e-01

Only 100 of 176 total aligned reads displayed.

CAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGA  >  NC_000913/4579775‑4579953

cAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCAc                                                                                        <  1:1388259/93‑1 (MQ=255)

cAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCAc                                                                                        <  2:2391887/93‑1 (MQ=255)

aaaaGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCAc                                                                                        <  1:1537824/92‑1 (MQ=255)

 aaaGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCCTCGCCAcc                                                                                       >  2:2252024/1‑93 (MQ=255)

 aaaGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGGTTTCAACGTCGCCAcc                                                                                       >  1:275306/1‑93 (MQ=255)

  aaGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCAccc                                                                                      <  2:2707165/93‑1 (MQ=255)

  aaGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCAccc                                                                                      >  2:132434/1‑93 (MQ=255)

   aGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCa                                                                                     >  1:207653/1‑93 (MQ=255)

   aGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCa                                                                                     >  1:2672908/1‑93 (MQ=255)

        cGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGGTTTCACCGTCGCCACCCagagag                                                                                >  2:336660/1‑91 (MQ=255)

         gggTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAgtgt                                                                               <  2:1906417/93‑1 (MQ=255)

         gggTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAgtgt                                                                               >  1:438591/1‑93 (MQ=255)

         gggTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAgtgt                                                                               >  2:4379711/1‑93 (MQ=255)

         gggTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAgtgt                                                                               <  1:3669923/93‑1 (MQ=255)

         gggTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGGTTTCACCGTCGCCACCCCGAgtgt                                                                               >  1:4021312/1‑93 (MQ=255)

          ggTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGtt                                                                              <  2:467690/93‑1 (MQ=255)

          ggTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGtt                                                                              <  2:4025789/93‑1 (MQ=255)

          ggTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGtt                                                                              <  2:2950189/93‑1 (MQ=255)

          ggTAAGGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGGTTTCACCGTCGCCACCCAGAGTGtt                                                                              >  2:2971451/1‑93 (MQ=255)

            tAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGGTTTCACCGTCGCCCCCCAGAGTGTTGc                                                                            >  2:274372/1‑93 (MQ=255)

             aaTGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCCGAGGGTTGcc                                                                           >  2:1010654/1‑93 (MQ=255)

             aaTGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGcc                                                                           <  2:140304/93‑1 (MQ=255)

              aTGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGGTTTCACCGTCGCCACCCAGAGTGTTGccc                                                                          >  1:259182/1‑93 (MQ=255)

                 ttGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGGCCAGa                                                                       >  2:2344485/1‑93 (MQ=255)

                 ttGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGa                                                                       <  2:3249049/93‑1 (MQ=255)

                   ggTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACg                                                                     >  1:3160047/1‑93 (MQ=255)

                     tGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGa                                                                   <  1:4316613/93‑1 (MQ=255)

                     tGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGa                                                                   <  2:3638711/93‑1 (MQ=255)

                       ccTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGGTTTCACCGTCGCCACCCAGAGCGCTGCCCAGACGGAtt                                                                 >  2:3871520/1‑93 (MQ=255)

                         tGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTgc                                                               <  1:2693172/93‑1 (MQ=255)

                         tGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTgc                                                               >  2:3981154/1‑93 (MQ=255)

                         tGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTgc                                                               <  2:1707682/93‑1 (MQ=255)

                           cggcgCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTgcgc                                                             <  1:1441397/93‑1 (MQ=255)

                            ggcgCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGcgcc                                                            >  1:3964200/1‑93 (MQ=255)

                            ggcgCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGACCAGACGGAGTGcgcc                                                            >  1:1230647/1‑93 (MQ=255)

                              cgcTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGcgccgc                                                          >  2:73272/1‑93 (MQ=255)

                               gcTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCgccgcc                                                         <  1:4667743/93‑1 (MQ=255)

                                 tGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGt                                                       >  1:1174202/1‑93 (MQ=255)

                                  gCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTg                                                      <  2:163442/93‑1 (MQ=255)

                                   ccAAACGGCGTGTTAGTGGCTACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTgg                                                     <  2:4468452/93‑1 (MQ=255)

                                    cAAACGGCTGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGt                                                    <  1:921899/93‑1 (MQ=255)

                                    cAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGt                                                    <  2:2424186/93‑1 (MQ=255)

                                      aaCGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCg                                                  >  1:3467987/1‑93 (MQ=255)

                                       aCGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGa                                                 <  1:4254518/93‑1 (MQ=255)

                                        cggcggGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAg                                                <  1:3297657/93‑1 (MQ=255)

                                        cggcggGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAg                                                >  1:1123497/1‑93 (MQ=255)

                                        cggcggGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAg                                                <  2:788504/93‑1 (MQ=255)

                                          gcggGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGt                                              <  2:4077917/93‑1 (MQ=255)

                                          gcggGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGt                                              <  2:891002/93‑1 (MQ=255)

                                           cggGTTAGTGGCGACAATATGCGCCTTCGGCAGTTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGtt                                             <  2:4659318/93‑1 (MQ=255)

                                           cggGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGtt                                             <  2:252221/93‑1 (MQ=255)

                                              gTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGGTTGAGCCCCCGTGGTCGAGGTTGcc                                          >  1:665268/1‑93 (MQ=255)

                                               ttAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCt                                         <  1:4552824/93‑1 (MQ=255)

                                               ttAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCt                                         >  1:4144370/1‑93 (MQ=255)

                                               ttAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCt                                         >  2:1995983/1‑93 (MQ=255)

                                                 aGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTc                                       <  1:18146/93‑1 (MQ=255)

                                                     gCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAata                                   <  1:2080182/93‑1 (MQ=255)

                                                     gCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAata                                   <  1:4064947/93‑1 (MQ=255)

                                                     gCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAata                                   >  2:1869465/1‑93 (MQ=255)

                                                       gACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATc                                 <  1:843515/93‑1 (MQ=255)

                                                       gACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATc                                 <  2:4460786/93‑1 (MQ=255)

                                                         cAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGt                               >  1:4639921/1‑93 (MQ=255)

                                                         cAATATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGt                               <  1:624151/93‑1 (MQ=255)

                                                          aaTATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTg                              >  1:856019/1‑93 (MQ=255)

                                                          aaTATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTg                              >  1:3727118/1‑93 (MQ=255)

                                                          aaTATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTg                              >  1:3302960/1‑93 (MQ=255)

                                                          aaTATGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCCTg                              >  1:1049516/1‑93 (MQ=255)

                                                           atatGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTgc                             >  1:2821277/1‑93 (MQ=255)

                                                            tatGCGCTTTCGGTAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCTGTGGTCGAGGTTGCCTTTAATATCGTgca                            <  2:606220/93‑1 (MQ=255)

                                                            tatGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTgca                            >  1:1048476/1‑93 (MQ=255)

                                                            tatGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTgca                            <  2:3985608/93‑1 (MQ=255)

                                                            aatGCGCCTTCGGCCGGTTTCCACCGTCGCCACCCAGAGTGTTGCCCAGACGAATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTgca                            <  2:2555218/92‑1 (MQ=255)

                                                             atGCGCCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTgcag                           <  1:2660164/93‑1 (MQ=255)

                                                                cgcCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGcagcag                        >  2:3456456/1‑93 (MQ=255)

                                                                 gcCTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAgcagg                       >  1:1354371/1‑93 (MQ=255)

                                                                  ccTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAgcaggc                      <  2:3929131/93‑1 (MQ=255)

                                                                   cTTCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAgcaggca                     <  1:4010891/93‑1 (MQ=255)

                                                                    ttCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAgcaggcag                    <  2:1626378/93‑1 (MQ=255)

                                                                    ttCGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAgcaggcag                    >  1:2950908/1‑93 (MQ=255)

                                                                      cGGCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGtt                  >  1:4106070/1‑93 (MQ=255)

                                                                        gCAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTtca                >  2:3513601/1‑93 (MQ=255)

                                                                         cAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTtcat               <  2:1040014/93‑1 (MQ=255)

                                                                         cAGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTtcag               >  1:4501519/1‑92 (MQ=255)

                                                                          aGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTtcatc              >  2:2694368/1‑93 (MQ=255)

                                                                          aGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTtcatc              <  1:1587180/93‑1 (MQ=255)

                                                                          aGGTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTtcatc              >  1:2028971/1‑93 (MQ=255)

                                                                           ggTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCTTGGTCGAGGTTGCCTTCAATATCGTGGAGGAGGGAGGtcatca             >  2:4066618/1‑93 (MQ=255)

                                                                           ggTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTtcatca             <  1:3269933/93‑1 (MQ=255)

                                                                            gTTTTCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGCGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGATCATCAg            >  2:4177139/1‑93 (MQ=255)

                                                                             ttttCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGt           >  2:3591912/1‑93 (MQ=255)

                                                                             ttttCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGt           <  2:575092/93‑1 (MQ=255)

                                                                              tttCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTg          >  2:2940338/1‑93 (MQ=255)

                                                                                tCACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGcc        <  2:3152048/93‑1 (MQ=255)

                                                                                 cACCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCa       <  1:3658197/93‑1 (MQ=255)

                                                                                  aCCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAg      >  1:3011008/1‑93 (MQ=255)

                                                                                  aCCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAg      >  1:4475310/1‑93 (MQ=255)

                                                                                  aCCGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAg      >  2:3110732/1‑93 (MQ=255)

                                                                                   ccGTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAga     <  2:2672908/93‑1 (MQ=255)

                                                                                     gTCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAgacg   >  2:3191556/1‑93 (MQ=255)

                                                                                      tCGCCACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAgacga  >  1:2369798/1‑93 (MQ=255)

CAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGA  >  NC_000913/4579775‑4579953

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