BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,554,118	+C	coding (367/396 nt)	yhgKJ ←	RNA phosphage cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,554,117	1	.	C	98.5%	607.4 / 1.5	194	coding (368/396 nt)	yhgKJ	RNA phosphage cyclase
Reads supporting (aligned to +/- strand): ref base . (2/1); new base C (98/92); total (100/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 210 total aligned reads displayed.

TGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑C‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGCG  >  NC_000913/3554026‑3554200

ttgTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                                                     <  1:1886217/91‑1 (MQ=255)

tGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                                                     >  1:1184793/1‑93 (MQ=255)

        tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACt                                                                             >  1:3960722/1‑93 (MQ=255)

        tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACt                                                                             <  1:2313471/93‑1 (MQ=255)

        tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCCGCGGATAAAATC‑CGCAGGCGGCGGCGACGGGCTATCGG‑CGCCGCC‑GCTCACt                                                                             >  1:582137/1‑93 (MQ=255)

         gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACtt                                                                            >  1:1456706/1‑93 (MQ=255)

          atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTc                                                                           <  2:3543343/93‑1 (MQ=255)

          atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTc                                                                           <  1:5104819/93‑1 (MQ=255)

           tgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCa                                                                          >  2:6176830/1‑93 (MQ=255)

             aaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAac                                                                        <  2:1636573/93‑1 (MQ=255)

              aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAaca                                                                       <  2:2784970/93‑1 (MQ=255)

              aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAaca                                                                       >  2:1523296/1‑93 (MQ=255)

              aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCC‑GCTCACTTCAaca                                                                       >  2:5598258/1‑93 (MQ=255)

               ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAacac                                                                      >  2:5004053/1‑93 (MQ=255)

                 ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACAcgc                                                                    <  2:6510905/93‑1 (MQ=255)

                  cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACAcgcg                                                                   <  1:5884196/93‑1 (MQ=255)

                  cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACAcgcg                                                                   >  1:3142799/1‑93 (MQ=255)

                  cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACAcgcg                                                                   >  1:2171326/1‑93 (MQ=255)

                  cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACAcgcg                                                                   >  1:5036827/1‑93 (MQ=255)

                    aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGaa                                                                 <  1:6286942/93‑1 (MQ=255)

                    aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGaa                                                                 <  1:466495/93‑1 (MQ=255)

                    aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGAC‑GCTCACTTCATCACGCGaa                                                                 <  1:6567894/93‑1 (MQ=255)

                     ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAg                                                                >  2:5463300/1‑93 (MQ=255)

                      tttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAgg                                                               <  2:3650628/93‑1 (MQ=255)

                      tttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAgg                                                               <  1:6119632/93‑1 (MQ=255)

                       ttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGt                                                              <  2:1313045/93‑1 (MQ=255)

                          gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCa                                                           <  2:6130207/93‑1 (MQ=255)

                           ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCAt                                                          <  2:3794896/93‑1 (MQ=255)

                            cagcagCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATc                                                         >  1:4596178/1‑93 (MQ=255)

                              gcagcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCgg                                                       <  1:5627625/93‑1 (MQ=255)

                              gcagcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCgg                                                       >  1:5257897/1‑93 (MQ=255)

                              gcagcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCgg                                                       >  2:4557130/1‑93 (MQ=255)

                              gcagcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGG‑CGCCGCC‑GCTCACCTCAACACGCGAAGGTCCATCgg                                                       >  2:2572462/1‑93 (MQ=255)

                               cagcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGc                                                      >  2:4428300/1‑93 (MQ=255)

                               cagcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGc                                                      >  2:4190803/1‑93 (MQ=255)

                               cagcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGc                                                      <  2:5096698/93‑1 (MQ=255)

                               cagcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGc                                                      >  2:304793/1‑93 (MQ=255)

                                 gcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCaa                                                    >  1:3386405/1‑93 (MQ=255)

                                 gcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCaa                                                    >  1:4402191/1‑93 (MQ=255)

                                 gcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCaa                                                    >  1:41072/1‑93 (MQ=255)

                                    gCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                                                 <  2:6706/93‑1 (MQ=255)

                                    gCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                                                 >  2:3133858/1‑93 (MQ=255)

                                     cTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTAACGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                                                >  1:5389108/1‑93 (MQ=255)

                                      tCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac                                               >  1:5183524/1‑93 (MQ=255)

                                      cccAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac                                               <  2:2435251/92‑1 (MQ=255)

                                       ctaGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGTTCCATCGGCAAACcaca                                              <  1:2374253/91‑1 (MQ=255)

                                        cAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                                             <  2:2043519/93‑1 (MQ=255)

                                         aGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                                            <  1:1304649/93‑1 (MQ=255)

                                         aGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                                            <  2:2349757/93‑1 (MQ=255)

                                         aGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                                            >  1:6454007/1‑93 (MQ=255)

                                            aCCCGGCGGATAAAATC‑CGCAGGCGGCGCGGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                                         <  1:2863473/93‑1 (MQ=255)

                                            aCCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                                         >  2:4937950/1‑93 (MQ=255)

                                            aCCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                                         <  2:1535742/93‑1 (MQ=255)

                                            aCCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                                         >  2:3888818/1‑93 (MQ=255)

                                            aCCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                                         >  2:271959/1‑93 (MQ=255)

                                             cccGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg                                        >  1:696881/1‑93 (MQ=255)

                                              ccGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCggg                                       <  1:8711/93‑1 (MQ=255)

                                               cggcggATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGTGTCGC‑C‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc                                      >  2:4699600/1‑93 (MQ=255)

                                                ggcggAGAAAAGC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca                                     <  1:2423795/93‑1 (MQ=255)

                                                cgcggATAAAATCGAC‑AGGCGCGGCGGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca                                     <  1:4351551/92‑1 (MQ=255)

                                                   ggATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                                  >  1:3665701/1‑93 (MQ=255)

                                                   ggATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                                  <  1:5230433/93‑1 (MQ=255)

                                                    gagaaaaTC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                                 <  1:6210310/90‑1 (MQ=255)

                                                    gATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                                 >  2:3850952/1‑93 (MQ=255)

                                                    gATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                                 <  1:2233488/93‑1 (MQ=255)

                                                     aTAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc                                >  1:1244931/1‑93 (MQ=255)

                                                     aTAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc                                >  1:4033202/1‑93 (MQ=255)

                                                     aTAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc                                >  2:56423/1‑93 (MQ=255)

                                                       aaaaTC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt                              <  2:3991033/93‑1 (MQ=255)

                                                        aaaTC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt                             <  1:2136774/93‑1 (MQ=255)

                                                        aaaTC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt                             >  2:3927438/1‑93 (MQ=255)

                                                         aaTC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt                            >  1:1102627/1‑93 (MQ=255)

                                                         aaTC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt                            >  2:4803481/1‑93 (MQ=255)

                                                             ccgcaGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgca                         <  1:23313/89‑1 (MQ=255)

                                                             ccgcaGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgca                         <  2:679439/89‑1 (MQ=255)

                                                             ccgcaGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgca                         <  1:3812398/89‑1 (MQ=255)

                                                                 aGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAc                     <  2:5797641/93‑1 (MQ=255)

                                                                  ggcggcGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAcc                    >  2:3258403/1‑93 (MQ=255)

                                                                   gcggcgCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCa                   <  1:2654625/93‑1 (MQ=255)

                                                                   gcggcgCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCa                   <  2:642512/93‑1 (MQ=255)

                                                                   gcggcgCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCa                   <  2:1251128/93‑1 (MQ=255)

                                                                    cggcgCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCACCACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAg                  <  1:1292108/93‑1 (MQ=255)

                                                                     ggcgCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGc                 >  1:327290/1‑93 (MQ=255)

                                                                      gcgcCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCg                >  1:3603208/1‑93 (MQ=255)

                                                                       cgcCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGt               <  2:2216865/93‑1 (MQ=255)

                                                                        gcCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTa              <  1:4318487/93‑1 (MQ=255)

                                                                         ccGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTAc             <  2:381816/93‑1 (MQ=255)

                                                                         ccGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTAc             >  2:6044558/1‑93 (MQ=255)

                                                                           gACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACaa           <  2:6657471/93‑1 (MQ=255)

                                                                            aCGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAAc          >  2:5189389/1‑93 (MQ=255)

                                                                              gggTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACtt        <  2:4599959/93‑1 (MQ=255)

                                                                               ggTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTc       <  1:5339689/93‑1 (MQ=255)

                                                                                gTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc      >  2:3385033/1‑93 (MQ=255)

                                                                                gTTATCGGT‑GCAG‑CGGCGCACTTCAACACGAGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc      <  1:4212396/93‑1 (MQ=255)

                                                                                 ttATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCg     <  1:6391904/93‑1 (MQ=255)

                                                                                 ttATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCg     <  1:6517096/93‑1 (MQ=255)

                                                                                 ttATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCg     >  1:471455/1‑93 (MQ=255)

                                                                                   aTCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGgc   >  2:85302/1‑93 (MQ=255)

                                                                                   aTCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGgc   >  1:1378989/1‑93 (MQ=255)

                                                                                     tcggcgccgcc‑gcTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGATTGCAGCACCAGCGTACAACTTCCGgcg  >  1:6660448/5‑93 (MQ=255)

TGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑C‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGCG  >  NC_000913/3554026‑3554200

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